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"Optoelectronic  lll-V  Heterostructures  by  Gas-Source  MBE" 

I.  Research  Objective 

The  objective  of  this  research  program  was  to  grow  by  molecular  beam  epitaxy  (MBE) 
quantum  well  heterostructures  of  the  alloy  InGaAsP  and  related  materials  for  use  in  high 
performance  devices.  The  first  task  was  to  grow  InGaAsP/GaAs  and  InGaAsP/inP  structures  for 
shaliow  quantum  well  (QW)  optical  modulators  and  other  optoelectronic  devices.  The  second  task 
was  to  develop  selective-area  regrowth  techniques  for  lateral  definition  of  complex  optoelectronic 
integrated  devices.  The  method  of  growth  was  gas-source  MBE,  a  technique  which  we  have 
previously  shown  to  produce  high  quality  QW  heterostructures  containing  alternating  layers  of  lll-V 
arsenide  and  phosphides. 


II.  Major  Accomplishments 

Under  this  contract,  we  have  concentrated  on  the  development  of  several  lll-V  material 
combinations  for  QW  heterostructures  in  order  to  optimize  optoelectronic  device  performance. 
Our  major  accomplishments  for  each  material  combination  were  as  follows: 

InGaP/lnAlP  and  InGaP/InGaP  Quantum  Wells: 

*  Combined  InGaP  as  the  well  material  with  InAiP  as  the  barrier  material  and 
reported  the  first  measurement  of  the  band  offset  energies  in  this  heterojunction 
system.  InGaP/lnAlP  and  InGaP/InGaAlP  are  grown  on  GaAs  substrates  and 
used  in  visible  emitters,  including  edge  emitting  lasers  (red)  and  LEDs  (red,  yellow, 
green). 

*  First  report  of  L-band  recombination  in  InGaP/lnAlP  QWs.  Demonstrated  that 
in  unstrained  narrow  wells  or  in  tensile-strained  wide  wells,  the  previously 
neglected  L-band  can  dominate  recombination  and  thus  adversely  affect  laser 
operation. 


*  Extended  our  work,  both  theoretical  and  experimental,  on  the  InAIP/InGaP 
system  to  include  the  introduction  of  controlled  amounts  of  strain  as  a  means  to 
further  tailor  the  heterostructure  to  meet  device  requirements. 

*  First  results  for  optical  modulators  using  lnxGa,.xP/lnyGa,.yP  (x^y)  MQWs.  In 
preliminary  devices,  transmission  modulators  on  GaP  substrates  exhibited  15% 
modulation  at  10V  and  594  nm,  and  asymmetric  reflection  modulators  on  GaAs 
substrates  showed  27%  modulation  at  20V  at  605  nm,  indicating  suitability  of 
InGaP  QW  MBE  modulators  for  plastic  optical  fiber  data  transmission  systems. 

InGaAs/InGaP  Quantum  Wells: 

*  Reported  the  first  strain-compensated  QW  heterostructures  of  the  material 
combination  InGaAs/InGaP,  and  obtained  fewer  defects  than  conventional 
InGaAs/GaAs  modulators  operating  at  a  wavelength  of  1  pm. 

*  Using  Stark-effect  red-shift  modulators,  obtained  42%  modulation  at  5V,  and  by 
using  a  superiattice,  fiirsf  observation  of  blue-shift  modulation  in  the  InGaAs/InGaP 
system. 

*  Careful  measurements  showed  that  the  switching  speed  was  slower  and  the 
modulation  saturates  at  lower  intensity  levels  in  InGaAs/InGaP  modulators  than  in 
conventional  AIGaAs/GaAs  structures,  primarily  because  of  the  larger  band  offset 
energies. 


InGaAsP/InGaAsP  Quantum  Wells: 

*  Demonstrated  that  GSMBE  can  produce  atomically  abrupt  QW  interfaces  in  the 
InGaAsP  system;  thus  the  randomness  associated  with  composition  does  not  limit 
the  abruptness  of  heterointerfaces  employing  quaternary  alloys. 
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Gas-Source  MBE  Technology 


*  Published  a  review  of  the  gas-source  MBE  technology  we  developed  under 

AFOSR  sponsorship  during  the  last  7  years. 

Most  of  the  work  performed  under  this  contract  addressed  the  firet  task  of  our  objective 
by  obtaining  a  better  understanding  of  the  material  properties  of  InGaP,  InAlP,  InGaAs,  and 
InGaAsP  heterostructures,  and  then  used  this  improved  understanding  to  design  optical 
modulators  for  operation  at  given  wavelength.  The  details  of  our  work  can  be  found  in  the 
reprints  of  the  published  papers  included  in  this  report. 

For  the  second  task  as  originally  proposed,  we  fabricated  lasers  at  1.3  pm  wavelength  with 
the  lowest  yet  reported  threshold  current  for  material  grown  by  gas-source  MBE.  The  lasers  used 
InAsP/InGaAsP  QWs,  a  material  developed  under  a  separate  AF/DARPA  contract  (Hanscom 
Laboratory).  The  buried  heterostructure  growth  process  needed  for  the  high  performance  lasers 
was  to  be  developed  under  this  contract,  but  we  found  the  selectivity  of  the  GSMBE  growth 
process  for  the  InGaAsP  material  was  insufficient  to  produce  the  desired  device  structures. 
Instead,  we  carried  out  the  growth  of  the  critical  InAsP/InGaAsP  active  region  by  GSMBE,  and 
then  had  the  buried  heterostructure  regrowth  done  by  MOCVD,  through  a  collaboration  at  Bell 
Laboratories.  The  resulting  wafers  were  fabricated  into  lasers  in  our  laboratories,  and  the  device 
performance  has  been  outstanding. 


III.  Publications  Under  AFOSR  Sponsorship 
(copies  of  publications  are  enclosed) 

1.  O.  Buccafusca,  J.A.L.  Chilla,  C.S.  Menoni,  J.J.  Rocca,  M.J.  Hafich,  L.M.  Woods,  and  G.Y.  Robinson, 
"Non-resonant  Tunneling  in  InGaP/lnAlP  Asymmetric  Double  Quantum  Wells",  Appl.  Phys.  Lett.  62, 399  (25 
Jan  1993). 

2.  J.W.  Kim,  C.W.  Chen,  T.J.  Vogt,  L.M.  Woods,  G.Y.Robinson,  and  D.L.  Lite,  "Strained  Layer  InGaAs/GaAs 
and  InGaAs/InGaP  Multiple  Quantum  Well  Optical  Modulators  Grown  by  Gas-Source  MBE,"  IEEE  Photonics 
Technol.  Letters  5  (9),  987  (September  1993). 

3.  C.W.  Chen,  J.W.  Kim,  P.  Silvestre,  M.J.  Hafich,  G.Y.  Robinson,  and  D.L.  Lite,  "A  Low  Drive  Voltage 
Electroabsorption  Modulator  Using  an  InGaAs/lnP  Superlattice,"  J.  Appl.  Phys.  74,  5895  (1  November 
1993). 

4.  D.  Patel,  M.J.  Hafich,  G.Y.  Robinson,  and  C.S.  Menoni,  "Direct  Determination  of  the  Band  Discontinuities 


3 


in  InGaP/lnAlP  Multiple  Quantum  Wells",  Phys.  Rev.  B  48  (24),  18,031  (15  December  1993). 

5.  LM.  Woods,  P.Silvestre,  P.  Thiagarajan,  G.A.  Patrizi,  and  G.Y.  Robinson,  "Photoluminescence  and 
Interface  Abruptness  in  InGaAsP/InGaAsP  Quantum  Wells",  J.  Electron.  Materials  23  (1 1)  1229,  (November 
1994). 

6.  M.E.  Watson,  J.LA.Chilla,  J.J.  Rocca,  J-W  Kim,  D.L.  Lite,  T.J.  Vogt,  and  G.Y.  Robinson,  "Saturation 
Intensity  and  time  response  of  InGaAs/InGaP  MQW  modulators"  IEEE  J.  Quantum  Electron.  31  (2),  254-260 
(February  1995). 

7.  G.Y.  Robinson,  "Gas-Source  Molecular  Beam  Epitaxy",  in  the  Handbook  of  Thin  Film  Process 
Technology:  A  Guide  to  Current  Methods,  D.A.  Glocker,  S.l.  Shah,  eds..  Institute  of  Physics,  Bristol  and 
Philadelphia,  1995,  pp  A.2.2.1-A.2.2.22. 

8.  D.  Patel,  K.  Interhozingher,  P.  Thiagarajan,  G.Y.  Robinson,  and  C.S.  Menoni,  "L-Band  Recombination 
in  InGaP/lnAlP  Multiple  Quantum  Wells",  Physical  Rev.  B,  53  (19),  12633-12636  (15  May  1996). 

9.  J.W.  Kim,  Y.J.  Lee,  T.J.  Vogt,  G.A.  Patrizi,  G.Y.  Robinson,  and  D.L.  Lite,  "Wannier-Stark  localization  in 
a  strained  InGaAs/InGaP  superlattice",  J.  Appl.  Phy.  79  (9),  7161-7163  (1  May  1996). 

10.  T.J.  Vogt,  P. Thiagarajan,  and  G.  Y.  Robinson,  "MBE  Growth  of  InGaP  MQW  Structures  on  GaP  for 
Optical  Modulators",  accepted  for  publication  in  J.  Vac.  Sc.  Technol.,  (May/June  1996). 

11.  C.M.  Menoni,  O.F.  Buccafusca,  M.C.  Marconi,  D.  Patel,  J.J  Rocca,  G.Y.  Robinson,  and  S.M.  Goodnick, 
"Effects  of  Indirect  G-L  and  G-X  Transfer  on  the  Carrier  Dynamics  of  InGaP/lnAlP  Multiple  Quantum  Wells", 
submitted  to  Appl.  Phy.  Letts.,  1996. 


IV.  Conference  Presentations  and  Invited  Talks  Under  AFOSR  Sponsorship 

1.  K.  Mahalingam,  N.  Otsuka,  M.J.  Hafich,  and  G.Y.  Robinson,  "Correlation  of  Spontaneous  Ordering  in 
InGaAlP  Alloys  with  Surface  Reconstruction",  APS  Meeting,  Seattle,  March  1993. 

2.  J.W.  Kim,  C.W.  Chen,  T.J.  Vogt,  LM.  Woods,  G.Y.Robinson,  and  D.L.  Lite,  "Strained  Layer  InGaAs/GaAs 
and  InGaAs/InGaP  Multiple  Quantum  Well  Optical  Modulators  Grown  by  Gas-Source  MBE,"  Fifth 
International  Conference  on  Indium  Phosphide  and  Related  Materials,  Paris,  France,  April  1993. 

3.  M.  Prasad,  O.E.  Martinez,  C.S.  Menoni,  J.J.  Rocca,  J.L.A.  Chilla,  M.J.  Hafich,  and  G.Y.  Robinson, 
"Transient  Grating  Measurements  of  Ambipolar  Diffusion  and  Carrier  Recombination  in  InGaP/lnAlP  MQWs 
and  InGaP  Bulk  Material",  Electronic  Materials  Conference,  Santa  Barbara,  June  1993. 

4.  T.J.  Vogt,  J.  K.  Kim,  L.M.  Woods,  G.A.  Patrizi,  D.L.  Lite,  and  G.Y.  Robinson,  "Strain-Compensated 
InGaP/InGaAs  Multiple  Quantum  Well  Modulators  Grown  by  Gas-Source  MBE",  Electronic  Materials 
Conference,  Santa  Barbara,  June  1993. 

5.  C.S.  Menoni,  D.  Patel,  M.J.  Hafich,  and  G.Y.  Robinson,  "Band  Offsets  in  InGaP/lnAlP  Quantum  Wells 
Using  High  Pressure",  presented  at  High  Pressure  Science  and  Technology  Conference,  Colorado  Springs, 
June  1993. 

6.  M.E.  Watson,  J.L.A.  Chilla,  J.J.  Rocca,  J-W  Kim,  D.L.  Lite,  T.J.  Vogt,  and  G.Y.  Robinson. 
"Modulation-Saturation  Measurements  on  InGaAs/InGaP  MQW  Modulators",  IEEE  LEOS  Conference  on 
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Lasers  and  Electro-optics,  Anaheim,  May  1994. 

7.  M.E.  Watson,  J.L.A.  Chilla,  J.J.  Rocca,  J-W  Kim,  D.L  Lite,  T.J.  Vogt,  and  G.Y.  Robinson,  "Saturation 
Intensity  and  Time  Response  of  InGaAs/InGaP  MOW  Modulators",  Electronic  Materials  Conference, 
Boulder,  June  1 994. 

8.  O.  Buccafusca,  J.E.  Fouquet,  G.A.  Patrizi,  LM.  Woods,  C.S.  Menoni,  J.J.  Rocca,  and  G.Y.  Robinson, 
"Optical  Properties  of  Lattice-Matched  InGaP/lnAlP  Multiple  Quantum  Wells",  Electronic  Materials 
Conference,  Boulder,  June  1994. 

9.  L.M.  Woods,  P.Silvestre,  P.  Thiagarajan,  G.A.  Patrizi,  G.Y.  Robinson,  K.M.  Jones,  and  M.  Al-Jassim, 
"Interface  Abruptness  in  InGaAsP/InGaAsP  Quantum  Wells  as  Determined  by  Photoluminescence",  North 
American  Conference  on  Molecular  Beam  Epitaxy,  Urbana,  IL,  October  1994. 

10.  G.  Y.  Robinson,  "Band  Offset  Energies  in  InGaP/lnAlP  Quantum  Wells"  (invited  talk)  Center  for 
Compound  Semiconductor  Microelectronics,  University  of  Illinois,  Champaign-Urbana,  November  1994. 

11.  K.  Interholzinger,  D.  Patel,  C.S.  Menoni,  O.  Buccafusca,  P.  Thiagarajan,  T.J.  Vogt,  L.M.  Woods,  G.Y. 
Robinson,  and  J.E.  Fouquet,  "Band  structure  and  band  alignments  of  strained  and  lattice-matched 
InGaP/lnAlP  heterostructures".  Electronic  Materials  Conference,  Charlottesville,  Virginia,  June  1995. 

12. 0.  Buccafusca,  M.C.  Marconi,  D.  Patel,  C.S.  Menoni,  M.  Prasad,  J.J.  Rocca,  and  G.Y.  Robinson,  "Effect 
of  G-X  coupling  on  carrier  lifetime  of  InGaP/lnAlP  multiple  quantum  wells".  Electronic  Materials  Conference, 
Charlottesville,  Virginia,  June  1995. 

13.  T.J.  Vogt,  P.  Thiagarajan,  and  G.Y.  Robinson,  "Visible  (X,=575nm)  Multiple  Quantum  Well  Modulators 
on  GaP  Substrates  by  Gas-Source  MBE",  North  American  Conference  on  Molecular  Beam  Epitaxy,  College 
Park,  MD,  September  1995. 

14.  D.  Patel,  K.  Interhozingher,  P.  Thiagarajan,  G.Y.  Robinson,  and  C.S.  Menoni,  "Direct  Identification  of 
the  L-Band  in  InGaP/lnAlP  Multiple  Quantum  Wells",  APS  Meeting,  St.  Louis,  MO,  March  1996. 

15.  D.  Patel,  K.  Interhozingher,  P.  Thiagarajan,  G.Y.  Robinson,  and  C.S.  Menoni,  "L-Band  Recombination 
in  ln^Gai.xP/lnAIP  Multiple  Quantum  Wells",  Electronic  Materials  Conference,  Santa  Barbara,  CA,  June 
1996. 

16.  D.  Patel,  K.  Interhozingher,  P.  Thiagarajan,  G.Y.  Robinson,  and  C.S.  Menoni,  "L-Band  Recombination 
in  In^Gai.^P/lnAlP  Multiple  Quantum  Wells",  High  Pressure  Physics  Conference,  Schwabisch-Gmung, 
Germany,  July  1996. 


V.  Awards 

G.  Y.  Robinson,  elected  Fellow  of  IEEE,  January  1996,  for  "contributions  to  the  field  of 
microelectronic  materials". 


VI.  Technology  Transfer 
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The  first  technology  transfer  under  this  contract  was  to  Charles  Cemy  and  Fritz 
Schuermeyer  at  Wright  Labs  and  consisted  of  custom  designed  epitaxial  InGaP  wafers  for  p- 
channel  FET  fabrication.  The  FETs  were  to  be  evaluated  for  a  use  in  a  CMOS  GaAs-based 
technology.  At  present,  the  wafers  are  being  processed. 

The  second  technology  transfer  was  to  Hewlett-Packard  Laboratories  where  Julie  Fouquet 
evaluated  our  InAIP/InGaP  QW  heterostructures  by  PLE.  From  our  combined  efforts  to  measure 
the  band  offsets,  HP  wanted  to  confirm  how  far  the  wavelength  could  be  decreased  in  InGaAlP 
LEDs  and  still  obtain  sufficient  light  output.  We  are  now  examining  the  use  of  strain  to  further 
decrease  the  wavelength.  Hewlett-Packard  is  the  only  American  manufacturer  of  commercial 
LEDs. 

VII.  Personnel 

Gary  Y.  Robinson,  Professor,  Principal  Investigator 

Gary  A.  Patrizi,  Research  Associate,  MSEE  from  Colorado  State  University,  US  citizen  (now  at 
Philips  Semiconductor,  Albuquerque,  NM). 

Prabhu  Thiagarajan,  Graduate  Student  Research  Assistant  and  Research  Associate,  Electrical 
Engineering,  non-US  citizen  with  green  card. 

Larry  M.  Woods,  Graduate  Student  Research  Assistant,  US  citizen,  received  MS  in  Electrical 
Engineering,  December  1994  (now  at  Golden  Photonics,  Golden  CO). 

Tim  J.  Vogt,  Graduate  Student  Research  Assistant,  Electrical  Engineering,  PhD  candidate,  US 
citizen  (to  be  employed  at  Honeywell,  Minneapolis,  MN). 

Joel  Fastenau,  Graduate  Student  Research  Assistant,  Electrical  Engineering,  PhD  candidate,  US 
citizen. 

Daryl  Pulver,  Student  Technician,  US  citizen. 
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VIII.  Research  Results 


The  majority  of  the  results  obtained  under  this  contract  have  been  published  or  have  been 
accepted  for  publication  in  the  near  future.  These  results  were  summarized  above  and  the  details 
can  be  found  in  previous  annual  contract  reports  and  in  the  reprints  of  our  published  papers 
provided  in  the  next  section.  Here  we  will  highlight  our  most  recent  research  results  which  have 
yet  to  be  published. 

We  have  fabricated  reflection  modulators  operating  at  near  600  nm,  the  wavelength  for 
low  loss  transmission  in  plastic  optical  fibers,  by  growing  InGaP  QWs  on  "GaAsP  substrates". 
The  "GaAsP  substrates"  are  commercially  available  GaAs  wafers  with  thick  epitaxial  layer  (about 
100  pm)  of  GaAS(j  sPo  4.  We  then  grow  a  thick  buffer  layer  of  Ino  jGao  yP  which  lattice  matches  to 
the  "GaAsP  substrate"  and  avoids  the  generation  of  misfit  dislocations  obtained  with  growth 
directly  on  GaP  and  GaAs  substrates.  A  MQW  region  is  then  grown  on  the  buffer,  with  tensile- 
strained  well  layers  of  lno4Gao6P  and  compressively  strained  barrier  layers  of  Ino  jGaogP-  VVith 
the  proper  choice  of  well  and  barrier  thickness,  operation  at  the  desired  wavelength  can  be 
obtained  with  a  MQW  stack  of  zero  net  strain.  This  combination  of  materials  has  proved  to  be 
the  most  successful  approach  we  have  yet  explored  in  attempting  to  obtain  optical  modulation  at 
short  wavelengths  with  the  InGaP  system. 

As  shown  in  the  figure,  the  device  consists  of  a  pin  diode,  with  the  InGaP/InGaP  MQW 
stack  as  the  i-region  and  p*  and  n*  InGaP 
layers  used  for  ohmic  contacts.  To  increase 
the  contrast  ratio  of  the  reflection  modulator, 
a  Fabry-Perot  optical  cavity  is  formed 
between  two  mirrors,  a  distributed  Bragg 
reflector  (DBR)  below  the  pin  structure  and 
the  air-semiconductor  interface  above.  The 
DBR  consists  of  19.5  pairs  of  X/4  layers  of 
InojGaojP  and  InojAlojP.  The  reflectivity 
spectra  of  the  two  mirrors  are  chosen  to 
match  the  change  in  MQW  absorption  with 
bias  resulting  from  the  quantum  confined 
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stark  effect.  The  device  is  an  asymmetrical  Fabry-Perot  reflection  modulator. 

The  measured  reflection  spectrum  is  shown  in  the  figure.  At  a  wavelength  of  603  nm,  the 
maximum  change  in  reflectivity  is  27%  and  the  contrast  ratio  is  2.0  for  a  change  in  applied  voltage 
of  20  V.  This  is  the  largest  modulation  yet  achieved  with  the  InGaP  materials  system,  and  is 
state-of-the-art  performance  for 
Stark-effect  modulators  in  the 


570-650  nm  range.  With 
improvements  in  the  mirrors,  we 
expect  to  achieve  substantially 
larger  contrast  ratios.  Such  a 
device  could  then  be  used  in 
short-haul  plastic  fiber  data 
transmission  systems  by  serving 
as  both  a  photodetector  for 
incoming  signals  and  a  modulator 
(by  using  some  of  the 
downstream  light)  for  outgoing 
signals.  Such  two-way 
transmission  networks  using  only 
one  optical  source  and  many 
optical  modulators  operating  in 
reflection  mode  would  minimize 
the  use  of  low  reliability  and 
expensive  components  in  robust 
networks. 


Reflectance  Spectra  (mbe1074,  Ga  rich  top) 


wavelength(nm) 
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The  band  structure  of  Ino.53Gao.47P/Ino.50Alo.soP  multiple  quantum  wells  grown  by  molecular- 

p-ssure-dependent-photownescenfe  measlT^t 

were  ^  t  A  otoluminescence  signals  from  the  direct-gap  well  and  the  inditect  barrier 

Tfr  “  pressure  up  to  4  GPa.  High  pressure  transformed  the  multiple 

bvThrT7f  “  7^^  I*'"  1  staggered  aligned,  type  II  at  1.1  GPa.  This  transition  was  evidenced 

L  mL  7  A  *  Pl^^tol^^i^escence  signal  due  to  the  recombination  of  carriers  separated 

in  momentum  and  space.  The  simultaneous  detection  of  this  transition  and  that  of  the  barrier 
material,  aUowed  the  direct  determination  of  a  valence-band  ofi^t  energy  of  (0.24±0.05)eV  without 
any  information  on  parameters  of  the  bulk  materials.  Considering  that  the  total  band-gap 
djscontinuijr  for  this  heterostructure  system  is  0.50  eV  at  20  K,  an  approximate  band-gap  spUtting 

the  pressmfcrefficW  InxGai_*P/In,AU.„P  interface.  Variations  in 

aj^lnm^nt  transitions  in  the  barrier  indicated  that  the  valence  band 

ght-hole  states  with  biaxial  strain  induced  in  the  epilayers  by  applying  pressure. 


I.  INTRODUCTION 

The  possibility  of  using  quantum  confinement  to  vary 
the  effective  band  gap  of  In.Gai_.P/In^Ali_„P  het¬ 
erostructures  m  the  wavelength  range  540-650  nm  makes 
this  system  suitable  for  the  engineering  of  optoelec- 
tromc  devices  which  operate  in  the  visible  region  of 
the  opticd  spectrum.  In,Gai_,P  heterostructures  have 
been  used  as  the  wide-band-gap  component  of  edge- 
emitting  and  vertical-cavity-surface-emitting  semicon¬ 
ductor  l^ers  2  and  have  also  been  employed  in  the  de- 
sign  of  heterojunction  bipolar  transistors.^  In  spite  of 
the  broad  application  of  this  material,  basic  knowledge  of 
bmdamental  parameters  such  as  the  band  discontinuities 
winch  me  required  for  the  design  of  efficient  devices  are 
still  unknown. 

We  report  the  first,  to  our  knowledge,  direct  de¬ 
termination  of  the  band  discontinuities  in  In*Gai_xP/ 
[nj,Ali_j,P  multiple  quantum  wells  (MQW’s)  from  low 
;emperature  photoluminescence  (PL)  measurements  per¬ 
formed  at  high  pressure.  The  valence  band  offset  of  this 
leterostructure  was  directly  deduced  from  the  PL  signals 
rorresponding  to  different  transitions  within  the  well  and 
Darrier  obtained  by  excitation  with  three  different  lines 
of  an  argon  ion  laiser. 

High  pressure  photoluminescence  measurements  offer 
;he  possibility  of  determining  the  band-gap  discontinu- 
ties  in  multiple  quantum  well  structures  which  show  a 
;ype-I  alignment  for  the  direct  gap  and  a  type-II  align- 
nent  for  the  ini^ect  X  minima.-*  This  method  uses  high 
pressure  to  modify  the  band  structure  of  the  heterostruc- 
•ure  materials  to  obtain  a  band  alignment  in  which  the 
owest  conduction  band  state  and  the  highest  valence 
land  state  are  separated  in  space.  Under  this  condition 

■163-1829/93/48(24)/18031(6)/$06.00 


indirect  recombination  of  the  spatially  separated  carriers 
is  possible.  Monitoring  the  variation  with  pressure  of  this 
indirect  transition  in  addition  to  the  recombination  of 
carriers  in  the  barrier  material  allows  the  determination 
of  the  energy  of  the  barrier  X  minima  at  atmospheric 
pressure  with  respect  to  the  well  and  barrier  valence 
b^d,  respectively.  The  valence  band  offset  can  then  be 
directly  calculated  as  the  energy  difference  between  these 
two  transitions.  The  mmn  advantage  of  the  high  pressure 
techmque  is  that  it  does  not  require  any  sissumptions  con¬ 
cerning  the  bulk  parameters  nor  does  it  require  theoret¬ 
ical  modeling  or  a  special  configuration  of  the  sample  as 
needed  in  some  of  the  other  techniques  which  have  been 
used  to  determine  the  band  offsets.®”*®  Using  this  tech¬ 
nique  the  band  alignment  of  several  heterostructure  sys¬ 
tems  such  as  GaAs/Al*Gai_,As,*  GaAs/In,Gai_,,P,*i 
and  more  recently  GaAs/Ga*Asi_^P  (Ref.  12)  have  been 
determined. 

The  rest  of  this  paper  is  orgeinized  cis  follows:  in  Sec.  II 
the  sample  characteristics  and  the  high  pressure  photo- 
luminescence  experiments  are  described,  the  results  are 
presented  and  discussed  in  Sec.  Ill  and  Sec.  IV  summa¬ 
rizes  the  main  findings  of  this  work. 


II.  EXPERIMENTAL  DETAILS 

The  Ino.53Gao.47P /Ino.soAlo.soP  MQW^  samples  were 
grown  by  gas-source  moleculan-beam  epitaxy  on  a  semi- 
insulating  (100)  GaAs  substrate.  The  MQW’s  were  com¬ 
posed  of  30  periods  of  8.5  nm  wells  and  22.5  nm  barriers. 
An  In-Al-P  capping  layer  of  210  nm  was  also  deposited. 
The  samples  were  unintentionally  n-doped  with  a  back¬ 
ground  concentration  of  10*®  cm”®.  Transmission  elec- 
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tron  microscopy,  and  photoluminescence  measurements 
indicated  that  the  Ina;Gai_a;P  material  is  mostly  disor¬ 
dered.  A  detailed  description  of  the  growth  and  charac¬ 
terization  of  these  samples  is  given  in  Ref.  13. 

For  the  high  pressure  photoluminescence  measure¬ 
ments,  the  GaAs  substrate  was  removed  by  chemical 
etching  and  a  small  portion  (100  x  150  /xm)  of  the  epilayer 
was  placed  in  the  gasketed  diamond  anvil  cell  (DAC).^^ 
A  small  piece  of  ruby  was  positioned  next  to  the  sam¬ 
ple  for  the  measurement  of  pressure,  and  a  second  ruby 
chip  placed  in  the  DAC  body  outside  the  pressure  cham¬ 
ber  was  used  as  the  zero  pressure  reference. A  sam¬ 
ple  of  bulk  Ina;Gai_a,P  was  also  placed  outside  the  pres¬ 
sure  chamber  for  the  mecisurement  of  the  band  gap  of 
this  material.  Argon  was  used  as  the  pressure  transmit¬ 
ting  medium.  The  choice  of  this  pressure  medium  along 
with  the  reduced  thickness  of  the  sample  were  required 
to  guarantee  quasihydrostatic  conditions  at  low  temper¬ 
ature  and  at  the  highest  pressure  of  4  GPa  achieved  in 
these  experiments.  The  DAC  was  cryogenically  cooled 
to  20  K  with  a  close-cycle  He  refrigerator.  Using  a  spe¬ 
cially  designed  cryostat  chamber  the  pressure  in  the  cell 
was  changed  while  still  at  low  temperature.  Our  abil¬ 
ity  to  vary  the  pressure  at  low  temperature  avoided  any 
hysteresis  problems  which  may  be  encountered  in  systems 
requiring  room  temperature  pressurization  subsequent  to 
cooling  and  measuring. 

The  PL  from  the  multiple  quantum  wells  was  excited 
using  three  different  lines  from  an  Ar"*"  laser,  514,  488, 
and  458  nm,  each  with  a  constant  power  of  3.5  mW.  This 
variation  in  the  wavelength  allowed  the  excitation  of  dif¬ 
ferent  transitions  in  the  well  and  barrier  materials.  The 
PL  was  dispersed  using  a  0.75  m  spectrometer  equipped 
with  a  1200  g/mm  grating  and  detected  by  a  thermoelec- 
trically  cooled  photomultiplier  tube.  Standard  chopping 
techniques  and  a  lock-in  amplifier  were  used  to  analyze 
the  signal. 

III.  RESULTS  AND  DISCUSSION 

Typical  PL  spectra  at  different  pressures  obtained  by 
excitation  with  the  different  Ar"^  lines  are  shown  in  Fig. 

1.  Curves  a,  6,  and  c  were  obtained  with  an  excitation 
energy  of  2.41  eV  (514  nm).  At  low  pressures  the  main 
PL  feature,  labeled  jE7i,  corresponds  to  recombination  in 
the  well  involving  the  n  =  1  lowest  confined  conduction 
band  state  and  the  Ihh  heavy  hole  valence  band  state. 
For  pressures  higher  than  1.1  GPa  a  second  much  weaker 
and  broader  peak  (£*2)  appears  at  an  energy  of  2.09  eV. 
For  pressures  higher  than  2.3  GPa  another  transition  of 
characteristics  similar  to  those  of  the  E2  peak  is  detected 
at  an  energy  of  2.2  eV.  This  peak  was  labeled  £3.  Both 
£2  and  £3  were  resolvable  even  at  the  highest  pressure 
of  4  GPa.  Excitation  of  the  carriers  with  2.54  eV  (488 
nm)  photons  resulted  in  PL  spectra  (not  shown)  with  the 
same  features  as  those  shown  in  curves  a-c. 

When  the  PL  was  excited  with  2,7  eV  (458  nm)  pho¬ 
tons  an  additional  broad  peak  at  2.35  eV  is  also  observed. 
This  peak,  which  was  labeled  £4,  is  shown  in  Fig.  1, 
curve  d,  along  with  the  £1  peak  and  a  much  weaker  peak 


FIG.  1.  Typical  PL  spectra  of  InxCai-zP/InyAli^yP  mul¬ 
tiple  quantum  wells  at  different  pressures  at  20  K.  Curves  a-c 
were  obtained  with  excitation  of  the  carriers  with  a  photon 
energy  of  2.41  eV  (514  nm)  and  curve  d  was  obtained  with  an 
excitation  of  2.7  eV  (458  nm)  photons. 


at  2.27  eV  which  corresponds  to  the  £3  transition.  For 
the  highest  excitation  energy  of  2.7  eV,  it  is  possible  to 
excite  carriers  in  the  well  with  sufficient  excess  energy  to 
scatter  to  the  indirect  X  level  as  well  as  to  excite  carriers 
within  the  barrier.  The  unique  identification  of  these  PL 
transitions,  as  shown  below,  requires  the  knowledge  of 
the  corresponding  pressure  coefficients. 

A  weak  but  resolvable  peak  at  1.94  eV  is  also  observed 
in  the  PL  spectra  near  zero  pressure,  for  all  three  different 
excitation  energies.  Olsthoorn  et  al}^  have  speculated  on 
the  nature  of  this  transition  in  In^Gai-xP  alloys  as  an 
isolated  center  whose  PL  was  both  temperature  and  exci¬ 
tation  independent.  The  fact  that  this  transition  is  only 
observable  at  very  low  pressures  suggests  that  this  tran¬ 
sition  is  possible  only  when  the  n  =  1  level  is  resonant 
with  the  center.  Regardless,  the  origin  of  this  PL  peak 
remains  speculative. 

The  measured  energy  variation  with  pressure  of  the 
£1  —  £4  transitions  is  shown  in  Fig.  2.  The  energy 
of  the  £1  transition  shows  the  characteristic  direct  gap 
behavior,  increasing  with  pressure  at  a  rate  of  (92=b3) 
meV/GPa.  The  remaining  transitions,  on  the  other 
hand,  show  negative  pressure  coefficients,  a  behavior 
which  is  characteristic  of  an  X-like  band.  The  high¬ 
est  energy  transition  £4  is  identified  as  that  correspond¬ 
ing  to  recombination  between  electrons  in  the  barrier  X 
level  (Xt)  and  holes  in  the  barrier  valence  band.  The 
atmospheric  pressure  value  of  the  energy  of  these  min¬ 
ima  agrees  very  well  with  previous  determination  of  the 
band  structure  of  bulk  luyAli-yP.^®  At  about  1.1  GPa 
the  onset  of  the  first  direct- to-indirect  crossover  (£1— £2) 
occurs.  The  transition  £2  is  identified  cls  that  involving 
states  from  the  Xi  level  in  the  barrier  and  the  Ihh  va¬ 
lence  band  states  in  the  well.  This  PL  is  the  result  of  the 
recombination  of  carriers  which  are  separated  in  momen- 
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FIG.  2.  Variation  with  pressure  of  the  measured  PL  tran¬ 
sition  energies  from  Ino.53Gao.47P/Ino. 5 Aio.sP  multiple  quan¬ 
tum  wells.  Results  of  different  excitations  are  indicated  with 
different  symbols:  square,  458  nm;  circle,  488  nm;  and  tri¬ 
angle,  514.5  nm.  Four  energy  bands  are  clearly  identified. 
The  lines  correspond  to  the  least  squares  fit  of  the  data  from 
which  the  pressure  coefficient  and  the  energy  at  atmospheric 
pressure  of  each  of  the  transitions  were  obtained. 


tmn  and  space,  that  is,  at  this  pressure  the  type-I  MQW 
is  transformed  into  a  type-II  staggered  aligned  system.  A 
second  direct-to-indirect  crossover  {Ei  —  Es)  is  observed 
at  2.3  GPa.  This  crossover,  which  is  accompanied  by  the 
quenching  of  the  Ei  PL  signal,  is  aissigned  to  the  cross¬ 
ing  of  the  n  =  1  level  and  the  well  X  minima  The 

corresponding  crossover  in  bulk  InxGai_a.P  at  low  tem¬ 
peratures  occurs  at  2.6  GPa.^"^  The  lower  transition  pres¬ 
sure  of  this  direct-to-indirect  crossover  in  the  MQW  is  a 
consequence  of  the  larger  pressxure  coefficient  of  the  direct 
gap  of  Ina;Gai_a;P  in  a  confined  structure  when  compcired 
with  bulk  material.  Thus  the  trcinsition  E3  is  assigned 
to  recombination  of  carriers  from  the  minima  into  the 

well  valence  band.  Figure  3  schematically  shows  the  dif¬ 
ferent  PL  transitions  observed  in  InarGai._*P/InyAli_yP 
multiple  quantum  wells  in  these  experiments. 

Corroborative  evidence  of  the  Ei  —  E2  direct-to- 
indirect  crossover  was  obtained  by  monitoring  the  PL 
intensity  as  a  function  of  pressure.  A  sharp  reduction 
in  the  Ei  PL  peak  intensity  is  observed  at  1.1  GPa,  as 
shown  in  Fig.  4.  This  behavior  in  the  PL  intensity  is  typ¬ 
ical  of  a  direct-to-indirect  crossover  and  has  been  previ¬ 
ously  observed  in  other  III-V  materials.^®’^®  Also  shown 
in  this  figure  is  the  PL  peak  intensity  corresponding  to 
the  E2  transition,  which  is  approximately  constant  in  the 
pressure  range  of  these  mecisurements. 

Another  feature  which  establishes  our  assignment  of 
the  E2  transition  was  obtciined  by  monitoring  the  depen¬ 
dence  of  the  PL  peak  intensity  of  the  indirect  transitions 
on  excitation  power.  Figure  5  shows  the  excitation  power 
dependent  PL  at  2.4  and  3.12  GPa  for  the  two  indirect 
transitions  E2  and  £3.  It  is  clear  that  the  PL  peak  in¬ 
tensity  from  the  indirect  tramsition  E^  within  the  well 
is  constant  with  excitation  power;  however,  the  indirect 


FIG.  3.  Schematic  diagram  of  the  band  structure  of 
lna:Gai_*P/InyAli_yP  MQW’s  showing  the  PL  transitions 
observed  at  20  K  at  different  pressures.  The  valence  band  off¬ 
set  can  be  directly  determined  from  E2  and  E4  with  a  small 
correction  due  to  the  heavy  hole  band  confinement  energy. 


transition  E2  from  the  bcirrier  into  the  well  is  excitation 
power  dependent.  The  observed  blueshift  in  the  PL  peak 
intensity  of  this  transition  is  attributed  to  band  bend¬ 
ing  due  to  an  electric  field  created  at  the  interface  by 
the  spatially  separated  canders  and  has  been  previously 
observed  in  the  GaAs/AlAs  quantum  well  system^®  and 
in  GaAs/In*Gai_a.P  MQW’s  at  high  pressure.^^  The  ab¬ 
sence  of  any  redshift  in  the  peak  intensity  of  the  Ez  tran¬ 
sition  indicates  negligible  lattice  heating. 

The  pressure  coefficients  of  the  different  transitions 
were  obtained  from  the  slope  of  the  least  square  fit  to  the 
data.  The  y  intercept  of  these  fits  yielded  the  energy  of 
the  different  electronic  levels  in  the  well  and  barrier  at  at- 


FIG.  4.  Pressure  dependence  of  the  PL  peak  intensity  from 
the  main  direct  transition  Ei  and  that  corresponding  to  the 
indirect  transition  E2  between  the  barrier  and  the  well.  A 
two-orders-of-magnitude  decrease  in  the  peak  intensity  is  ob¬ 
served  at  about  1  GPa  due  to  the  direct-to-indirect  transitions 
between  the  barrier  and  the  well.  Lines  through  the  data  are 
only  a  guide  to  the  eye. 
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FIG.  5.  Excitation  power  dependent  PL  at  2.4  and  3.1 
GPa.  The  E2  PL  pecJc  energy  from  the  barrier  is  blueshifted 
with  excitation  power  whereas  the  indirect  transition  within 
the  well  {E3)  is  relatively  constant  with  power.  The  blueshift 
is  a  consequence  of  band  bending  due  to  the  electric  field  cre¬ 
ated  by  the  spatially  separated  carriers. 

mospheric  pressure.  These  parameters  and  the  crossover 
pressures  of  the  two  direct-to-indirect  transitions  have 
been  summarized  in  Table  I. 

Relevant  information  on  the  variation  with  pressure  of 
the  valence  band  offset  can  be  obtained  by  comparing  the 
pressure  coefficients  of  the  E2  and  E4,  transitions.  Since 
both  transitions  originate  at  the  same  level  X5,  their  pres¬ 
sure  coefficient  should  be  the  same.  However,  our  results 
show  that  the  E2  transition  decreases  with  pressure  at 
about  twice  the  rate  of  the  E?4  transition.  This  difference 
in  the  pressure  coefficients  indicates  that  either  the  Ihh 
well  valence  band  state  or  the  barrier  valence  band  are 
changing  with  pressure,  possibly  due  to  the  presence  of 
biaxial  strain. 

In  two-dimensional  systems  the  effect  of  hydrostatic 
pressure  results  in  the  generation  of  a  biaxial  strain  in 
the  growth  plane  due  to  the  difference  in  the  elastic  con¬ 
stants  of  the  well  and  barrier  materials.^^  This  biaxial 
strain  is  responsible  for  the  shifting  of  the  heavy  and 
light  hole  bands  and  the  splitting  of  the  indirect  conduc¬ 
tion  band  minima. We  calculated  the  pressure  induced 
strain  for  our  sample  using  elastic  theory  and  the  corre¬ 
sponding  shift  of  the  heavy  and  light  hole  bands.^"^  The 
elastic  constants  of  both  ternary  materials  were  interpo¬ 
lated  from  the  binary  constituents^^  and  were  assumed 
to  be  constant  in  the  pressure  range  of  these  experiments. 
These  calculations  show  that  as  pressure  is  increased  a 
biaxial  tensile  strain  is  induced  in  the  In^Gai-xP  well 
and  a  compressive  strain,  less  than  half  the  value  of  that 
in  the  well,  is  generated  in  the  luyAli^yF  barrier.  This 


FIG.  6.  Energy  shift  of  the  heavy  and  light  hole  bands  in 
InzGai_iP  as  a  function  of  pressure  The  shifts  are  calculated 
with  respect  to  the  average  valence  band  energy  using  the 
solid-model  approach  (Ref.  24).  Both  bands  were  assumed  to 
be  separated  by  16  meV  at  atmospheric  pressure. 

difference  in  the  magnitude  of  the  strain  in  the  weU  and 
barrier  indicates  that  the  changes  in  the  band  alignment 
are  mostly  due  to  the  variation  with  pressure  of  the  well 
valence  band  states.  The  corresponding  shift  in  energy 
with  pressure  of  the  heavy  and  light  hole  well  valence 
band  states  is  shown  in  Fig.  6.  The  heavy  hole  band 
shifts  with  pressure  at  a  rate  of  —13.5  meV/GPa  while 
the  light  hole  energy  increases  with  pressure  at  a  rate  of 
-1-21.5  meV/GPa.  Taking  into  account  a  separation  of  16 
meV  between  the  heavy  and  light  hole  band  as  calculated 
from  the  envelope-approximation,^®  it  is  found  that  the 
two  bands  cross  at  0.5  GPa.  Therefore,  in  almost  the 
entire  pressure  range  of  these  experiments  the  light  hole 
band  is  the  top  valence  band  state.  Since  this  band  has 
a  positive  pressure  coefficient  and  the  indirect  transition 
E2  has  a  negative  pressure  coefficient,  then  the  mecisured 
rate  is  larger  than  that  of  the  E4  transition.  A  pressure 
rate  for  the  light  hole  band  of  18  meV/GPa  is  obtained 
from  the  difference  between  the  measured  pressure  rates 
of  the  E2  and  E4  transitions,  in  good  agreement  with 
the  rate  obtained  in  the  calculations.  The  shifting  of  the 
light  hole  band  is  also  responsible  for  the  increased  pres¬ 
sure  coefficient  of  the  E3  transition.  The  lower  crossover 
pressure  of  this  transition  when  compared  to  bulk^^  is 
due  to  the  splitting  of  the  Xb  minima  with  strain.  In 
calculating  the  splitting  of  the  well  valence  band  states 
with  strain  we  have  ignored  the  variation  in  the  confine¬ 
ment  energies  due  to  the  change  of  the  effective  masses 


TABLE  I.  Zero  pressure  energy  band  position,  pressure  coefficient,  and  crossover  pressure  of  the 
different  transitions  observed  in  In*Gai_*P/InyAii-yP  multiple  quantum  weUs. 


E  (eV) 

[at  P  =  0  GPa] 

dEgjdP 
.  (meV/GPa) 

Transition  pressure 
(GPa) 

In-Ga-P  {El) 

1.984±0.004 

92±3 

In-Ga-P/In-Al-P  (£2) 

2.12db0.01 

-28±5 

1.1±0.2  {El  -  Ei) 

In-Ga-P  (£3) 

2.272±0.007 

-31±3 

2.3±0.5  {El  -  Ez) 

In-Al-P  {Ei) 

2.35±0.01 

-10±4 

14 
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with  pressure  and  that  due  to  the  change  in  the  well 
width  with  pressure.  These  corrections  are  very  small, 
have  opposite  pressure  coefficients,  and  therefore  their 
effect  on  the  confinement  energy  change  with  pressure 
is  negligible.  Prom  our  experimental  results  and  calcu¬ 
lations  we  find  that  in  the  InxGai-arP/InyAli^yP  sys¬ 
tem  the  valence  band  alignment  cannot  be  assumed  to 
remain  constant  with  pressure  as  has  been  the  case  in 
the  GaAs/AlxCai-xAs  system."^  While  this  variation  in 
energy  of  the  valence  band  states  modifies  the  pressure 
coefficients  of  the  PL  transitions,  it  does  not  alter  the  en¬ 
ergy  difference  between  the  conduction  and  valence  band 
states  at  atmospheric  pressure. 

The  determination  of  the  valence  band  offset  of  the 
Lia,Gai_a;P/InyAli_yP  heterointerface  is  straightforward 
knowing  the  energy  of  the  different  indirect  transitions  at 
atmospheric  pressure.  In  reference  to  Fig.  3,  the  valence 
band  offset  is  given  by 


AEy  =  E4  --  E2  +  (1) 

where  Eihh  is  the  confinement  energy  of  the  heavy  hole 
band.  In  the  envelope  approximation,  E^ihh  is  found  to  be 
equal  to  9  meV  when  an  In^Gai^^P  hole  effective  mass 
of  0.46mo  is  used.^^  With  this  value  of  Ei^h  a  valence 
band  offset  of  0.24  eV  is  obtained  from  Eq.  (1).  This 
value  of  the  band  offset  is  smaller  than  that  predicted  by 
the  model  solid  approach  of  Van  de  Walle.^^ 

The  calculation  of  the  valence  band  offset  using  Eq. 
(1)  assumes  that  the  acceptor  concentration  is  negligible 
in  both  the  well  and  barrier  materials,  since  our  samples 
are  n-type.  However,  in  the  presence  of  acceptor  levels  in 
either  the  well  or  barrier  material,  this  expression  needs 
to  be  modified  to  account  for  the  binding  energy  of  the 
acceptor  level.  We  have  calculated  the  valence  band  off¬ 
set  in  this  case  considering  two  limiting  situations,  one 
in  which  an  acceptor  level  with  a  binding  energy  of  50 
meV  is  located  in  the  well  material,  and  the  other  in 
which  the  same  acceptor  level  is  located  in  the  barrier. 
In  the  first  case,  the  valence  band  offset  is  foimd  to  be 
0.19  eV  while  in  the  second  Ccise  a  value  of  0.29  eV  for 
the  valence  band  offset  is  obtained.  Including  the  accep¬ 
tor  binding  energy  and  the  experimental  uncertainties 
we  find  AEv  =  (0.24  ±  0.05)  eV.  To  our  knowledge,  the 
only  other  reported  measurement  of  the  band  offset  ener¬ 
gies  for  the  In-Ga-P/In-Al-P  interface  has  been  made  by 
Watanabe  and  Ohba  using  capacitance-voltage  {C  —  V) 
carrier  profiling. They  obtained  a  value  at  300  K  of 
0.11  eV  for  AEc  for  an  Ino.sGao.sP/Ino.sAlo.sP  hetero¬ 
junction  grown  by  metal-organic  chemical  vapor  epitaxy. 
Using  a  band-gap  difference  of  0.50  eV,  their  result  indi¬ 
cates  AE^  of  0.39  eV,  considerably  higher  than  our  value 
of  AEv  =(0.24±0.05)  eV.  The  reason  for  this  discrepancy 
may  be  the  precision  of  the  measurement  techniques.  The 
AEc  obtained  by  the  C  —  V  method  is  now  known  to 
be  strongly  influenced  by  the  structure  (i.e.,  interfacial 
charge,  doping,  and  compositional  nonuniformities)  and 
measurement  conditions  (i.e.,  frequency,  temperature), 
and  thus  large  errors  in  AEc  can  easily  occur In  fact, 
Watanabe  and  Ohba  note  that  a  simulated  carrier  pro¬ 
file  based  on  their  AEc  value  did  not  fit  their  measured 


profile.^^  In  comparison,  the  pressure-dependent  PL  tech¬ 
nique  employed  here  measures  energy  levels  directly  and 
does  not  require  additional  knowledge  of  the  bulk  or  in¬ 
terface  properties.  Thus,  the  value  of  AE^  =  0.24±0.05 
eV  is  free  of  the  errors  inherent  in  the  C  —  V  technique. 

The  conduction  band  offset  can  be  calculated  from  the 
measured  value  of  the  valence  band  offset  and  the  band- 
gap  energy  difference  of  the  well  and  barrier  materials. 
The  value  of  the  direct  band-gap  of  bulk  In^Gai-xP  was 
determined  from  the  maximum  of  the  photoluminescence 
signal  from  the  bulk  sample  corrected  for  the  exciton 
binding  energy.  A  value  of  1.95  eV  was  obtained  at  20  K. 
The  band-gap  of  the  barrier  materials  was  taken  equal  to 
2.45  eV  from  previous  measurements  on  the  band  struc¬ 
ture  of  In-Al-P.^®  The  total  band  discontinuity  for  this 
heterostructure  system  is  calculated  to  be  0.50  eV  at  20 
K.  Therefore,  the  conduction  band  offset  energy  is  0.26 
eV  and  a  52:48  band-gap  splitting  is  determined  to  be 
the  band  lineup  in  the  InxGai«a:P/l%Ali.-3;P  MQW.  A 
schematic  of  the  band  structure  diagram  of  this  MQW  is 
shown  in  Fig.  7.  The  only  assumption  made  in  this  band 
diagram  involves  the  energy  of  the  barrier  T  minimum; 
however,  this  assumption  does  not  play  a  role  in  deter¬ 
mining  the  valence  band  offset  as  discussed  above.  The 
band  diagram  also  shows  the  position  of  the  indirect  X 
levels  in  the  well  and  barrier  materials  with  respect  to  the 
well  and  barrier  valence  bands  as  obtained  from  the  mea¬ 
surements.  The  136  meV  separation  between  the  barrier 
X  level  and  the  lowest  confined  conduction  band  state  in 
the  well  sets  a  lower  limit  for  the  well  width  of  16  A  for 
a  type-I  structure.  This  result  agrees  well  with  the  mea¬ 
surements  of  Hafich  et  who  observed  the  quench¬ 

ing  of  the  photoluminescence  in  Ina.Gai_.a.P/I%Ali-j;P 
for  well  widths  lower  than  18  A,  Therefore  the  shortest 
wavelength  that  can  be  achieved  with  this  heterostruc¬ 
ture  system  is  563  nm.  Shorter  wavelengths  would  result 
in  a  type-II  MQW  which  has  been  shown  to  have  low 
radiative  efficiency  and  radiative  lifetimes  up  to  several 
These  characteristics  along  with  the  low  mobility 
of  the  carriers  associated  with  the  X  valley  are  detrimen¬ 
tal  in  device  performance. 
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FIG.  7.  Energy  band  diagram  of  InajGai-xP/InyAli-^/P  at 
atmospheric  pressme  and  20  K  .  The  values  for  the  band-gap 
energies  and  energies  of  the  barrier  and  well  X  minima  are 
all  obtained  from  the  experiments. 
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IV.  CONCLUSIONS 

The  band  structure  of  Ino.5Gao.5P/I1io.5Alo.5P  niuh 
tiple  quantum  wells  was  determined  from  low  temper¬ 
ature  photoluminescence  measurements  at  high  pres¬ 
sure.  Selective  excitation  of  different  transitions  in  the 
well  and  barrier  materials  allowed,  to  our  knowledge, 
the  first  direct  determination  of  the  valence  band  off¬ 
set  at  atmospheric  pressure.  This  heterostructure  sys¬ 
tem  transformed  from  a  type-I  alignment  to  a  type-II 
staggered-aligned  system  at  1.1  GPa.  This  direct-to- 
indirect  crossover  was  accompanied  by  a  sharp  decrease 
in  the  intensity  of  the  direct  gap  photoluminescence  sig¬ 
nal.  The  valence  band  alignment  was  foxmd  to  vary  with 
press  me  due  to  pressure  induced  strains  in  the  epilayers. 


ACKNOWLEDGMENTS 


The  authors  kindly  acknowledge  Professor  M.  Chan- 
drcisekar  and  B.  Weinstein  for  stimulating  discussions. 
Professor  J.  Mahan  for  lending  us  the  close-cycle  He  re¬ 
frigerator,  and  Professor  D.  Lile  for  the  photomultiplier 
tube  and  cooled  housing.  This  work  was  supported  by 
the  Center  for  Optoelectronic  Computer  Systems  spon¬ 
sored  by  the  National  Science  Foundation/Engineering 
Research  Center  Grant  No.  ECD  9015128,  the  Colorado 
Advanced  Material  Institute,  the  State  of  Colorado,  and 
the  Air  Force  OflSce  of  Scientific  Research  (Contract  No. 
F49620-93-1-0021). 


^  G.  Hatakoshi,  K.  Itaya,  M.  Ishikawa,  M.  Okajima,  and  Y. 

Uematsu,  IEEE  J.  Quantum  Electron  QE-27,  1476  (1991). 
^  R.  P.  Schneider,  Jr.,  R.  P.  Bryan,  J.  A.  Lott,  G.  R.  Olbright, 
Appl.  Phys.  Lett.  60,  1830  (1992). 

^  M.  J.  Mondry  and  H.  Kroemer,  IEEE  Electron.  Devices 
Lett.  EDL-6,  175  (1985). 

^  D.  J.  Wolford,  T.  F.  Kuech,  and  J.  A.  Bradley,  J.  Vac.  Sci. 
Technol.  B  4,  1043  (1986). 

^R.  C.  Miller,  A.  C.  Gossard,  D.  A.  Kleinman,  and  O, 
Munteamu,  Phys.  Rev.  B  29,  1329  (1984). 

®  W.  1.  Wang,  Solid  State  Electron.  29,  133  (1986). 

^  M.  A.  Haase,  N.  Pan,  and  G.  E.  Stillman,  Appl.  Phys.  Lett. 
54,  1457  (1989). 

®  J.  Batey  and  S.  L.  Wright,  J.  Appl.  Phys.  59,  200  (1986). 
®  M.  A.  Rao,  E.  J.  Caine,  H.  Kroemer,  S.  I.  Long,  and  D.  1. 
Babic,  J.  Appl.  Phys.  61,  643  (1987). 

T.  J.  Drummond  and  I.  J.  Fritz,  Appl,  Phys,  Lett.  47,  284 
(1985). 

J.  Chen,  J.  R.  Sites,  I.  L.  Spain,  M.  J.  Hafich,  and  G.  Y. 
Robinson,  Appl.  Phys.  Lett.  58,  744  (1991). 

W.  Shan,  S.  J.  Hwang,  and  J.  J.  Song,  H.  Q.  Hou,  and  C. 
W.  Tu,  Appl.  Phys.  Lett.  62,  2078  (1993), 

M.  J.  Hafich,  H.  Y.  Lee,  G.  Y.  Robinson,  D.  Li,  and  N. 
Otsuka,  J.  Appl.  Phys.  69,  752  (1991). 

A.  Jayaraman,  Rev.  Sci.  Instrum.  57,  1013  (1986),  and 
references  therein. 

S.  M.  Olsthoorn,  F.  A.  J.  M.  Driessen,  and  L.  J.  Giling,  J. 
Appl.  Phys.  71,  2423  (1992). 


D.  P.  Bout,  J.  R.  Shealy,  G.  W.  Wicks,  and  W.  J.  Schaff, 
Appl.  Phys.  Lett.  50,  615  (1987). 

Chen,  Ph.  D.  dissertation,  Colorado  State  University, 
1990. 

C.  S.  Menoni,  H.  D.  Hochheimer,  and  I.  L.  Spain,  Phys, 
Rev.  B  33,  5896  (1986). 

H.  Miiller,  R.  Troemmer,  and  M.  Cardona,  Phys.  Rev.  B 
21,  4879  (1980). 

=^°G.  Olbright,  W.  S.  Fu,  J.  T.  Klem,  J.  M.  Gibbs,  G. 
Khitrova,  R.  Pon,  B.  Fluegel,  K.  Meissner,  N.  Peygham- 
bmian,  R.  Binder,  I.  Galbraith,  and  S.W.  Koch,  Phys.  Rev. 
B  44,  3043  (1991). 

J.  Chen,  D.  Patel,  J.  R.  Sites,  I.  L.  Spain,  M.  J.  Hafich, 
and  G.  Y.  Robinson,  Solid  State  Commun.  75,  693  (1990). 
J.  A.  Tuchman  and  I.  P.  Herman,  Phys.  Rev.  B  45,  11 929 
(1992). 

T.  P.  Pearsall,  Semiconductor  and  Semimetals  (Academic 
Press,  Boston,  1990),  p.  1. 

Chris  G.  Van  de  WaUe,  Phys.  Rev.  B  39,  1871  (1989). 

T.  Y.  Wang,  A.  W.  Kimball,  G.  S.  Chen,  D.  Birkedal,  and 
G.  B.  Stringfellow,  J.  Appl.  Phys.  68,  3356  (1990). 

G.  Bastard,  Wave  Mechanics  Applied  to  Semiconductor 
Heterostructures  (Halsted  Press,  New  York,  1988),  p.  63. 
M.  O.  Watanabe  and  Y.  Ohba,  Appl.  Phys.  Lett.  50,  906 
(1987). 

L.  Y.  Leu  and  S.  R.  Forrest,  J.  Appl.  Phys.  64,  5030  (1988). 
J.  Nunenkamp,  K.  Reimann,  J.  Kuhl,  and  K.  Ploog,  Phys. 
Rev.  B  44,  8129  (1991). 


16 


PHYSICAL  REVIEW  B 


VOLUME  53,  NUMBER  19 


15  MAY  1996-1 


L-band  recombination  in  InjcGai_jP/Ino^AIo.sP  multiple  quantum  wells 
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We  report  the  direct  observation  of  recombination  from  the  band  in  In^Gai-;,P/Ino.5Alo.5P  multiple 
quantum  wells.  The  indirect  transition  is  observed  in  unstrained  structures  with  narrow  wells  and  in  tensile 
strained  structures,  using  high-pressure  photoluminescence  measurements.  Li^  recombination  is  characterized 
by  a  pressure  coefficient  of  60±  5  meV/GPa,  considerably  smaller  than  that  of  the  direct  gap  states.  In  the 
same  experiments  we  also  identify  heterostructure  states  associated  with  T and  Xi^  from  w^ch  we  determine 
the  separation  among  the  conduction  minima  in  unstrained  bulk  In^Gaj.^P  for  x^OAS.  [S0163- 
1829(96)06620-9] 


In^Gai__^P  is  one  of  the  most  attractive  materials  for  the 
development  of  semiconductor  lasers  emitting  in  the  yellow- 
red  region  of  the  optical  spectrum.^  Current  visible  laser  di¬ 
ode  technology  uses  multiple  quantum  wells  (MQW’s),  with 
In^Gai.^P  wells  and  Ino.5Alo.5P  or  In^(Al^.Gai_p^P  bar¬ 
riers.  This  heterostructure  combination  coupled  with  the  pos¬ 
sibility  of  varying  well  width  and  composition  offers  large 
flexibility  in  the  selection  of  the  operating  wavelength  as 
well  as  for  tailoring  of  the  output  characteristics.  Successful 
design  and  optimization  of  advanced  laser  structures  is  based 
on  a  complete  understanding  of  the  band  structure  and  band 
alignments  of  the  heterostructure  materials  composing  the 
laser.  Normally  knowledge  of  the  direct  energy  band  gap  is 
sufficient  for  the  selection  of  the  operating  wavelength. 
However,  for  specific  requirements  such  as  the  reduction  of 
the  laser  threshold  current,  it  is  necessary  to  know  the  sepa¬ 
ration  among  the  conduction-band  minima,  so  as  to  avoid 
any  detrimental  influences  in  the  laser  operation  from  the 
higher  effective  mass,  indirect  and  Xi^  valleys. 

The  band  structure  of  bulk  In^Gai^^^P  has  been  investi¬ 
gated  previously.^"*^  It  is  known  that  In^Ga  1  ^^P  is  a  direct- 
gap  material  for  In  compositions  jc>0.32  and  that  it  becomes 
indirect,  with  being  the  lowest  conduction-band  minima, 
for  In  compositions  a:  <0.32.^’'^  Decreasing  the  In  composi¬ 
tion  further  reveals  the  Xi^  extrema  that  become  the 
conduction-band  minima  for  x<0.2.^’^’^  The  proximity  of 
the  and  Xi^  transitions  has  made  difficult  the  identifica¬ 
tion  of  Lie  -  Previous  optical  studies  in  alloys  with  varying 
compositions  did  not  show  any  evidence  of  the  hand} 
The  L-like  behavior  was  also  absent  in  optical  measurements 
at  high  pressure  conducted  on  In-rich  alloys."^"^  Evidence  of 
Lie  was  obtained  from  modulated  piezoreflectance 
measurements^  and  high-pressure  mobility  measurements^ 
from  which  the  energy  separation  among  the  conduction- 
band  extrema  of  In^^Ga  1  _^P  alloys  was  not  obtained  directly. 
The  proximity  of  the  Tj^,  and  Xi^  extrema  for 

.r<0.48  places  restrictions  in  the  design  of  In^Gai-^P 
MQW’s,  as  indirect  well  structures  can  be  obtained  by  tai¬ 
loring  the  separation  among  the  conduction-band  states  by 
means  of  varying  the  composition  and  the  well  width. 

In  this  paper  we  report  the  observation  of  recombination 
from  the  states  in  In^Gai«^P/Ino.5Alo.5P  MQW’s.  The 
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L-like  behavior  of  the  conduction-band  minima  was  identi¬ 
fied  in  narrow  lattice  matched  MQW’s,  and  also  confirmed 
on  tensile  strained  MQW’s,  from  pressure-dependent  photo¬ 
luminescence  (PL)  measurements  at  low  temperature.  Carrier 
recombination  from  Lj^  was  characterized  by  a  reduced 
pressure  coefficient  compared  to  that  measured  for  the  lowest 
confined  F state.  The  identification  of  the  MQW  states 
associated  with  F and  Xi^  in  the  same  experiments  made 
possible  the  calculation  of  the  separation  among  the 
conduction-band  extrema  in  In^Gai.^^^?  for  compositions 
jc^O.48. 

The  combination  of  using  hydrostatic  pressure  and 
MQW’s  with  different  well  width  and  compositions  was  es¬ 
sential  for  the  unique  determination  of  the  Lj^  transition  in 
In;fGai_;fP/Ino.5Alo.5P  MQW’s.  With  hydrostatic  pressure 
the  relative  position  of  the  conduction-band  extrema  was 
modified  and  the  indirect  states  were  revealed  as  they  be¬ 
came  the  lowest  conduction-band  minima.  The  Lj^,, 
and  Xie  minima  were  distinctly  identified  as  they  shifted 
with  pressure  at  rates  of  about  1(X),  60,  and  -20  meV/GPa, 
respectively,  typical  of  other  semiconductor  materials.^  Since 
both  r  and  Lj^  have  positive  pressure  coefficient,  obtain¬ 
ing  a  band  alignment  in  which  Lj^  becomes  the  conduction- 
band  minimum  at  high  pressures,  requires  the  F u-Lj^  sepa-. 
ration  to  be  considerably  smaller  than  that  of  at 

atmospheric  pressure.  This  condition  was  achieved  in  the 
In^Gai>.jpP/Ino.5Alo.5P  MQW’s,  by  decreasing  the  well 
width  and  in  a  separate  experiment  by  decreasing  the  In  com¬ 
position  X. 

The  In^Ga|_^P/Ino.5AIo.5P  MQW  samples  used  in  these 
studies  were  grown  by  gas-source  molecular  beam  epitaxy 
on  (100)  GaAs  substrates  with  nominally  lattice  matched 
composition  for  the  barrier  material.  The  growth  of  the 
In^Ga  I  .j.P/Ino 5AI0.5P  MQW’s  was  carried  out  at  530  ®C  at 
1.0  /LLvn/h,  conditions  that  produced  a  disordered,  random 
alloy  in  bulk  samples  as  previously  determined  by  photolu¬ 
minescence,,  photoluminescence  excitation  (PLE),  and  trans¬ 
mission  electron  diffraction  measurements.^  Each  sample 
contained  a  buffer  layer  of  GaAs  of  0.5-/um  thickness,  fol¬ 
lowed  by  an  Ino.5Alo.5P  buffer  layer  0.035  /xm  thick,  the 
MQW  region  and  an  Ino.5Alo.5P  cap  layer  0.1  /xm  thick.  All 
layers  were  unintentionally  doped.  The  MQW’s  consisted  of 
50  periods  of  50-A  wells  and  150-A  barriers.  The  well  com¬ 
position  X  was  selected  equal  to  0.48,  0.41,  and  0.37  corre- 
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FIG.  1.  (a)  PL  spectra  from  the  -0.57%  ten¬ 
sile  strain  MQW,  corresponding  to  transitions  as¬ 
sociated  with  Tic,  Xic^  (b)  Schematic 

diagram  of  the  band  structure  of  In^^Gai^^^P/ 
Ino.5Alo.5P  MQW’s  showing  the  various  photo¬ 
luminescence  transitions  observed  in  the  experi¬ 
ments. 


spending  to  0,  -0.57%,  and  ~0.87%  tensile  strain.  For  the 
unstrained  composition  a  second  structure  consisting  of  50 
periods  of  30- A  wells  and  150- A  barriers  was  also  grown. 
The  high-pressure  photoluminescence  measurements  were 
conducted  at  50  K.  Details  of  these  experiments  have  already 
been  described  in  Ref.  10. 

Typical  PL  spectra  originating  from  well  states  associated 
with  Ticj  Lie,  ^ic  shown  in  Fig.  1(a).  The  full 
width  at  half  maximum  of  the  PL  signal  corresponding  to  the 
direct  transition  was  measured  to  be  17  meV.  The  PL  sig¬ 
nals  of  the  indirect  transitions  were  broader  and  weaker  com¬ 
pared  to  the  direct  transition.  We  notice  that  since  the 
states  are  confined  in  the  well,  their  PL  is  more  intense  and 
narrower  than  that  of  .  As  shown  next,  the  identification 
of  the  PL  transitions  was  made  strictly  on  their  pressure  rate 
of  change.  Their  labeling  is  shown  in  the  band  diagram  of 
Fig.  1(b)  and  is  consistent  with  our  previous  work  on  this 
material  system.^® 

Figure  2  shows  the  pressure  dependence  of  the  PL  transi¬ 
tions  measured  in  the  unstrained  Ino.48Gao.52P/Ino.5Alo.5P 
MQW’s  investigated.  The  pressure  behavior  of  both  the 
50-A  andthe30-A  MQW’s  is  qualitatively  similar.  Four  PL 
transitions  were  identified,  one  with  a  positive  pressure  co¬ 
efficient  (dEfdp)  and  three  PL  transitions  with  negative 
slope,  characteristic  of  an  X  band.  Among  the  three  X-like 
transitions,  two  were  observed  over  the  whole  pressure  range 
of  the  experiments.  These  PL  transitions,  labeled  E3  and 
E4,  correspond  to  the  indirect  recombination  of  photoexcited 
carriers  from  the  Xj^  minima  in  the  well  (Xj^-  Ihh)  and  in 
the  barrier,  respectively  (Xicb^r  i^).  The  third  indirect  tran¬ 
sition  (E2)  was  only  observed  at  high  pressures  after  the 
lowest  confined  conduction-band  state  becomes  resonant 
with  the  barrier  X I  cb  states.  E2  corresponds  to  the  indirect,  in 
real  and  k  space,  (Xicb^lhh)  recombination.  As  expected,  the 
onset  of  E2  occurred  at  lower  pressures  in  the  30- A  MQW’s 
as  the  separation  between  the  lowest  confined  state 
(1^)  and  Xicb  is  reduced  with  increased  earner  confinement. 

The  PL  transitions  characterized  by  a  positive  dEidp 
showed  a  distinct  behavior.  In  the  50-A  MQW,  this  transi¬ 
tion  labeled  El  was  characterized  by  dE! dp— 92±3  meV/ 
GPa,  a  value  typical  of  a  F-like  behavior.^  We  associate 
El  with  the  recombination  between  the  conduction-  and 
valence-band  ground  states  (le-lhh).  By  contrast,  dEfdp  in 
the  30-A  MQW’s  was  measured  to  be  60±5  meV/GPa.  Al¬ 


though  an  '-“6%  reduction  in  the  pressure  coefficient  of  El  is 
expected  from  the  variation  in  the  well  width,  this  effect 
cannot  explain  the  almost  50%  decrease  measured  in  the  nar¬ 
row  MQW’s.  The  lower  value  of  dEfdp  indicates  that  in  the 
30-A  MQW’s  recombination  takes  place  from  a  conduction- 
band  state  with  an  L-like  behavior.  Li^-lhh  recombination, 
labeled  El  in  Fig.  1(b),  is  observed  in  the  narrow  well  struc¬ 
tures  since  El  is  shifted  above  Lie  increased 

earner  confinement.  The  calculated  pressure  dependence  of 
El  in  the  30-A  MQW  is  represented  by  the  solid  line  in  Fig. 
2(b).  In  the  50-A  MQW  the  El  transition  was  not  observed, 
as  it  occurred  at  higher  pressures  where  the  indirect  E2  tran¬ 
sition  becomes  dominant.  The  El  behavior  in  the  50- 
A  MQW,  calculated  assuming  the  same  pressure  coefficient 


FIG.  2.  Pressure  dependence  of  the  photoluminescence  transi¬ 
tions  of  unstrained  Ino.48Gao.52P/Ino.5Alo.5P  multiple  quantum  wells 
of  (a)  50-A  well  width  and  (b)  30-A  well  width.  The  pressure  be¬ 
havior  of  all  transitions  has  been  fitted  using  a  least-square  routine 
from  which  their  pressure  rate  of  change  and  the  atmospheric  pres¬ 
sure  values  are  determined. 
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FIG.  3.  Energy-pressure  diagram  for  tensile  strained  MQW’s 
with  (a)  -0.57%  strain  and  (b)  -0.87%  strain. 

as  that  measured  in  the  30-A  MQW,  is  shown  by  the  dash- 
dotted  line  in  Fig.  2(a). 

To  confirm  our  assignment  of  we  examined  the  re¬ 
sults  of  the  tensile  strained  MQW’s.  The  L-like  behavior 
should  also  be  apparent  in  these  structures  as  the  separation 
between  Lj-  and  Ti^  decreases  with  decreasing  In 
composition.^  Figure  3  shows  the  energy-pressure  data  of  the 
In;cGai_j.P/Ino.5Alo.5P  MQW’s  with  well  compositions 
x=0,4i  and  0.37,  corresponding  to  —0.57%  and  —0.87% 
strain,  respectively.  Hereafter,  we  discuss  the  behavior  of  the 
transitions  with  positive  dEfdp,  as  the  indirect  transitions 
are  equivalent  to  those  observed  in  the  unstrained  MQW’s. 
In  the  -0.87%  strain  MQW’s,  dEfdp  was  measured  to  be 
60±5  meV/GPa,  identical  to  that  of  LI  in  the  30-A  un¬ 
strained  MQW’s.  In  the  —0.57%  strain  MQW’s  the  PL  tran¬ 
sition  with  positive  dEfdp  showed  at  first  glance  a  distinct 
behavior.  Below  0.5  GPa  it  shifted  at  a  rate  dEldp=\\0±5 
meV/GPa  and  above  0.5  GPa  at  a  slower  rate  dEfdP 
=60±5  meV/GPa,  as  shown  by  the  solid  and  dash-dotted 
lines  of  Fig.  3(a).  This  change  in  the  slope  indicates  that  a 
pressure-induced  £1-L1  crossover  takes  place  at  0.5  GPa. 
The  low  pressure  at  which  this  crossover  occurred  shows  that 
E\  and  LI  are  separated  approximately  35  meV  at  ambient 
conditions.  With  increased  tensile  strain,  the  £1 -LI  separa¬ 


tion  was  decreased  and  L 1  became  the  lowest  conduction- 
band  state,  as  observed  in  the  —0.87%  strained  MQW, 

The  identification  of  the  MQW  states  associated  with  the 
r ic ,  Lie »  ^i^d  Xie  extrema  in  the  In^Ga  i  _^P  wells  allowed 
the  calculation  of  the  energy  separation  among  the 
conduction-band  extrema  in  In^cGai-^^P  forx^  0.48.  In  ref¬ 
erence  to  Fig.  1(b),  the  separation  was  calculated 

from  the  difference  of  LI  and  £1  at  atmospheric  pressure 
and  by  adding  the  le  confinement  energy.^**  We  neglected 
the  LI  electron  confinement  energy  since  it  is  small  due  to 
the  large  effective  mass  of  the  L]^  band.^^  The  T  sepa¬ 
ration  was  determined  from  the  difference  between  £3  and 
£1  at  atmospheric  pressure  and  by  adding  the  le  confine¬ 
ment  energy^  and  the  strain-induced  splitting  of  Xi^,  calcu¬ 
lated  using  the  model  solid  theory.  The  values  of  £1  in  the 
MQW’s  with  composition  x=0.37  and  .r==0.48  (30  A  well 
width),  in  which  this  transition  was  not  observed,  were  de¬ 
termined  from  PLE  measurements.^^ 

In  unstrained  bulk  Ino.48G^o.52P>  was  calculated  to  be 
0.1  ±0,02  eV  above  T and  0.1 8± 0.04  eV  below  .  De¬ 
creasing  the  In  composition  x  decreased  the  F  i^-Li^  separa¬ 
tion  to  0.08±0.02  and  0.03±0.02  eV  for  A'  =  0.4i  and  0.37, 
respectively.  The  separation  was  detremined  to  be 

0.21  ±0.05  eV  for  0,48.  For  the  lower  In  compositions 
its  value  was  equal  to  0.175±0.05  eV.  The  insensitivity  of 
our  results  for  determining  the  separation  of  Xi^  relative  to 
the  other  conduction-band  extrema  arises  from  the  uncertain¬ 
ties  involved  in  calculating  the  strain-induced  shift  of 
From  the  measured  variation  of  the  Fj^-Li^  energy  separa¬ 
tion  with  In  composition  and  assuming  a  linear  behavior  in 
the  range  of  compositions  investigated,  we  determined  that 
Lie  becomes  the  conduction-band  minima  atx  =  0.3.  A  simi¬ 
lar  value  was  predicted  by  Bugaski,  Kontkiewicz,  and 
Mariette^^  and  also  determined  from  Hall  measurements  at 
high  pressure^  and  piezomodulated  measurements.^  Al¬ 
though  the  agreement  in  the  F  le-Lie  crossover  composition 
with  previous  work^’^  is  very  gooi  the  values  of  the  Fj^ 
band  gap  used  as  reference  in  this  work  are  found  to  be 
approximately  60  meV  smaller  than  that  reported  by  Merle 
et  al}  While  *^20  meV  can  be  accounted  for  by  the  localiza¬ 
tion  energy,  the  rest  is  likely  to  arise  from  differences  in  the 
growth  conditions. 

In  conclusion,  we  have  identified  the  L^^  band  in  narrow 
unstrained  and  in  tensile  strained  In;pGai_^P/Ino.5Alo.5P 
MQW’s  and  determined  from  these  results  the  separation  of 
the  conduction-band  extrema  in  In^Gaj^^P  for  In  composi¬ 
tions  x^O.48. 
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Abstract 

The  epitaxial  growth  and  characterization  of  multiple  quantum  well  (MQW)  optical  modulators  on  GaP  substrates 
is  reported.  The  devices  were  fabricated  using  graded  composition  InGaP  strain  relieving  buffer  layers  and  strain- 
compensated  InGaP  MQWs.  The  devices  operate  in  the  570  nm  to  620  nm  wavelength  range  with  a  maximum 
transmission  change  of  15%  occurring  at  594  nm  for  an  applied  bias  of  -10  V.  The  composition  range  required 
for  the  target  wavelength  range  is  near  the  indirect-direct  crossover  point  for  the  InGaP  material  system  and  the 
corresponding  effects  on  photoluminescence  and  absorption  data  is  presented. 


The  current  interest  in  plastic  optical  fibers, 
displays,  and  optical  storage  has  spurred  interest  in 
optical  devices  operating  in  the  visible  wavelength 
range.  The  wavelength  range  from  550  nm  to  670  nm 
is  of  particular  interest  for  plastic  fiber 
commxinications  since  the  lowest  fiber  absorption 
occurs  at  570  and  650  nm.  In  addition,  two-way 
transmission  networks  using  only  one  optical  source 
and  many  optical  modulators  operating  in  reflection 
mode  have  been  proposed  to  minimize  the  use  of  low 
reliability  and  costly  components  in  robust  networks 
[1].  To  date  visible  modulators  have  been  fabricated 
in  both  the  III-V  and  II-VI  semiconductors  at 
wavelengths  ranging  from  506  nm  to  660  nm  [2-9]. 
Recently,  Blum  et  al,  reported  observation  of  the 
quantum  confined  Stark  effect  (QCSE)  in  surface 
normal  InGaAlP/InGaP  multiple  quantum  well 
structures  grown  on  GaAs  by  metal-organic  vapor 
phase  epitaxy  and  operating  at  645  nm  [2].  In 
addition,  Goosen  et  al.  reported  MQW  modulators 
operating  at  575  nm  using  indirect  AlAs/AlGaAs  in 
the  active  region  grown  on  GaAs  [10].  Here  we 
report  the  fabrication,  by  gas  source  molecular  beam 
epitaxy  (GSMBE),  of  multiple  quantum  well  (MQW) 
optical  modulators  for  operation  near  570  mn.  The 
modulators  consist  of  strain-compensated  In^Gai. 
2P(weli)/InyGai.yP(barrier)  MQWs  grown  on  a  thick 
buffer  layer,  where  y<x<z,  all  of  which 
were  grown  on  a  transparent  GaP  substrate  using  a 


buffer  layer  of  graded  composition. 

A  graded  buffer  layer  was  used  to  reduce  the 
density  of  misfit  dislocations  in  the  MQW  region. 
Several  graded  buffer  schemes  were  evaluated,  on 
both  GaP  and  GaAs  substrates,  using  double  crystal 
X-ray  diffraction  (DCXD)  and  photoluminescence 
(PL).  Constant-composition,  step-graded,  linear- 
graded,  and  superlattice-graded  buffers  were  tested 
with  the  step-graded  scheme  producing  marginally 
better  results  than  the  linear-graded  scheme  and 
significantly  better  than  the  superlattice-graded 
approach.  The  graded  buffer  layers  grown  on  GaP 
were  consistently  of  higher  structural  and  optical 
quality  than  those  on  GaAs. 

The  MQW-PIN  device  structure  xmder  study 
is  shown  in  Fig.  1  and  was  grown  on  a  (100)  n'^GaP 
substrate  using  thermally  cracked  PH3.  The  In  mole 
fraction  in  the  thick  buffer  layer  of  constant 
composition,  denoted  as  x,  was  varied  systematically 
in  a  series  of  growths.  The  composition  of  the  wells 
z  and  the  barriers  y  were  varied  along  with  x  such 
that  z  =  x+0.15  and  y  =  x-0.10.  The  GaP 
substrates  were  degreased  and  etched  in  a  4:4:5 
HCL:HN03:H20  solution  at  25  ®C  for  2  minutes  and 
then  soaked  for  10  minutes  in  HjO  to  prepare  the 
surface  prior  to  growth.  Surface  oxide  was  desorbed 
at  700  °C  for  30  minutes  under  a  PH3  flow  of  8  seem. 
Reflection  high  energy  electron  diffraction  (RHEED) 
gave  a  characteristic  (2x4)  pattern.  An  improved 
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(i.e.,  more  streaked)  RHEED  pattern  was  typically 
observed  after  depositing  a  few  monolayers  of  Ga  and 
annealing  at  660°C  for  30  seconds  under  the 
phosphorous  flux.  The  step-graded  buffer  layer 
consisted  of  three  layers  of  InGaP  of  different 
composition,  each  200-nm  thick,  followed  by  a  2.0- 
|xm  thick  constant  composition  In^Gaj,^  layer.  The 
MQW  active  region  consisted  of  50  periods  of  7.5- 
nm  InyGaj.yP  wells  and  11 -nm  In^Ga^.^?  barriers. 
Compensating  strain  in  the  wells  and  barriers  was 
used  to  avoid  introduction  of  additional  defects  that 
may  result  from  the  strain  in  the  MQWs.  The 
devices  were  etched  mesas  of  400x400)im^  with  Au- 
Sn-Au  and  Ti-Au  ohmic  contacts  for  the  n"^  and  p"^ 
layers,  respectively. 

One  of  the  issues  in  using  the  alloy  InGaP 
for  devices  operating  in  the  550-670  nm  range  is  the 
crossover  in  the  conduction  band  from  indirect  to 
direct  for  In  mole  fractions  greater  than  about  0.30 
[12].  For  the  range  of  x  used  in  this  study 
(0.25<x<0.35),  the  wells  were  direct  (0.40<z<0.50) 
and  the  barriers  were  indirect  (0.15<y<0.25).  The 
well  composition  and  thickness  was  chosen  to  match 
the  target  wavelength,  including  any  strain  and 
confinement  effects  as  determined  from  the  model 
solid  theory  [11].  The  barrier  composition  was 
chosen  to  balance  the  strain  in  the  well,  giving  an 
average  lattice  constant  for  the  MQW  stack  equal  to 
that  of  the  buffer  layer.  The  buffer  layer  composition 
was  chosen  so  as  to  allow  the  least  amount  of  strain 
required  in  the  well  while  at  the  same  time 
minimizing  absorption  in  the  buffer  at  the  wavelength 
of  interest. 

A  typical  PL  spectra  is  shown  in  Fig.  2(a) 
for  the  modulator  structure.  Two  peaks  are  observed: 
one  at  565  nm  corresponding  to  emission  from  the 
cap  and  buffer  layers  with  x=0.28  and  one  at  595  nm 
for  emission  from  the  QWs.  The  peak  locations  are 
in  agreement  with  the  calculated  energies.  The  effect 
of  the  direct-indirect  crossover  was  evident  in  the 
variation  of  PL  intensity  as  x  was  varied  from  sample 
to  sample.  Fig.  2(b)  shows  a  strong  decrease  in  PL 
intensity  with  decreasing  wavelength  as  the 
composition  of  the  buffer  layer  approaches  the 
crossover  point.  This  effect  is  consistent  with  a 
decrease  in  the  number  of  electrons  available  for 
radiative  recombination  due  to  the  indirect  valley 
becoming  more  populated  as  x  decreases.  Likewise, 
the  PL  intensity  from  the  QWs  decreases  at  shorter 
wavelengths.  The  variation  in  the  QW  PL  intensity 
could  result  from  the  electrons  in  the  well  occupying 
the  indirect  minima  in  the  well.  Since  our 


calculations  suggest  that  the  InyGai,yP/In2Gai.JP 
interface  is  a  type-II  heterojunction,  electrons  in  the 
well  could  also  transfer  to  the  indirect  minima  in  the 
barrier.  Based  on  our  calculations,  the  slopes  of  the 
curves  in  Fig.  2(b)  for  the  bulk  and  QW 
luminescence  are  expected  to  differ  if  transfer  to  the 
barrier  indirect  minima  is  the  dominating  effect  for 
the  QWs.  Since  the  slopes  are  the  same  within  the 
experimental  imcertainty ,  it  appears  that  an  increasing 
occupation  of  the  indirect  minima  in  the  well  is  the 
dominant  mechanism  for  the  decrease  of  PL  intensity 
as  X  decreases. 

Optical  transmission  data  was  taken  to 
estimate  the  absorption  coefficient  in  the  MQW 
region  and  to  determine  the  feasibility  of  InGaP  on 
GaP  for  normal-incident  optical  modulators.  Since 
the  GaP  substrate  is  transparent  for  wavelengths  >545 
nm,  the  transmission  measurement  was  made  without 
removal  of  the  substrate.  The  absorption  spectra  at 
room  temperature  for  three  different  modulators  are 
shown  in  Fig.  3(a).  Each  exhibits  features 
characteristic  of  QW  excitonic  absorption,  vrith  the 
peak  at  594  nm  in  the  sample  with  x==0.30  due  to  the 
le-hh  transition,  as  based  on  our  calculations.  A 
shoulder  in  the  spectra  at  ^12  nm  below  the  le-hh 
transition  corresponds  to  the  expected  position  of  the 
le-lh  transition.  The  strong  absorption  edge  observed 
--30  nm  below  the  le-hh  peak  corresponds  to  that 
expected  for  the  buffer  layer.  Fig.  3(b)  shows 
absorption  coefficient  data  for  the  same  samples  at 
low  temperature  (10  K).  Measurements  at  low 
temperature  confirmed  carrier  confinement  in  the  well 
and  provided  estimates  of  the  half  width  of  the 
excitonic  resonance.  Half  widths  for  the  le-hh 
resonance  at  10  K  were  typically  15  to  20  meV  for 
the  range  of  wavelengths  covered  by  the  samples 
grown,  with  the  narrowest  peak  found  for  the  sample 
operating  at  594  nm  at  room  temperature. 

The  smaller  le-hh  peak  amplitude  in  Fig.  3 
for  the  sample  with  x=0.25  compared  to  that  for 
x=0.30  may  result  from  the  direct-indirect  crossover. 
Not  only  is  the  well  material  closer  to  indirect  at 
shorter  wavelengths,  but  calculations  suggest  that  the 
indirect  barrier  conduction  band  edge  is  near  or  below 
the  direct  confined  energy  level  in  the  well.  Such  a 
condition  would  result  in  a  reduction  of  exciton 
lifetime  [3]  and  a  broadening  of  the  excitonic 
absorption  peak,  in  agreement  with  our  observations. 
We  would  then  expect  the  sample  with  x=0.35  to 
show  a  more  pronounced  exciton  absorption  peak 
than  for  x=0.30;  as  is  evident  from  Fig.  3,  this  is  not 
the  case.  Alternatively,  the  material  quality  could 
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suffer  as  x  is  increased  because  of  a  larger  defect 
density  produced  by  the  increased  mismatch  with  the 
GaP  substrate.  We  suspect  that  the  material  quality 
is  not  the  dominant  effect  causing  the  weakening  and 
broadening  of  the  excitonic  resonance  at  longer 
wavelengths,  since  additional  efforts  to  improve  the 
material  quality  by  changing  the  buffer  layer 
produced  improved  DCXD  linewidths  but  gave  no 
significant  improvement  in  the  absorption  data. 

Voltage  dependent  absorption  and 
transmission  data  are  shown  in  Fig.  4(a)  for  a  device 
with  x=0.30,  which  corresponds  to  room  temperature 
modulation  at  594  nm.  The  excitonic  absorption 
peak  shows  a  weakening,  broadening,  and  a  red  shift 
with  increasing  reverse  bias.  The  shift  in  the  peak  to 
longer  wavelengths  was  most  distinct  for  the  594-nm 
device  and  was  consistent  with  that  expected  for  the 
QCSE.  A  change  in  the  absorption  coefficient  of 
1700  cm"^  was  measured  for  devices  operating  at  594 
nm.  Differential  transmission  data  is  shown  in  Fig. 
4(b).  A  maximum  change  in  transmission  of  15%  at 
-lOV  bias  was  measured  at  594  nm. 

In  summary,  the  growth  and  operating 
characteristics  of  MQW  modulators  in  the  InGaP 
material  system  fabricated  on  transparent  GaP 
substrates  is  reported  for  the  first  time.  Using 
GSMBE  the  devices  were  grown  on  graded 
composition  buffer  layers  and  employed  strain 
compensation  in  the  MQW  region.  PL  intensity 
variation  with  composition  showed  clear  evidence  of 
the  direct-to-indirect  crossover  in  the  InGaP  QWs 
designed  to  operate  in  the  570-620  nm  range.  The 
best  devices  produced  a  transmission  change  of  15% 
with  a  modest  10  V  bias  at  a  wavelength  of  594  nm. 
Operation  at  shorter  wavelengths  was  observed  but 
with  weaker  excitonic  effects. 

The  authors  gratefully  acknowledge  the 
support  of  the  Air  Force  Office  of  Scientific  Research 
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FIG.  1.  Schematic  representation  of  the  p-i-n 
diode  layer  structure  used  in  this  study.  An  InGaP 
step  composition  graded  buffer  was  used  to 
reduce  dislocation  densities  in  the  MQW  active 
region. 
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FIG.  2.  Measured  (PL)  data  for  the  MQW 
samples  in  this  study,  (a)  PL  spectra  for  the 
device  with  a  nominal  buffer  composition  of 
x=0.30.  Luminescence  from  the  bulk  cap  layer 
and  MQW  is  observed,  (b)  Variation  in  PL 
intensity  at  different  wavelengths  for  both  the 
cap  layer  and  MQW.  The  decrease  in  intensity 
with  decrease  in  wavelength  is  expected  as  the 
composition  approaches  the  indirect-direct 
crossover. 
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Nonresonant  tunneling  rates  have  been  measured  in  InGaP/InAlP  asymmetric  double 
quantum-well  structures  for  which  optical  phonon  assisted  tunneling  is  energetically  forbidden. 

For  an  initial  photoexcited  carrier  density  of  2.4  X 10^*  cm‘“^,  tunneling  times  of  220,  60,  and 
less  than  9  ps  have  been  measured  in  samples  with  barrier  thickness  4.5,  3.0,  and  1.5  nm, 
respectively.  The  tunneling  times  were  found  to  be  strongly  dependent  on  carrier  density.  The 
measured  tunneling  times  and  their  dependence  on  carrier  density  are  compatible  with  impurity 
scattering  being  the  dominant  mechanism  assisting  the  tunneling. 


Tunneling  phenomena  in  semiconductor  structures  is 
being  actively  investigated  for  fundamental  reasons  and 
possible  device  applications.*’^  Recently,  resonant  and  non¬ 
resonant  carrier  tunneling  has  been  studied  in  GaAs  based 
asymmetric  double  quantum-well  (ADQW)  structures^”** 
in  part  motivated  by  their  potential  use  in  optical  switching 
and  computing.^’ *°  Herein,  we  report  the  first  tunneling 
measurements  in  InGaP/InAlP  asymmetric  double  quan¬ 
tum  wells  which  are  suitable  for  the  design  of  modulators 
that  operate  in  the  visible  region  of  the  spectrum.  In  this 
work,  nonresonant  carrier  tunneling  was  investigated  in 
ADQW  structures  where  optical  phonon  assisted  tunneling 
is  energetically  forbidden.  Steady  state  and  time  resolved 
photoluminescence  (PL)  measurements  were  conducted 
on  samples  with  dilferent  barrier  thicknesses  as  a  function 
of  carrier  density  and  temperature. 

The  samples  studied,  grown  by  gas  source  MBE,  were 
unintentionally  n  doped  and  consisted  of  15  periods  of 
ADQW  structures  separated  by  an  InAlP  layer  of  15  nm. 
Each  period  consists  of  InGaP  wells  of  6  and  12  nm 
widths,  separated  by  an  InAlP  barrier  having  a  different 
width  for  each  sample:  1.5,  3.0,  4.5,  and  15  nm.  The  layer 
thicknesses  are  nominal  values  with  an  error  of  less  than 
10%.  The  structure  with  a  15  nm  barrier,  for  which  no 
carrier  tunneling  is  expected  to  occur,  was  used  as  a  con¬ 
trol  sample. 

The  epitaxial  structures  were  grown  on  (100) 
72 -h  GaAs  substrates  at  approximately  1  /zm/h  at  500 
using  thermally  cracked  ASH3  and  PH3,  conditions  that 
have  previously  been  shown  to  produce  abrupt  interfaces 
and  high  quality  quantum  wells.*^’*^  Each  sample  consisted 
of  a  GaAs  layer  (0.25  /zm  thickness)  followed  by  15  peri¬ 
ods  of  the  ADQW  structure  and  capped  by  an  InAlP  layer 
(0.5  /zm  thickness).  At  each  heterointerface  of  the 
ADQWs,  the  growth  was  paused  for  18  s  without  inter¬ 
ruption  of  the  PH3  flow.  X-ray  diffraction  indicated  the 
alloy  compositions  to  be  Ino4gGao  52P  and  Ino.47Alo.53P, 
with  lattice  match  to  the  GaAs  substrate  to  within  500 
ppm.  Previous  measurements  indicate  an  impurity  concen¬ 
tration  of  (2-5)  X 10*^  cm“^  in  the  InGaP  wells.*^ 

In  the  PL  measurements,  the  samples  were  excited  by 
a  train  of  2  ps  pulses  produced  by  a  Rhodamine  560  dye 
laser  synchronously  pumped  by  a  frequency  doubled 


mode-locked  Nd-YLF  laser  at  a  repetition  rate  of  76  MHz. 
The  excitation  energy,  around  2.14  eV,  was  selected  to 
excite  both  wells  simultaneously  and  to  be  sufficiently  sep¬ 
arated  from  the  photon  energies  corresponding  to  the  PL 
of  the  two  quantum  wells  to  facilitate  detection.  The  laser 
beam  was  focused  on  the  sample  with  a  spot  diameter  of 
approximately  36  /zm.  Steady  state  PL  spectra  were  taken 
with  a  0.75  m  focal  length  spectrometer  and  a  Hamamatsu 
R2228  photomultiplier  using  standard  lock-in  techniques. 
The  time  resolved  PL  was  measured  by  a  synchroscan 
streak  cameras  (Hamamatsu  Ml 955)  coupled  to  a  ther- 
moelectrically  cooled  two-dimensional  CCD  array  detec¬ 
tor.  The  streak  camera  was  attached  to  the  exit  of  a  0.3  m 
focal  length  spectrometer  provided  with  a  150  g/mm  grat¬ 
ing.  The  time  resolution  of  the  system  was  determined  to 
be  9  ps.  The  samples  were  mounted  in  a  closed  cycle  re¬ 
frigerated  cryostat  for  low-temperature  measurements. 

Figure  1  shows  the  steady  state  PL  signal  correspond¬ 
ing  to  the  four  ADQW  structures  with  different  barrier 
thicknesses  obtained  at  1 1  K.  All  spectra  show  two  peaks: 
the  stronger  one,  at  lower  energy,  corresponds  to  radiative 
recombination  in  the  wide  well  and  the  weaker  peak,  at 
higher  energy,  to  the  same  process  in  the  narrow  well.  The 
decreasing  PL  intensity  for  the  narrow  well  with  decreas¬ 
ing  barrier  thickness  suggests  that  carriers  tunnel  from  the 
narrow  well  to  the  wide  well.  This  is  confirmed  by  the  time 
resolved  PL  measurements  which,  as  discussed  below, 
show  a  strong  dependence  of  the  time  decay  of  the  narrow 
well  PL  with  barrier  thickness. 

From  the  steady  state  PL  data  of  Fig.  1  and  the  values 
for  effective  masses,  the  band  gap  from  kef  13  and  the 
band  offset  given  in  Ref  14,  we  estimated  the  energy  dif¬ 
ference  between  levels  in  the  wells.  For  all  the  samples  in 
which  tunneling  can  occur,  these  calculations  yield  a  dif¬ 
ference  between  the  first  energy  levels  in  the  conduction 
band  of  the  two  wells  of  less  than  34  meV.  This  energy 
difference  is  more  than  10  meV  below  the  LO  phonon 
energy  for  InGaP  having  the  same  composition  as  in  our 
samples.*^  This  implies  that,  in  these  samples,  LO  phonon 
emission  is  not  energetically  available  to  assist  the  tunnel¬ 
ing  process.  Next,  we  discuss  the  results  of  time  resolved 
PL  measurements. 

Time  resolved  PL  signals  obtained  at  1 1  K  are  shown 
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FIG.  1.  CW  PL  at  1 1  K  of  ADQWs  with  different  barrier  thickness,  d.  A 
significant  drop  in  the  intensity  of  the  PL  peak  associated  with  the  narrow 
well  is  observed  when  the  barrier  thickness  is  reduced. 
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in  Fig.  2  as  a  function  of  barrier  thickness.  Figures  2(a) 
and  2(b)  illustrate  the  temporal  evolution  of  the  PL  peak 
corresponding  to  the  narrow  and  the  wide  well,  respec¬ 
tively.  For  the  sample  with  the  15  nm  barrier,  the  time 
decay  of  the  PL  from  the  narrow  well  is  only  slightly 
shorter  than  that  corresponding  to  the  wide  well,  which  is 
about  300  ps.  As  the  barrier  thickness  is  decreased,  the 
decay  time  of  the  narrow  well  PL  is  observed  to  be  consis¬ 
tently  shorter.  Values  of  130,  50,  and  less  than  9  ps  (the 
time  resolution  of  our  system)  were  measured  for  the 
ADQWs  with  barriers  of  4.5,  3.0,  and  1.5  nm,  respectively. 
In  contrast,  the  decay  time  of  the  wide  well  PL  was  found 
to  be  independent  for  all  the  samples,  except  for  the 
ADQW  with  a  1.5  nm  barrier  that  showed  a  shorter  decay, 
probably  due  to  a  variation  in  sample  quality. 

The  decay  time  of  the  wide  well  PL  is  due  to  both 
radiative  and  nonradiative  recombination.  In  the  narrow 
well,  however,  carriers  can  also  relax  through  tunneling  to 
the  wide  well.  The  tunneling  time  can  then  be  calculated 
from  the  measured  values  of  the  narrow  well  PL  decay 
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FIG.  2.  Time  resolved  PL  signals  corresponding  to  the  PL  peak  wave¬ 
length  of  the  narrow  well  (a)  and  the  wide  well  (b)  for  ADQWs  with 
different  barrier  corresponding  to  an  initial  photoexcited  sheet  carrier 
density  of  2.4 X  lO”  cm“^.  The  measurements  were  conducted  at  1 1  K. 


FIG.  3.  Low  temperature  (II  K)  time  resolved  PL  signals  for  the 
ADQW  with  a  3  nm  barrier  as  a  function  of  excitation  intensity  for  both 
the  narrow  well  (a)  and  the  wide  well  (b). 


time  and  the  narrow  well  recombination  time.  To  calculate 
the  narrow  well  recombination  time,  we  measured  the  de¬ 
cay  rate  of  the  corresponding  wide  well  and  corrected  this 
value  for  the  difference  in  well  thickness  according  to  the 
ratio  of  the  PL  decay  times  from  both  wells  in  our  control 
sample.  For  an  initial  photoexcited  carrier  concentration  of 
2.4 X  lO"  cm“^,  the  tunneling  times  were  calculated  to  be 
approximately  220,  60,  and  less  than  9  ps  for  the  samples 
with  4.5,  3.0  and  1.5  nm  barriers,  respectively.  These  rel¬ 
atively  long  tunneling  times  are  consistent  with  the  fact 
that  LO  phonon  assisted  tunneling  is  not  available  in  these 
samples.  Acoustic  phonon  scattering  can  also  be  ruled  out 
as  the  dominant  mechanism  assisting  the  tunneling.  Anal¬ 
ysis  of  the  time  resolved  PL  spectra  shows  that  the  intra¬ 
band  carrier  relaxation  time  due  to  acoustic  phonon  scat¬ 
tering  is  of  tens  of  picoseconds.  The  associated  interwell 
transfer  time  is  expected  to  be  several  orders  of  magnitude 
larger  due  to  the  reduced  overlap  of  the  wave  functions. 

In  the  absence  of  phonon  assisted  tunneling,  the  scat¬ 
tering  mechanisms  that  may  contribute  to  tunneling  are 
electron-electron,  acoustic  phonon,  impurity,  interface- 
roughness,  and  alloy  disorder.  To  gain  additional  informa¬ 
tion  on  the  role  of  different  scattering  mechanisms,  we 
measured  the  tunneling  rate  dependence  on  carrier  density 
and  temperature.  Time  dependent  PL  measurements  were 
performed  for  laser  intensities  ranging  from  0.08  to  10.5 
mW  corresponding  to  sheet  carrier  densities  between 
8x10^  and  1x10*^  cm”^  Figure  3  shows  the  measured 
dependence  of  the  PL  decays  for  the  ADQW  with  a  3  nm 
barrier  as  a  function  of  excitation  intensity.  For  the  wide 
well,  the  PL  decay  was  relatively  independent  of  carrier 
density.  In  contrast,  the  narrow  well  PL  decay  time  in¬ 
creased  approximately  four  times  as  the  intensity  was  in¬ 
creased  in  the  range  mentioned  above.  A  similar,  but 
weaker,  dependence  on  excitation  intensity  was  measured 
in  the  narrow  well  PL  decay  of  the  sample  with  the  4.5  nm 
barrier.  For  the  sample  with  the  1.5  nm  barrier  the  time 
decay  of  the  narrow  well  PL  remained  below  the  temporal 
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resolution  of  our  detection  system  for  intensities  below  6.5 
mW  and  for  intensities  above  this  value  the  PL  decay  was 
also  observed  to  increase.  Local  heating  of  the  samples 
caused  by  the  increased  intensity  may  affect  PL  decay  mea¬ 
surements.  However,  this  effect  can  be  ruled  out  as  a  pos¬ 
sible  cause  of  an  increase  in  carrier  lifetime  because  mea¬ 
surements  done  as  a  function  of  temp>erature  showed  the 
opposite  behavior:  the  PL  decay  of  the  wide  well  was  ob¬ 
served  to  decrease  slightly  as  a  function  of  temperature. 
Hence,  we  conclude  that  the  scattering  processes  that  assist 
tunneling  are  carrier  density  dependent. 

A  possible  cause  of  the  observed  decrease  of  the  tun¬ 
neling  rate  with  the  increased  photoexcited  density  is 
space-charge  buildup  caused  by  the  significantly  larger  tun¬ 
neling  rates  for  electrons  compared  to  that  of  holes,  due  to 
their  large  difference  in  effective  mass.  This  difference  in 
tunneling  rates  can  cause  a  charge  buildup  between  the  two 
wells,  producing  an  electric  field  which  increases  the  po¬ 
tential  barrier  between  the  two  wells,  opposing  tunneling. 
We  have  calculated  that  the  maximum  potential  change  for 
our  largest  carrier  density  power  to  be  5  meV.  In  this 
calculation  we  assumed  the  extreme  case  in  which  all  of  the 
electrons  and  none  of  the  holes  have  tunneled.  Since,  in 
practice,  the  change  in  the  potential  barrier  is  much 
smaller  than  this  limit  value,  which  is  small  compared  with 
the  energy  separation  between  the  initial  and  final  states, 
we  conclude  that  space-charge  effects  are  not  responsible 
for  the  observed  density  dependence.  An  analysis  of  the 
shape  of  the  measured  PL  decay  as  a  function  of  time  for 
difference  excitation  intensities  also  shows  no  evidence  of 
space-charge  effects. 

The  increase  of  the  tunneling  time  with  carrier  density 
rules  out  electron-electron  scattering  a  the  dominant  tun¬ 
neling  mechanism.  The  remaining  three  processes  have  a 
dependence  on  carrier  density  which  is  compatible  with 
our  observations.  However,  calculated  tunneling  times  due 
to  alloy  scattering,  using  the  potential  given  by  Mukhopad- 
hyay  and  Nag,^^  resulted  in  times  larger  than  1  ns.  Tun¬ 
neling  times  due  to  interface  roughness,  calculated  by  Fer¬ 
reira  and  Bastard,*^  were  found  to  be  between  one  and  two 
orders  of  magnitude  larger  than  for  impurity  scattering 
when  the  concentration  of  scatterers  is  the  same  as  that  of 
impurities.  We  have  calculated  the  nonresonant  tunneling 
rate  for  impurity  scattering  following  the  procedure  out¬ 
lined  by  Ferreira  and  Bastard.'^  Assuming  a  density  of 
impurities  of  5X10^^  cm”^,  and  that  these  impurities  are 
localized  mostly  at  the  interfaces,  we  calculated  nonreso¬ 
nant  electron  tunneling  times  of  450,  70,  and  12  ps  for 
barrier  widths  4.5,  3.0,  and  1.5  nm,  respectively.  These 
calculated  values  are  very  similar  to  the  measured  tunnel¬ 
ing  times.  This  suggests  that  the  mechanism  assisting  the 
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tunneling  is  likely  to  be  impurity  scattering.  A  similar  de¬ 
crease  of  the  nonresonant  tunneling  rate  as  a  function  of 
the  density  of  photoexcited  carriers  was  recently  observed 
in  GaAs/AlGaAs  ADQW’s.  *  *  The  measured  tunneling  de¬ 
pendence  in  that  experiment  was  found  to  be  compatible 
with  impurity  assisted  interwell  scattering,  and  it  was  also 
suggested  that  excitons  may  tunnel  as  a  single  entity. 

In  summary,  we  have  measured  for  the  first  time  non¬ 
resonant  tunneling  times  in  InGaP/InAlP  ADQW  in 
which  LO  phonon  assisted  tunneling  is  not  allowed.  The 
tunneling  times  were  found  to  increase  as  a  function  of  the 
photoexcited  carrier  density.  The  magnitude  of  the  tunnel¬ 
ing  times  and  their  dependence  on  carrier  density  are  com¬ 
patible  with  impurity  scattering  being  the  dominant  tun¬ 
neling  mechanism. 

The  authors  wish  to  acknowledge  Professor  J.  Mahan 
for  lending  the  closed  cycle  refrigeration  system.  Dr.  D. 
Patel  and  M.  Prasad  for  their  assistance  in  the  experiment, 
and  Dr.  K.  Bhattacharyya  for  helpful  discussions.  This 
work  was  supported  by  the  Center  for  Optoelectronic 
Computer  Systems  through  National  Science  Foundation 
Grant  ECD  9015128  and  by  the  Colorado  Advanced  Tech¬ 
nology  Institute.  C.  S.  Menoni  acknowledges  here  support 
by  the  Engineering  Foundation  AFOSR  Research  Initia¬ 
tion  Grant  RI-B-91-16. 


*  F.  Capasso,  K.  Mohammed,  and  A.  Y.  Cho,  IEEE  J.  Quantum  Elec¬ 
tron.  QE-22,  1853  (1986). 

^S.  Luryi,  Superlattices  Microstructures  5,  475  (1989). 

^T.  Tada,  A.  Yamagushi,  T.  Ninomiya,  H.  Uchiki,  T.  Kobayashi,  and  T. 
Yao,  J.  Appl.  Phys.  63,  5491  (1988). 

*S.  Muto,  T.  Inata,  A.  Tacheuchi,  Y.  Sugiyama,  and  T.  Fujii,  Appl. 
Phys.  Lett.  58,  2393  (1991). 

^M.  Nido,  M.  G.  W.  Alexander,  W.  W.  Ruble,  T.  Schweizer,  and  K. 
Kohler,  Appl.  Phys.  Lett.  56,  355  (1990). 

*D.  Y.  Oberli,  J.  Shah,  T.  C.  Damen,  C.  W.  Tu,  T.  Y.  Chang,  D.  A.  B. 
Miller,  J.  E.  Henry,  R.  F.  Kopf,  N.  Sauer,  and  A.  E.  DiGiovanni,  Phys. 
Rev.  B40,  3028  (1989). 

’K.  Leo,  J.  Shah,  E.  O.  Gobel,  T.  C.  Damen,  S.  Schmitt-Rink,  W.  Scha¬ 
fer,  and  K.  Kohler,  Phys.  Rev.  Lett.  66,  201  (1991). 

*N.  Sawaki,  R.  A.  Hopfel,  E.  Gomik,  and  H.  Kano,  Appl.  Phys.  Lett. 
55,  1996  (1989). 

’C.  Tanguy,  B.  Deveaud,  A.  Regreny,  D.  Hulin,  and  A.  Antonetti,  Appl. 
Phys.  Lett.  58,  1283  (1991). 

’®A.  Tackeuchi,  S.  Muto,  T.  Inata,  and  T.  Fujii,  Appl.  Phys.  Lett.  58, 
1670  (1991). 

D.  H.  Levi,  D.  R.  Wake,  M.  V.  Klein,  S.  Kumar,  and  H.  Morko^,  Phys. 
Rev.  B45,  4274  (1992). 

M.  J.  Hafich,  J.  H.  Quigley,  R.  E.  Owens,  G.  Y.  Robinson,  D.  Li,  and 

N.  Otsuka,  Appl.  Phys.  Utt.  54,  2686  (1989). 

’^M.  J.  Hafich,  H.  Y.  Lee,  G.  Y.  Robinson,  D.  Li,  and  N.  Otsuka,  J.  Appl. 
Phys.  69,  752  (1991). 

‘^M.  O.  Watanabe  and  Y.  Ohba,  Appl.  Phys.  Lett.  50,  906  (1987). 

‘^M.  Zachau  and  W.  T.  Masselink,  Appl.  Phys.  Lett.  60,  2098  (1992). 
*^S.  Mukhopadhyay  and  B.  R.  Nag,  Appl.  Phys.  Lett.  60,  2897  (1992). 
'^R.  Ferreira  and  G.  Bastard,  Phys.  Rev.  B  40.  1074  ( 1989). 


Buccafusca  et  al.  401 


29 


30 


IEEE  PHOTONICS  TECHNOLOGY  LETTERS,  VOL.  5,  NO.  9,  SEPTEMBER  1993 


987 


Strained  Layer  In^Gaj_^As/GaAs  and 
In  ^Ga  1  As/In  ^Ga  1  P 
Multiple-Quantum-Well  Optical  Modulators 
Grown  by  Gas-Source  MBE 

J.  W.  Kim,  C.  W.  Chen,  T.  J.  Vogt,  L.  M.  Woods,  G.  Y.  Robinson,  and  D.  L.  Lile 


Abstract — The  first  results  on  low-power  p-i-n  diode  modula¬ 
tor  structures  using  strained  multiple  quantum  wells  (MQWs) 
of  InGaAs /InGaP  grown  by  gas-source  molecular  beam  epitaxy 
(MBE)  on  GaAs  are  presented.  A  comparison  of  transmission, 
reflection,  and  photocurrent  spectra  for  these  nonresonant  de¬ 
vices  with  those  fabricated  from  InGaAs  /GaAs  indicates  larger 
modulation,  with  a  maximum  change  in  reflection  of  >  42% 
observed  at  5-V  bias  at  a  wavelength  of  0.96  |xm. 


I.  Introduction 

ITH  advances  in  crystal  growth  techniques  such  as 
molecular  beam  epitaxy  (MBE)  and  metalorganic 
deposition  (MOCVD),  multiple-quantum-well  (MQW) 
structures  based  on  a  variety  of  III-V  semiconductor 
materials  have  attracted  considerable  attention  for  their 
potential  use  in  high-speed  optoelectronic  devices.  In  par¬ 
ticular,  and  starting  with  the  work  of  Wood  et  al  in 
GaAs/AlGaAs,  [1]  many  groups  have  used  MQW  struc¬ 
tures  to  investigate  optical  modulation  based  on  the  quan¬ 
tum  confined  Stark  effect  (QCSE).  To  date  most  of  this 
research  has  relied  upon  GaAs/AlGaAs  MQW  p-i-n  diode 
structures,  lattice  matched  on  GaAs  substrates,  which 
operate  at  wavelengths  around  0.85  ^tm  [2].  Since  the 
GaAs  is  opaque  at  these  operating  wavelengths,  cascad- 
ability  for  optical  interconnects,  switching,  or  signal  pro¬ 
cessing  applications  requires  removal  of  the  substrate. 
Ino.53Gao47As/InP  structures,  operating  at  1.5  (im, 
avoid  this  problem,  but  at  the  cost  of  reduced  constrast 
[3].  Another  possibility  is  strained-layer  In^Ga^.^As/GaAs 
QWs  on  GaAs  [4],  which,  in  addition  to  having  a  transpar¬ 
ent  substrate,  also  possess  a  relatively  low  potential  bar¬ 
rier  in  the  conduction  and  valence  bands,  with  the  result¬ 
ing  potential  for  improved  speed  and  reduced  saturation 
via  faster  carrier  sweep-out  times  and  smaller  voltage, 
lower  power  operation  [5].  Similar  advantages  would  be 
expected  to  attend  the  In^Gai_^As/In^Gai_^P  system 
except  that,  in  addition,  and  like  In^Ga^.^As/AlGaAs 
[6],  the  barrier  composition,  strain,  and  related  barrier 
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height  can  also  be  varied.  Strain  compensation  can  also  be 
provided  through  the  use  of  the  ternary  barrier — as 
demonstrated  by  Cunningham  et  al.  [7]  in  InGaAs/GaAsP 
structures — thereby  permitting  the  growth  of  thicker 
modulators  with  a  minimum  of  strain  induced  defects. 
Lasers  have  already  been  demonstrated  in  these  materials 
[8]  and,  operating  in  the  0.95-  to  l.l-/xm  range,  these 
strained  materials  combinations  are  compatible  with  oper¬ 
ation  at  wavelengths  at  which  Si  is  absorbing,  hence 
allowing  for  optical  interconnects  in  hybrid  Ill-V/Si  sys¬ 
tems  [9]. 

In  this  paper  we  present  results  obtained  on  strained- 
layer  MQW  modulators  made  from  In^Ga^^^As/GaAs 
and  In^Gaj.^As/In^Gai^^P  grown  by  gas  source  MBE. 
These  are  the  first  results  for  the  In^Gaj^^As/In^Gaj.^P 
materials  system  by  any  growth  technique  and  suggest 
that,  in  addition  to  providing  greater  material  flexibility, 
the  ternary  barrier  may  lead  to  an  improvement  in  modu¬ 
lator  performance. 

II.  Device  Structure 

The  Inj.Gai.^As/GaAs  and  In^Gai^^As/In3,Gai^^P 
MQW  p-i-n  structures  were  grown  by  gas  source  MBE  at 
500°C  on  n"^  (100)  GaAs  substrates  [10].  Two  samples, 
one  for  each  material  system,  were  grown  with  essentially 
the  same  structure.  This  consisted  of  a  0.5-^tm  n"^  GaAs 
buffer  layer  followed  by  an  n"^  InGaP  etch-stop  layer,  and 
24.5  periods  of  undoped  strained  In^Gai.^^As  quantum 
wells,  separated  by  150-A-thick  barriers  of  either  GaAs 
or  InGaP,  The  structure  was  completed  by  growing  an 
undoped  GaAs  or  InGaP  buffer  layer  and  a  p"^  GaAs  or 
InGaP  cap.  The  nominal  thickness  and  In  mole  fraction  x 
of  the  In^Gaj.jpAs  wells  were  100  A  and  0.15,  respec¬ 
tively,  and  the  In  mole  fraction  y  of  the  In^Gaj.^P 
barriers  was  0.5.  Individual  devices  were  fabricated  in 
both  samples  by  etching  mesas  of  size  360  ptm  X  510  /xm 
and  e-beam  evaporating  Au/Ti  to  make  ohmic  contact  to 
the  n-  and  p-type  GaAs  and  InGaP  regions.  The  resulting 
devices  exhibited  good  electrical  behavior  with  breakdown 
voltages  >  20  V  for  both  types  of  samples.  For  optical 
measurements  the  back  side  of  the  devices  was  chemome- 
chanically  polished,  and  then  both  surfaces  were  antire¬ 
flection  coated  using  Si02.  For  reflection  measurements 
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the  SiOj  on  the  rear  surface  was  replaced  by  an  evapo¬ 
rated  Au  layer. 

III.  Results  and  Discussion 

Photocurrent,  transmission,  and  reflection  spectra  of 
the  individual  modulators  were  measured  using  a 
monochromator  and  conventional  lock-in  techniques.  Fig. 
1  shows  absorption  spectra  versus  reverse  bias  for  an 
In^Gaj.^As/GaAs  device.  Photocurrent  spectra  were 
veiy  similar  in  shape.  The  absorption  coefficient  values 
were  calculated  from  the  transmission  spectra  using  a 
total  QW  layer  thickness  of  0.61  ^tm.  From  these  data  we 
see  that  the  excitonic  resonances  in  this  sample  are  rela¬ 
tively  broad,  resulting  presumably  from  strain-induced 
defects  as  well  as  alloy  disordering.  Visual  inspection 
using  a  Nomarski  microscope  revealed  relatively  weak, 
though  certainly  discemable,  cross  hatchings  across  the 
entire  wafer,  indicating  the  presence  of  misfit  dislocations 
resulting  from  strain  relief.  On  application  of  bias,  the 
shift  in  the  absorption  peaks  to  larger  wavelength  result¬ 
ing  from  the  QCSE  is  clearly  evident.  In  particular,  the 
change  in  the  absorption  coefficient  at  the  zero-bias  heavy 
hole  peak  at  9880  A  is  2743  cm"^  for  7  volts  of  reverse 
bias  corresponding  to  a  contrast  ratio  (CR)  of  1.16.  Due 
to  the  small  difference  in  the  bandgap  energy  of  the 
InGaAs  well  material,  estimated  in  our  case  to  be  1.27 
eV,  and  the  GaAs  (1.42  eV)  barrier,  we  would  expect 
rapid  dissociation  of  the  excitons  at  low  applied  voltages 
due  to  field-assisted  tunneling.  This  is  evidenced  in  Fig.  1 
by  both  the  rapid  distortion  of  the  spectra  at  low  biases 
and  the  associated  marked  broadening  of  the  excitons. 
Noteworthy  for  low-power  modulator  applications  is  the 
almost  complete  extinguishing  of  the  exciton  resonance  at 
voltages  as  low  as  5  V,  resulting  presumably  from  the 
shallow  wells  in  this  structure  [7]. 

In  Fig.  2  we  show  the  photocurrent  and  equivalent 
absorption  data  to  Fig.  1,  but  for  an  In^Gai_,As/ 
InyGai_yP  device,  where  we  have  normalized  the  pho¬ 
tocurrent  curves  to  the  same  value  at  0.90  fim  to  more 
clearly  show  the  progression  of  the  exciton  peaks  with 
bias.  The  inset  shows  the  fit  between  the  experimental 
and  theoretical  shift  of  the  heavy  hole  peak  with  bias.  The 
theoretical  curve  in  this  figure  was  calculated  using  a 
second-order  perturbation  theory  [11],  using  values  for  the 
well  width  of  90A,  obtained  from  TEM  cross  sections  [12], 
and  effective  mass  values  of  0.0627  and  0.444o  for  the 
electron  and  heavy  hole,  respectively.  These  mass  values 
were  chosen  to  effect  the  best  fit  with  the  data  and  are 
similar  to  those  reported  earlier  by  others  [13],  [14].  It  is 
clear  in  both  the  sets  of  data  in  Fig.  2,  but  most  pro¬ 
nounced  in  the  photocurrent  curves,  that  in  contrast  to 
the  GaAs  barrier  sample,  where  weak  transitions  at  higher 
energies  completely  disappear  with  applied  field,  increas¬ 
ing  reverse  bias  leads  to  the  appearance  of  strong  higher 
order  peaks.  These  are  somewhat  reminiscent  of  the  spec¬ 
tra  observed  by  Pezeshki  et  al.  [6]  on  InGaAs/ AlGaAs , 
and  may  result  in  both  cases  from  strain.  Although  uncer¬ 
tainties  in  barrier  heights  and  strain  in  this  material 
system  make  the  calculations  of  exciton  energies  some¬ 


Fig.  1.  Absorption  spectra  of  an  InGaAs/ GaAs  device. 


(a) 


Fig.  2.  (a)  Photocurrent  spectra  of  an  InGaAs/InGaP  device.  The 
arrows  show  the  calculated  positions  of  the  first  three  exciton  peaks,  (b) 
Absorption  spectra  of  an  InGaAs/InGaP  device.  The  inset  shows  exper¬ 
imental  (dotted)  and  theoretical  (solid)  Stark  shifts  of  the  \ehh  exciton 
versus  reverse  bias. 

what  imprecise,  reasonable  estimates  for  these  quantities 
lead  to  resonance  values  not  inconsistent  with  those  in  the 
data  in  Fig.  2,  as  indicated  by  the  arrows  which  show  the 
calculated  positions  of  the  \ehh,  lelhh,  and  lelh  peaks 
using  the  same  mass  values  as  above.  Presumably  the 
absence  of  higher  energy  excitons  in  the  spectra  of  Fig.  1 
results  from  the  extreme  shallowness  of  the  wells  and 
their  inability  to  confine  other  than  the  n  =  1  level. 
Exciton  full  width  at  half  maximum  (FWHM),  obtained 
from  a  Lorentzian  fit  with  the  lehh  transition  at  zero  bias 
in  Fig.  2(b),  gives  a  value  of  10.5  meV,  slightly  smaller 
than  the  value  of  14  meV  reported  for  the  InGaAs/ Al¬ 
GaAs  MQW  system  [15]. 

Reflection  spectra  and  the  associated  transmission  and 
reflection  modulation  behavior  of  the  InGaP  barrier  de¬ 
vice  are  shown  in  Figs.  3-5.  From  the  last  two  figures  we 
see  that  the  majority  of  the  modulation  takes  place  for  the 
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Wavelength  ( jum  ) 

Fig.  3.  Reflection  spectra  of  an  InGaAs/InGaP  MQW  modulator 
measured  at  various  values  of  reverse  bias. 


Wavelength  ( /im  ) 

Fig.  4.  Transmission  modulation  spectra  of  an  InGaAs/InGaP  MQW 

device. 


Wavelength  ( Aim ) 

Fig.  5.  Reflection  modulation  spectra  of  an  InGaAs/InGaP  MQW 

device. 


first  7  volts  of  applied  bias.  From  the  transmission  data, 
and  at  the  peak  wavelength  of  0.959  ptm,  we  obtain  an 
absorption  coefficient  change  of  3511  cm"^  correspond¬ 
ing  to  a  9%  change  in  absolute  transmission  and  a  frac¬ 
tional  change  of  24%  for  7-V  bias.  This  gives  a  CR  of 
1.24,  which  should  be  compared  to  the  value  of  1.16 
obtained  under  the  same  conditions  in  the  GaAs  barrier 
structures.  In  reflection,  because  of  the  double  pass 
through  the  optically  active  QW  stack,  we  observe  a  larger 
effect  with  >  42%  change  in  the  optical  signal.  These 
results  suggest  that  InGaP  may  result  in  some  significant 
advantages  over  GaAs  when  employed  with  InGaAs  for 
MQW  modulator  applications.  Certainly,  the  variable 
composition  of  barrier  and  well  means  that  band  offsets 
and  strain  can  be  varied  independently  to  optimize  the 
electrooptic  effect.  From  our  results  it  is  evident  that  the 


increased  barrier  height,  compared  to  that  with  GaAs,  has 
enhanced  the  excitonic  effects  at  zero  bias  and  increased 
the  resistance  to  rapid  deterioration  of  the  excitonic  reso¬ 
nance  with  increasing  applied  biases.  This  is  similar  to 
results  achieved  by  adding  aluminum  to  the  barrier  re¬ 
gions  in  the  InGaAs/AlGaAs  material  system  to  enhance 
electroabsorption  [6].  Another  advantage  is  the  flexibility 
the  ternary  barrier  allows  for  the  growth  of  a  sample  with 
a  balanced  strain,  where  biaxial  compressive  strain  in¬ 
duced  in  the  InGaAs  well  due  to  a  larger  lattice  constant 
can  be  balanced  with  a  tensile  strain  occurring  in  the 
InGaP  barrier  [7],  [16].  Thus,  in  principle,  a  sample  with 
InGaP  barriers  (y  <  0.48)  and  InGaAs  wells  can  be  grown 
with  no  net  strain  accumulation,  resulting  in  no  lattice 
relaxation  and  a  low  density  of  dislocations.  In  this  way  a 
greater  number  of  wells  can  be  grown  so  as  to  increase 
the  optical  interaction  distance  and  thereby  enhance  the 
modulation. 
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We  present  an  observation  of  the  Wannier-Stark  effect  in  a  strained  InGaAs/InGaP  superlattice 
grown  on  a  GaAs  substrate.  A  blueshift  of  the  effective  absorption  edge  is  observed  in  room  and 
low-temperature  photocurrent  and  transmission  measurements.  A  ~2000  cm“'  absorption  change 
due  to  the  transition  of  the  absorption  edge  from  a  broad  to  a  sharp  quantum  well-like  excitonic 
shape  was  obtained  for  as  little  as  a  2  V  change  in  bias  voltage.  ©  1996  American  Institute  of 
Physics.  [80021-8979(96)02409-2] 


INTRODUCTION 

In  recent  years,  there  has  been  growing  interest  in  elec- 
troabsorption  modulators  based  on  Wannier-Stark  localiza¬ 
tion  in  superlattice  (SL)  structures  because  of  their  potential 
for  low  voltage,  and,  hence,  low  power,  operation  in  optical 
interconnect  and  switching  systems  applications.'  In  contrast 
to  a  quantum  well  (QW)  structure,  a  superlattice  consists  of 
alternating  semiconductor  layers  where  the  barrier  layer  is 
sufficiently  thin  to  allow  appreciable  tunneling  and  a  strong 
resonant  coupling  between  adjacent  quantum  wells.  Due  to 
this  coupling,  the  discrete  energy  levels  found  in  QWs 
broaden  into  minibands  of  widths  iSE  and  LH  in  the  con¬ 
duction  and  valence  bands,  respectively.  By  applying  an 
electric  field  the  resonant  condition  of  aligned  energy  levels 
in  adjacent  wells  can  be  easily  destroyed  so  that  the  initially 
delocalized  SL  wave  functions  become  localized  over  a  few 
adjacent  QWs,  leading  to  a  blueshift  of  the  effective  absorp¬ 
tion  edge  of  the  order  (A£+A/7)/2.^  Since  this  Wannier- 
Stark  localization  can  occur  at  a  relatively  small  electric 
field,  typically  in  the  range  of  a  few  tens  of  kV/cm,  a  low 
drive  voltage  electroabsorption  modulator  can  in  principle  be 
realized  using  the  SL  structure.  To  date  the  Wannier-Stark 
effect  has  been  mostly  observed  in  the  GaAs/AlGaAs  and 
InP-based  material  systems  lattice  matched  on  GaAs  and  InP 
substrates,  respectively.^*’  A  22  000  cm"'  absorption  change 
for  6.8  V  voltage  swing  in  a  GaAs/AlGaAs  SL^  and  a  reflec¬ 
tion  change  of  32%  for  4  V  of  bias  in  an  InGaAs/InP  SL 
modulator’  are  examples  of  results  reported  to  date.  More 
recently,  with  advances  in  crystal  growth  techniques, 
strained-layer  InGaAs/GaAs  SLs  have  also  been 
fabricated.*’®  In  this  paper,  we  report  for  the  first  time  obser¬ 
vation  of  a  blueshift  in  a  superlattice  made  from  strained 
lUj-Gai  _^As/Inj,Gai  a  materials  system  which  has  been 
previously  employed  for  low  power  multiple  QW 
modulators.'® 

EXPERIMENT 

The  In^Gaj  _jAs/InyGaj _yP  SL  p-i-n  structures  used 
in  this  study  were  grown  by  gas  source  molecular  beam  ep¬ 
itaxy  (MBE)  on  n'''(001)  G^s  substrates  at  500  “C.  The 
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sample  consisted  of  a  0.6  ^tm  thick  GaAs  buffer  layer 
followed  by  an  undoped  300  A  GaAs  layer,  and  49.5  periods 
of  an  undoped  60  A  Ino.ijGao.gyAs  well/20  A  Ino.48Gao.52P 
barrier  SL  of  residual  carrier  concentration  /i~l6'®  cm"^. 
The  structure  was  completed  by  growing  an  undoped  300  A 
thick  GaAs  layer  and  a  /?'*’  GaAs  cap  layer  as  shown  sche¬ 
matically  in  Fig.  1 .  Individual  devices  were  fabricated  by  wet 
etching  mesas  of  size  360X510  /im’  Ohmic  contacts  were 
made  by  e-beam  evaporating  Ti/Au  on  both  the  n-  and 
p-type  GaAs  regions.  Typical  I-V  characteristics  of  the 
p-i-n  diodes  showed  low  reverse  leakage  currents  of  <100 
nA  and  breakdown  voltages  >15  V.  For  optical  measure¬ 
ments  the  back  side  of  the  substrate  was  polished  and  anti- 
reflection  coated  using  Si02.  Photocurrent  and  transmission 
measurements  were  then  performed  at  room  and  low  tem¬ 
peratures  (~20  K)  using  a  monochromator  and  conventional 
lock-in  techniques. 


RESULTS 

Figure  2  shows  typical  photocurrent  (PC)  spectra  ob¬ 
tained  as  a  function  of  reverse  bias  voltage  at  300  K  and  20 
K.  At  low  bias,  the  spectra  show  no  significant  excitonic 
features,  consistent  with  the  presence  of  minibands  and  de- 
localization  of  the  earners.  In  Fig.  2(a)  we  observe  two 
bumps,  identified  by  the  arrows,  near  1.319  eV  and  1.366 
eV,  which  are  presumably  associated  with  miniband  transi¬ 
tions  from  the  ground-state  electron  to  heavy  hole  (HH)  and 
light  hole  (LH)  excited  states,  respectively.  For  a  quantitative 
analysis,  the  SL  band  structure  was  calculated  using  a  simple 
Kronig-Penney  model.  Values  of  the  material  parameters 
were  estimated  by  linearly  interpolating  between  the  binary 
parent  compounds."  The  band  offsets  at  the  heterojunction 
in  the  presence  of  strain  were  determined  using  the  model- 
solid  theory  of  Van  de  Walle.'’  Although  x-ray  diffraction 
indicates  dislocations  are  distributed  throughout  the  superlat¬ 
tice,  for  purposes  of  calculation,  we  assumed  that  misfit  dis¬ 
locations  were  confined  only  to  the  interface  between  the 
multilayer  and  the  substrate  and  thus  that  the  superlattice 
tends  to  the  average  lattice  constant  of  the  unsupported 
multilayer  structure.'*  The  calculated  band-offset  values  are 
197,  418,  and  352  meV  for  the  electron,  HH,  and  LH,  re- 
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^ _ p^-  GaAs  (SOOnm) 

^ —  i  -  GaAs  (30nm) 

^  SL  InGaAs/InGaP 

(6nm/2nm) 

4 _ i  -  GaAs  (30nm) 

^ _ ri*"-  GaAs  buffer 

(600nm) 

n  -  GaAs  substrate 

FIG.  1.  Schematic  diagram  of  the  strained  InGaAs/InGaP  superlattice  struc¬ 
ture. 

spectively.  Corresponding  fractional  conduction  band  offsets 
for  the  electron-heavy-hole  and  the  electron-light-hole  sys- 
tem  are  0.32  and  0.36,  respectively.  With  these  values,  the 
carrier  miniband  widths  were  computed;  35,  0.1,  and  16 
meV  for  the  electron,  HH,  and  LH,  respectively.  As  the  ap¬ 
plied  reverse  bias  increases,  the  absoiption  edge  evolves 
from  that  appropriate  to  a  3D  SL  to  a  2D  QW-like  excitonic 
line  shape  due  to  the  Wannier— Stark  localization.  A  rela¬ 
tively  small  electric  field  induces  carrier  wave  function  lo¬ 
calization  and  thus  increases  the  excitonic  binding  energy 
and  the  excitonic  oscillator  strength.  This  can  be  seen  in  Fig. 
2(a)  where  a  well-resolved  excitonic  peak  appears  at  1.310 
eV,  in  good  agreement  with  a  calculated  value  of  1.308  eV 


Energy  (eV) 


Energy  (eV) 

FIG.  2.  Photocurrent  spectra  taken  at  (a)  300  K  and  (b)  20  K  for  various 
values  of  reverse  bias. 


Energy  (eV) 


FIG.  3.  Low  temperature  absorption  spectra  of  sample  in  the  presence  of 
various  applied  biases. 

we  obtained  using  a  QW  exciton  binding  energy  of  8  meV.^ 
The  20  K  PC  spectra  shown  in  Fig.  2(b)  are  similar  to  those 
at  300  K,  but  show  more  clearly  a  blueshift  of  the  effective 
absorption  edge  at  ~  1.364  eV.  In  addition,  a  weak  absorp¬ 
tion  variation  is  also  observed  below  the  SL  band  edge  in 
Fig.  2(b),  possibly  caused  by  the  appearance  of  the  low- 
energy  oblique  transitions  due  to  the  carriers  localized  in 
different  wells.  Nevertheless,  in  this  sample  a  pronounced 
Stark  ladder  was  not  observed,  possibly  as  a  result  of  gen¬ 
eration  of  misfit  dislocations  due  to  strain  relaxation  and  the 
electric  field  inhomogeneity  due  to  the  large  residual  doping. 

Transmission  measurements  were  performed  to  directly 
study  electroabsorption  properties  of  the  SLs.  Figure  3  shows 
absorption  spectra  at  20  K  in  the  presence  of  various  applied 
biases.  The  absorption  coefficient  values  were  calculated 
from  the  transmission  spectra  using  a  total  SL  layer  thickness 
of  0.4  fim.  We  obtain  a  positive  absorption  coefficient 
change  of  -2000  cm~“  at  the  QW-like  excitonic  peak  at 
1.393  eV  for  only  2  V  bias  swing  and,  due  to  the  blueshift  of 
the  effective  absorption  edge,  a  negative  change  of  1565 
cm~*  at  wavelengths  close  to  1.375  eV.  At  high  reverse  bi¬ 
ases  S3  V,  the  oscillator  strength  of  the  QW  exciton  tends  to 
saturate,  corresponding  to  full  localization  of  the  carriers, 
and  the  exciton  peak  starts  to  redshift  due  to  the  quantum 
confined  Stark  effect  (QCSE),  giving  rise  to  an  absorption 
change  of  2770  cm"‘  for  applied  voltages  of  - 1  V  to  -7V 
at  1.393  eV.  The  field-induced  Wannier-Stark  shifts  of  the 
SL,  corresponding  to  a  blueshift  of  the  absorption  edge  can 


FIG.  4.  Transmission  modulation  spectra  measured  at  20  K.  The  relative 
transmissions  were  calculated  with  respect  to  the  transmission  Tq  at  1  V. 
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be  used  in  a  “normally  off’  electroabsorptive  modulator  to 
achieve  high  modulation  ratios  at  moderately  low  bias  volt¬ 
ages.  A  contrast  ratio  greater  than  130:1  from  this  mecha¬ 
nism  was  reported  using  GaAs/AlGaAs  SL  asymmetric 
Fabry-Perot  modulators.*"^ 

Figure  4  shows  the  transmission  modulation  response 
(Ar/fo)  where  Tq  is  the  transmission  at  +1  V.  From  the 
data  we  see  that  the  majority  of  the  modulation  occurs  for  the 
first  2  volts  of  applied  bias.  Maximum  fractional  modulation 
at  -2  V  is  ^8%  at  1.375  eV.  Our  results  suggest  that  the 
InGaAs/InGaP  SL  can  have  application  as  a  low-drive  volt¬ 
age  optical  modulator.  The  relatively  low  values  of  modula¬ 
tion  could  presumably  be  enhanced,  as  has  been  accom¬ 
plished  in  other  materials  systems,  using  reflection  or 
resonant  structures.  In  addition,  the  use  of  the  ternary  binary 
InGaP  provides  the  flexibility  to  tailor  the  electronic  proper¬ 
ties  of  the  SL  layer  through  changes  in  the  barrier  composi¬ 
tion,  and  further  allows  for  the  growth  of  samples  with  bal¬ 
anced  strain  in  both  the  well  and  barrier,  possibly  resulting  in 
an  enhancement  in  modulator  performance. 
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Abstract — We  report  modulation  saturation  and  time  response 
measurements  on  InGaAs-InGaP  MQW  modulators.  The  mea¬ 
surements  yield  a  saturation  intensity  of  (3.7  ±  0.1)  kW/cm- 
for  a  0-10  V  swing  and  switching  times  between  10  and  90 
ns,  depending  on  the  bias  voltage  and  incident  light  intensity. 
The  observed  dependence  indicates  that  field  screening  due  to 
carrier  build-up  is  the  dominant  physical  mechanism  determining 
both  the  speed  and  the  saturation  intensity.  This  conclusion  is 
supported  by  results  of  theoretical  calculations. 

L  Introduction 

INCREASING  demands  for  high-speed  optical  switching 
devices  has  sparked  considerable  interest  in  the  area  of 
addressable  optical  modulators.  Much  success  has  come  from 
work  on  multiple  quantum  well  modulators  utilizing  the  quan¬ 
tum  confined  Stark  effect  (QCSE).  Material  systems  which 
have  been  studied  for  this  purpose  include  GaAs-AIGaAs 
[1H3]  InGaAs-AlGaAs  [4],  [5],  InGaAs-GaAs  [6H8], 
and  InGaAs-GaAsP  [9].  Recently,  p-i-n  modulators  based 
on  strained  layers  of  InGaAs-InGaP  grown  on  GaAs  have 
also  been  demonstrated  [10].  The  InGaAs— InGaP  modulators, 
which  offer  the  advantage  of  a  transparent  GaAs  substrate, 
were  shown  to  require  low  bias  voltages  (~  5  V)  for  operation, 
and  hence  reduced  power  dissipation.  Modulation  of  the 
absorption  at  the  wavelength  of  the  excitonic  peak  up  to  25% 
has  been  obtained  for  a  single  transmission  pass.  However,  two 
important  device  parameters  for  its  utilization  in  applications: 
the  saturation  intensity  and  time  response  had  not  been  studied. 

In  this  paper,  we  report  the  results  of  a  study  of  saturation 
intensity  and  time  response  for  nonresonant  InGaAs-InGaP 
MQW  modulators  and  analyze  the  dominant  physical  mech¬ 
anisms  responsible  for  the  observed  device  performance.  The 
MQW  modulators  used  in  the  experiments  are  p-i-n  diodes 
with  the  intrinsic  region  comprised  of  24.5  periods  of  undoped, 
strain  compensated  Ina,Gai_;cAs-InyGai..yP  quantum  wells 
fabricated  using  gas  source  MBE  [10].  The  mole  fractions 
used  were  a:  =  0.15  and  y  =  0.5.  The  top  and  bottom  of  the 
modulators  have  InGaP  buffer  layers  0.5  thick.  The  quan¬ 
tum  well  thickness  is  10  nm  and  the  barrier  thickness  is  15  nm. 
The  barrier  thickness  was  chosen  to  provide  adequate  tensile 
strain  to  compensate  the  compressive  strain  of  the  mismatched 
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well  layers,  and  thus  avoid  the  formation  of  misfit  dislocations. 
The  devices  were  grown  with  360  /zm  x  510  /zm  active  areas, 
Au-Ti  contacts,  and  Si02  anti-refiection  coatings.  The  samples 
used  were  from  an  array  of  modulators  on  the  same  region  of 
a  wafer  which  was  configured  for  transmission  measurements. 

The  next  section  discusses  the  saturation  intensity  mea¬ 
surements  and  Section  III  the  results  of  the  time  re,sponse 
measurements.  Finally,  Section  IV  presents  a  simple  model  for 
the  carrier  sweep-out  that  describes  well  the  observed  intensity 
and  voltage  dependence  of  the  transmission  of  the  devices  and 
gives  values  of  the  saturation  intensity  and  time  response  of  the 
modulation  which  are  similar  to  those  obtained  experimentally. 

II.  Measurements  of  the  Saturation  Intensity 

The  performance  of  modulators  at  high  intensity  is  affected 
by  several  phenomena  caused  by  the  high  density  of  pho¬ 
togenerated  carriers  including  phase  space  filling  [3],  field 
screening  due  to  charge  build-up  [1 1],  [12],  and  thermal  effects 
[13].  To  determine  the  saturation  intensity,  we  measured  the 
transmission  of  the  modulator  at  the  excitonic  wavelength  as 
a  function  of  incident  light  intensity.  For  this  measurement 
we  used  a  tunable  CW  TirSapphire  laser  pumped  by  a  CW 
argon  ion  laser  as  the  light  source.  The  light  detection  was 
done  using  p-i-n  silicon  photodiodes  and  standard  lock-in 
detection  techniques.  To  minimize  heating  of  the  sample,  the 
output  of  the  laser  was  modulated  using  a  specially  constructed 
optical  chopper  to  achieve  pulses  of  ~4  ^s  duration  with  a 
repetition  rate  of  ~3.5  kHz.  The  beam  was  focused  on  the 
sample  utilizing  a  5  cm  focal  length  lens  and  the  sample  was 
positioned  at  the  focal  plane.  The  beam  profile  was  measured 
to  be  approximately  Gaussian  and  circular,  with  a  1/e  intensity 
radius  of  10  ^zm  at  the  focus.  The  power  density  on  the  sample 
was  adjusted  by  selecting  the  power  of  the  laser  using  a 
variable  neutral  density  filter. 

Fig.  I  shows  the  transmission  spectra  of  the  InGaAs-InGaP 
modulators  for  applied  voltages  of  0  and  6  V  of  reverse  bias 
for  two  different  optical  powers.  For  low  intensity  (75  W/cm-) 
and  no  bias  applied  to  the  device,  the  excitonic  resonance  is 
clearly  visible  as  an  absorption  peak  at  960  nm.  Applying  6 
V  to  the  device,  the  excitonic  peak  shifts  and  broadens  due  to 
the  quantum  confined  Stark  effect,  resulting  in  a  net  change  in 
transmission  at  the  excitonic  wavelength  that  allows  the  use 
of  this  device  as  a  modulator  at  that  wavelength.  With  a  10  V 
swing  in  applied  bias,  the  contrast  ratio  is  1.25  which  could 
be  increased  by  configuring  the  device  for  reflection  operation 
[10].  Fig.  1(b)  shows  the  effect  of  high  intensity  (2.5  kW/cm^) 
on  the  transmission  spectrum  of  the  device.  The  0  V  data  still 
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Fig.  I.  Transmission  spectra  of  the  modulator  measured  at  two  different 
intensities.  75  W/cm-  (a)  and  2.5  kW/cm'  (b),  and  applied  voltages,  0  (solid 
line)  and  6  V  (dashes). 


shows  the  excitonic  peak,  although  not  as  distinct  as  in  the 
low  intensity  case.  The  change  in  transmission  with  applied 
voltage  is  also  smaller  at  high  intensity,  thus  reducing  the 
contrast  ratio  of  the  device. 

From  a  technological  point  of  view,  and  regardless  of  the 
physical  origin  of  this  saturation,  it  is  necessary  to  define  a 
parameter  that  characterizes  the  dependence  of  the  contrast 
ratio  CR  with  incident  intensity  I.  As  done  in  [14],  we  will 
define  the  modulation  saturation  intensity  /gat  according  to 
the  following  expression 


CR(/)  =  1  + 


CR(0)  -  1 

^  +  1/ -fsat 


(1) 


The  measured  dependence  of  the  contrast  ratio  with  incident 
power  density  was  fitted  with  expression  (I)  for  voltage  swings 
ranging  from  (0-3  V)  to  (O-IO  V).  Fig.  2  shows  the  change 
in  contrast  ratio  corresponding  to  a  0-10  V  bias  swing  for 
four  different  devices  for  input  intensities  between  0. 1  and  60 
kW/cm".  As  can  be  observed  in  the  figure,  the  dependence  of 
the  contrast  ratio  as  a  function  of  the  light  intensity  is  well 
described  by  expression  (1)  (solid  line)  for  all  four  devices. 
The  fitted  value  of  saturation  intensity  for  a  0-10  V  swing  is 
(3.7  ±0.1)  kW/cm^. 

Fig.  3  shows  the  resulting  dependence  of  modulation  and 
saturation  intensity  with  voltage.  At  low  power,  the  modulation 
saturates  at  about  6  V  (Fig.  3(a)).  This  is  consistent  with  the 
fact  that  the  excitonic  absorption  is  washed  out  for  voltages 
greater  than  ~6  V  [lOj.  For  higher  power,  the  modulation 
remains  linear  with  voltage  within  the  0-10  V  range  studied, 
but  its  maximum  decreases  due  to  saturation.  The  saturation 


Fig.  2.  Contrast  ratio  versus  incident  intensity  for  a  10  V  change  in  applied 
voltage.  The  symbols  represent  different  devices  in  the  same  wafer,  the  solid 
curve  is  a  least  squares  fitting. 


Fig.  3.  Dependence  of  modulation  (a)  and  saturation  intensity  (b)  as  a 
function  of  operating  voltage. 


intensity  as  a  function  of  applied  voltage  is  shown  in  Fig.  3(b). 
The  saturation  intensity  is  relatively  constant  with  a  value  of 
(1.3  ±  0.1)  kW/cm^  up  to  ~6  V  where  the  saturation  begins  to 
increase  linearly  with  voltage  to  reach  (3.7  ±  0.1)  kW/cm“  at 
10  V.  This  saturation  intensity  value  is  relatively  small  when 
compared  with  that  of  similar  devices  constructed  utilizing 
other  material  systems  [3],  [14],  [15]. 

A  physical  mechanism  affecting  the  performance  of  mod¬ 
ulators  at  high  light  intensities  is  the  saturation  of  the  ex¬ 
citonic  absorption  [3].  Our  measurements  indicate  that  in 
these  devices,  this  is  not  the  dominant  physical  mechanism 
responsible  for  the  observed  saturation  of  the  modulation.  This 
can  be  concluded  from  Fig.  4.  Fig.  4(a)  shows  the  differential 
transmission  spectra  for  a  0-6  V  swing  both  at  low  and  high 
intensity.  The  operating  wavelength  is  indicated  by  a  vertical 
dotted  line  in  the  figure.  It  is  evident  that  degradation  of 
the  modulation  occurs  at  high  intensity.  Fig.  4(b)  shows  the 
corresponding  transmission  spectra  at  0  and  6  V  for  low  and 
high  intensity.  The  transmission  at  the  operation  wavelength 
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Fig.  4.  (a)  Differential  transmission  spectra  at  low  and  high  intensities  for 
a  0--6  V  voltage  swing,  (b)  Transmission  spectra  at  0  and  6  V  for  the  same 
intensities  as  in  (a). 

drops  noticeably  at  high  intensity  for  the  6  V  data,  contrary 
to  what  would  be  produced  by  saturation  of  the  excitonic 
absorption.  We  note  that  most  of  the  saturation  observed  in 
the  modulation  is  due  to  the  change  in  transmission  at  6  V, 
suggesting  that  field  screening  due  to  charge  build-up  [11], 
[12]  and/or  thermal  effects  caused  by  ohmic  heating  associated 
with  the  photocurrent  [13]  are  responsible  for  the  observed 
saturation  of  the  modulation.  The  role  of  these  mechanisms 
will  be  discussed  in  the  following  sections. 

As  discussed  in  [3],  there  is  a  close  relationship  between  the 
sweep  out  time  of  the  carriers  and  the  value  of  the  saturation 
intensity.  For  that  reason,  time  response  measurements  on  the 
modulators  not  only  have  a  value  of  their  own  in  characterizing 
device  performance,  but  also  provide  further  evidence  that 
allows  the  attribution  of  the  observed  saturation  to  field 
screening  effects.  The  results  of  time  response  measurements 
are  discussed  in  the  following  section. 

III.  Time  Response  Measurements 

The  time  response  of  the  device  was  characterized  by  mea¬ 
suring  the  change  in  transmission  at  the  excitonic  wavelength 
produced  by  fast  (<4  ns  risetime)  electrical  bias  pulses  applied 
to  the  device.  Simultaneous  measurement  of  the  current  shows 
that  the  capacitance  of  the  device  ('^30“50  pF)  is  not  a 
limitation  for  its  speed  in  this  range.  The  transmitted  light 
was  detected  with  a  silicon  photodetector  having  a  response 
time  of  approximately  1  ns.  The  signal  was  amplified  by  a  fast 
(3  ns  risetime)  amplifier  and  was  measured  with  a  250  MHz 


Fig.  5.  Switching  lime  (I09f-909^-)  versus  operating  voltage  measured  for 
50  W/cm“  and  2.5  kW/cm'.  The  lines  are  only  intended  to  guide  the  eye. 

oscilloscope.  The  overall  resolution  of  the  system  thus  allows 
us  to  measure  changes  in  transmission  occurring  in  times  of 
approximately  4  ns. 

We  measured  the  switching  of  the  transmission  for  different 
voltages  and  incident  light  intensities.  The  results  are  summa¬ 
rized  in  Fig.  5,  where  we  plot  the  switching  times  (for  10%  to 
90%  change  in  transmission)  for  applied  voltages  from  2  to  10 
V,  at  low  and  high  incident  intensities  of  50  W/cm-  and  2.5 
kW/cm“,  respectively.  For  low  intensity  there  is  a  significant 
reduction  in  the  switching  time  from  90  ns  at  2  V  to  10  ns  at 
10  V  as  applied  bias  increases.  In  contrast,  for  high  intensity 
the  switching  time  remains  relatively  constant  for  all  applied 
voltages.  As  discussed  below,  the  ob.served  dependence  with 
voltage  and  intensity  is  consistent  with  the  presence  of  field 
screening  and  rules  out  heating  as  the  dominant  mechanism. 
Let  us  note  here  that  the  measured  switching  times  are  slower 
than  other  devices  previously  reported  in  the  literature  [1], 
[11],  which  is  consistent  with  a  lower  saturation  intensity  [3]. 

Heating  effects  were  studied  for  AlGaAs/GaAs  p-i-n  mod¬ 
ulators  [13].  As  thermal  properties  are  similar  for  all  III-V 
materials,  we  would  expect  comparable  results  for  a  similar 
structure  in  the  InGaP/lnGaAs  material  system.  As  shown 
in  [13],  when  heating  is  the  dominant  mechanism  the  time 
response  of  the  transmission  has  two  components,  one  fast 
corresponding  to  the  onset  of  the  applied  field,  and  one  slow 
corresponding  to  the  change  in  transmission  due  to  heating. 
As  heat  is  produced  mainly  by  ohmic  dissipation  caused  by 
the  photocurrent  [13],  the  power  to  the  device  increases  with 
applied  voltage,  making  the  slow  component  of  the  signal 
more  noticeable  as  voltage  increases.  The  decay  time  of  the 
slow  component  was  determined  [13]  to  be  close  to  50  ns  and 
should  be  independent  of  the  applied  voltage. 

Although  the  measured  switching  times  for  our  devices  are 
in  that  range,  our  experiments  do  not  show  any  of  the  features 
that  can  be  attributed  to  heating,  as  can  be  observed  in  Fig.  6, 
which  shows  the  temporal  evolution  of  the  transmission  of 
the  device  for  different  conditions  of  incident  intensity  and 
bias  voltage.  All  curves  have  been  normalized  to  the  steady 
state  modulation  for  ease  of  comparison.  Curves  a  and  b  were 
measured  with  50  W/cm“  of  incident  intensity,  while  curves  c 
and  d  correspond  to  2.5  kW/cm“.  The  applied  voltage  changes 
from  0  to  10  V  (a,c)  and  0  to  4  V  (b,d).  The  two  component 
rise  of  the  transmission  is  not  observed  in  any  of  the  curves, 
while  the  change  in  switching  speed  is  evident. 
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Fig.  6.  Measured  temporal  evolution  of  the  transmission  of  the  device 
normalized  to  the  steady  state  modulation.  The  applied  voltage  changes  from 
0  to  the  voltage  in  the  figure  with  a  risetime  of  4  ns,  the  intensity  of  the  laser 
is  also  indicated  in  the  figure. 


Fig.  4  provides  further  evidence  against  the  hypothesis  of 
thermal  effects  because  heat  absorbed  in  the  sample  would  pro¬ 
duce  a  localized  increase  in  lattice  temperature  which  would 
shift  the  band  gap  and  accompanying  excitonic  resonance 
towards  lower  energy.  This  shift  would  move  the  excitonic 
absorption  out  of  the  operation  wavelength  of  the  device,  thus 
increasing  the  transmission,  contrary  to  the  behavior  of  the 
transmission  shown  in  Fig.  4. 

Field  screening  effects  [1 1],  [12]  occur  due  to  charge  build¬ 
up  of  one  of  the  photogenerated  carriers:  either  electrons  or 
holes.  Carriers  may  escape  from  the  quantum  wells  through 
thermionic  emission  or  tunneling  through  the  barrier.  As  it 
is  likely  that  electrons  escape  much  more  rapidly  than  holes 
due  to  the  difference  in  effective  masses,  accumulation  of 
holes  results  in  charge  build-up.  Tunneling  escape  probabil¬ 
ity  is  negligible  because  the  barriers  are  15  nm  thick  [3], 
leaving  thermionic  emission  as  the  only  escape  process.  The 
thermionic  emission  rate  for  holes  [16]  is  given  by 


1  (  IcbT  y/-  /  H{F)\ 

r~[27rTnnuLlJ  ^bT;' 


(2) 


In  this  expression  T  is  the  lattice  temperature,  rnhw  is  the 
effective  mass  for  holes  in  the  well,  ks  is  Boltzmann  constant 
and  Lu-  is  the  well  width.  The  barrier  height  H  depends  on 
the  field  F  through 


H{F)  =  Q^Eg  -  Eh  -  |e|FL„72  (3) 


where  Q  is  the  fraction  of  discontinuity  in  the  valence  band. 
AFj  is  the  difference  in  band  gaps  between  the  well  and 
barrier  material,  e  is  the  electron  charge  and  Eu  is  the 
energy  of  the  n  =  1  level  for  holes.  Although  measurements 
in  AlGaAs/GaAs  [12]  showed  that  the  escape  time  is  not 
accurately  predicted  by  expression  (2),  it  gives  the  correct 
qualitative  behavior  and  is  a  good  approximation  for  a  simple 


model.  Calculations  a.ssuming  60%  [17]  of  the  ~0.67  eV  of 
bandgap  difference  [18]  occurs  in  the  valence  band  offset 
reveal  that  the  hole  thermionic  lifetime  varies  between  a  few  ns 
and  400  ns  depending  on  the  applied  field.  A  similar  time  scale 
is  to  be  expected  for  changes  in  the  transmission.  From  (2) 
we  can  see  that  reduction  of  barrier  height  (through  increased 
applied  bias)  results  in  faster  escape  times,  in  agreement  with 
the  low  intensity  data  in  Figs.  5  and  6.  Also,  for  high  intensity, 
the  accumulated  carriers  would  be  expected  to  screen  the  field 
so  that  no  change  in  speed  would  result,  again  in  agreement 
with  the  observed  behavior  of  the  InGaAs-lnGaP  modulators. 
Similar  effects  have  been  observed  in  InGaAs-InP  MQW 
modulators  [11],  on  a  much  faster  time  scale,  which  agrees 
with  estimates  from  (2). 

We  performed  simulations  using  a  simple  model  of  the 
modulator  that  support  the  above  qualitative  reasoning  that 
field  screening  is  the  dominant  saturation  mechanism  in  these 
devices. 

IV.  Model 

The  model  we  used  to  simulate  charge  build-up  is  similar  to 
that  presented  by  Wood  er  a/.  [  1 1  ].  We  consider  the  reverse  bias 
p-i-n  structure  as  a  capacitor  with  the  quantum  wells  located 
between  the  plates.  The  differential  equations  governing  the 
areal  densities  of  electrons  n  and  holes  p  in  a  quantum  well 
illuminated  with  an  intensity  I  and  subjected  to  a  field  F  are 

~  =  a(F)I  -  Bnp  -  n/r,,  (4) 

at 

^  =  a{F)I-  Bnp-plTp.  (5) 

at 

Here  a(F)  is  the  field  dependent  absorption  of  one  well,  B  is 
the  recombination  constant  and  (rp)  is  the  electron  (hole) 
escape  time.  The  fact  that  the  escape  times  of  electrons  and 
holes  are  different  leads  to  accumulation  of  one  kind  of  carrier 
in  the  well  with  the  corresponding  loss  of  charge  neutrality. 
We  will  simulate  this  effect  by  the  introduction  of  a  layer  of 
charge  in  the  central  plane  of  the  quantum  well  with  a  charge 
density  given  by 

cr  =  e{p  —  n),  (6) 

For  high  incident  light  intensity,  the  net  positive  charge  can 
significantly  distort  the  electric  field  in  the  depletion  region. 
When  the  positive  charge  is  large  enough,  the  potential  at 
one  particular  well  may  be  reduced  until  it  approaches  the 
Fermi  level  of  the  n-type  material.  Under  these  conditions, 
electrons  can  diffuse  from  the  bulk  into  that  well,  effectively 
reducing  the  width  of  the  depletion  region.  In  each  step  of  the 
simulation,  we  determine  the  field  using  Gauss’s  law  and  the 
depletion  width  d  from  the  boundary  condition: 


where  Vh  and  are  respectively  the  built-in  and  external 
voltages.  The  background  doping  level  in  the  intrinsic  material 
of  wells  :and  barriers  is  assumed  negligible  for  simplicity  of 
the  model,  which  is  completed  with  charge  layers  at  both  ends 
of  the  depletion  region.  We  number  the  wells  starting  from  the 
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TABLE  1 

Parameters  Used  in  the  Simulation 


Temperature 

r  =  300  K 

Bandgap  difference"' 

AE  =  673  meV 

Fraction  of  discontinuity  in  valence  band*' 

Q  =  0.6 

Effective  mass  electrons  in  the  well"’  *^ 

m,  ir  =  0,0604 X/»o 

Effective  mass  electrons  in  the  barrier"'-^ 

ni,i,  =  0.1 19x/i^o 

Effective  mass  holes  in  the  welF'^^ 

=  0.471  X/;/o 

Effective  ma.ss  holes  in  the  barrier^' 

ffihh  =  0.62  xn/n 

Recombination  constant^' 

D  -  10"'^  cm-/s 

Diffusion  time  constant'* 

T,i  =  400  ns 

Well  width 

La-  =  10  nm 

Barrier  width 

Li,  =  15  nm 

Energy  of  /?  =  1  level  for  electrons 

E.  =  30.8  meV 

Energy  of  //  =  1  level  for  holes 

El,  =  26.3  meV 

Dielectric  constant 

f  =  12 

Area  of  the  device 

.4  =  360  X  510 //m- 

Number  of  wells 

=  24 

Built  in  voltage 

IL  =  I  V 

Focal  spot  radius 

r  =  10  pm 

Notes: 

a)  From  [18]. 

b)  From  calculations  following  the  procedure  in  [171. 

c)  Results  in  lifetime  around  I  ns. 

d)  Computed  from  a  diffusion  constant  of  10  cm”/s. 

e)  Order  of  magnitude  from  GaAs-AlGaAs  [19]  0  =  electron  mass. 

one  closest  to  the  p  region  of  the  diode  and  assume  the  light 
is  traveling  from  p  to  ii. 

The  numerical  values  of  the  parameters  used  in  the  simu¬ 
lation  are  listed  in  Table  I.  The  field  dependent  absorption  is 
obtained  from  low  intensity  transmission  measurements  under 
conditions  so  as  to  make  charge  build-up  effects  negligible. 
Besides  thermionic  emission  time,  computed  using  expression 
(1),  the  escape  time  for  both  electrons  and  holes  also  includes 
tunneling,  and  diffusion  of  the  carriers  in  the  plane  of  the  wells 
out  of  the  illuminated  area.  Although  it  is  negligible  for  our 
structure,  tunneling  was  included  in  order  to  be  able  to  evaluate 
different  device  designs  and  was  computed  using  formulas 
given  in  [3].  The  ambipolar  diffusion  time  for  a  diffusion 
coefficient  of  10  cm^/s  [19]  is  comparable  to  thermionic 
emission  when  the  field  is  low  but,  because  of  its  ambipolar 
nature,  it  does  not  contribute  to  the  charge  build-up. 

We  first  performed  steady  state  calculations,  in  which  we 
computed  the  mutually  dependent  field  and  charge  distribu¬ 
tions  in  a  self  consistent  manner.  Fig.  7  shows  the  steady 
state  carrier  density  (a)  and  the  field  distribution  (b)  for  the 
different  wells  in  the  structure.  These  results  were  obtained  for 
intensities  of  16  W/cm'  and  2.5  kW/cm“  while  the  external 
voltage  was  kept  constant  at  5  V.  In  Fig.  7(a)  the  circles 
(triangles)  represent  the  hole  (electron)  population  while  the 
intensity  is  denoted  by  the  use  of  empty  symbols  for  high 
intensity  and  solid  symbols  for  low  intensity.  In  Fig.  7(b)  the 
result  for  low  intensity  is  shown  by  the  solid  line  and  that  for 
the  high  intensity  by  the  dotted  line. 

In  the  depletion  region,  where  the  field  is  not  zero 
(Fig.  7(b)),  the  electrons  are  efficiently  swept  out  of  the  wells, 
and  there  is  accumulation  of  holes  (Fig.  7(a))  due  to  their 


(a) 


Fig.  7.  (a)  Computed  steady  state  distribution  of  electron  (triangles)  and 

holes  (circles)  area!  densities  with  an  applied  voltage  of  5  V  for  intensities  of 
16  W/cm"  (solid  symbols)  and  2.5  kW/cm“  (empty  symbols).  The  p  region 
is  located  to  the  left  of  the  wells.  The  light  travels  from  left  to  right,  (b) 
Computed  steady  state  electric  field  for  the  same  conditions  as  in  (a).  Low 
intensity  data  is  plotted  with  solid  line,  and  high  intensity  data  with  dashed 
line. 

longer  escape  time.  For  5  V  of  external  bias,  and  in  absence 
of  illumination,  all  the  wells  would  be  subjected  to  an  electric 
field  around  10^  V/cm.  Due  to  the  presence  of  net  charge  in 
the  wells,  the  field  is  not  the  same  for  all  the  wells  and  those 
where  the  field  is  lower  have  higher  absorption  and  longer 
escape  times,  leading  to  higher  densities  to  the  right  end  (n 
side)  of  the  depletion  region  in  Fig  7.  The  accumulation  of 
carriers  is  not  very  important  when  the  intensity  is  low,  and 
most  of  the  wells  experience  a  value  of  the  field  close  to  that 
in  the  absence  of  light.  Notable  changes  occur  in  the  charge 
distribution  when  the  intensity  is  increased.  As  the  carrier 
density  increases,  the  charge  accumulates  and  the  width  of  the 
depletion  region  varies  accordingly.  Once  a  well  is  out  of  the 
depletion  region,  the  field  goes  to  zero  with  a  corresponding 
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change  in  transmission  and  with  equal  populations  of  electrons 
and  holes.  As  shown  in  Fig.  7(b),  at  high  intensity  only  few 
of  the  wells  experience  the  field,  resulting  in  a  reduction 
of  the  modulation  with  respect  to  the  low  intensity  case. 
These  qualitative  results  of  the  model  show  the  same  behavior 
that  was  observed  in  the  experiment  (Fig.  4).  In  spite  of 
uncertainties  in  certain  material  parameters  and  its  simplicity, 
the  model  also  predicts  the  modulation  with  less  than  20% 
error  respect  to  the  experimental  values  for  intensities  up  to  the 
measured  /sat  (~4  kW/cm“).  At  higher  intensities,  however, 
the  discrepancies  are  larger  and  it  does  not  account  for  all 
the  reduction  of  modulation  that  is  observed.  This  difference 
might  be  due  to  effects  not  included  in  our  model  that  are 
nevertheless  suggested  by  the  shape  of  the  spectra  (Figs.  1 
and  4),  such  as  slight  saturation  of  the  excitonic  absorption 
or  unintentional  doping  of  the  material  in  the  intrinsic  region. 
Another  possible  source  of  error  is  inaccuracy  associated  with 
the  use  of  (2)  to  estimate  the  thermionic  escape  times. 

The  simulation  of  the  time  response  of  the  modulators  starts 
from  the  steady  state  situation  for  14  =  0  V,  and  follows  the 
changes  produced  by  an  increase  of  applied  voltage  with  a 
risetime  of  4  ns.  In  each  step  of  the  simulation,  the  following 
tasks  are  performed. 

1)  An  estimate  of  the  external  charge  flow  needed  to 
establish  the  new  voltage  across  the  device  is  computed. 
As  the  flow  of  carriers  external  to  the  MQW  structure  is 
much  faster  than  the  time  scale  of  the  escape  from  the 
wells,  we  assume  that  the  external  carriers  flow  so  as  to 
accommodate  the  internal  changes  in  accordance  to  the 
boundary  condition  (7). 

2)  The  rate  of  change  of  the  carrier  population  in  all  the 
wells  is  computed  by  means  of  expression  (4)  and  (5) 
and  is  used  to  calculate  the  new  carrier  densities  after  a 
short  time  interval  dt  The  total  transmission  coefficient 
of  the  structure  is  also  computed. 

3)  The  external  charge  is  recomputed  to  meet  the  boundary 
condition  (7)  taking  into  account  the  new  internal  charge 
distribution.  This  is  usually  a  small  correction  of  the 
values  calculated  in  step  1.  Also  in  this  step  the  width 
of  the  depletion  region  is  modified  when  the  new  internal 
charge  distribution  makes  it  necessary. 

Typical  results  of  the  simulation  are  plotted  in  Fig.  8,  which 
shows  the  simulated  temporal  evolution  of  the  transmission  of 
the  device  for  different  conditions  of  incident  intensity  and 
bias  voltage.  All  curves  have  been  normalized  to  the  steady 
state  modulation  for  ease  of  comparison.  Curves  a  and  b 
correspond  to  50  W/cm“  of  incident  intensity,  while  curves  c 
and  d  were  computed  using  2.5  kW/cm^.  The  applied  voltage 
changes  from  0  to  10  V  (a,c)  and  0  to  4  V  (b,d)  with  a 
risetime  of  4  ns.  As  can  be  seen  by  comparing  these  results 
with  the  corresponding  experimental  curves  shown  in  Fig.  6, 
the  model  reproduces  qualitatively  the  temporal  changes  in 
transmissivity  observed  in  the  experiment.  At  low  intensity 
(a  and  b),  an  increase  in  bias  voltage  is  accompanied  by  a 
decrease  in  switching  time,  in  agreement  with  the  data  in 
Figs.  5  and  6,  while  at  high  intensity  (c  and  d)  that  change  is 
much  less  noticeable,  also  in  agreement  with  the  experiment. 


Fig,  8.  Simulated  temporal  evolution  of  the  transmission  of  the  device 
normalized  to  the  steady  stale  modulation.  The  conditions  of  the  simulation 
are  the  same  as  those  of  the  experiment  in  Fig,  6. 


The  switching  time  for  10  V  bias  at  high  intensity  is  slower 
than  that  of  10  V  at  low  intensity.  The  origin  of  this  effect 
is  the  change  in  depletion  width  with  voltage.  As  can  be  seen 
from  the  steady  state  data  in  Fig.  7(a),  most  of  the  carriers  are 
accumulated  near  the  end  of  the  depletion  region.  A  sudden 
increase  of  voltage  moves  the  boundary  of  the  depletion  region 
closer  to  the  n  region,  leaving  several  wells  in  a  situation  out 
of  equilibrium  with  low  field  and  consequently  ((1))  long  time 
to  reach  the  new  steady  state  situation. 

While  the  simulation  gives  numerical  values  for  the  risetime 
which  do  not  match  exactly  those  obtained  experimentally,  this 
difference  is  again  to  be  expected  given  the  simplicity  of  the 
model  and  the  fact  that  some  of  the  parameters  used  are  only 
estimates.  Nevertheless,  the  results  of  the  model  give  solid 
support  to  the  qualitative  reasoning  presented  in  the  previous 
section  leading  to  the  conclusion  that  field  screening  due  to 
carrier  build-up  is  the  most  important  effect  determining  the 
observed  saturation  of  the  modulation  in  the  InGaAs-InGaP 
modulators  investigated.  We  also  ran  our  simulation  for  a 
different  structure,  with  barriers  of  7.5  nm  thickness.  This 
change  makes  tunneling  the  dominant  escape  mechanism, 
and  produces  approximately  a  tenfold  increase  in  saturation 
intensity  and  switching  speed  limited  by  the  risetime  of  the 
electrical  pulse.  Although  the  simplicity  of  the  model  does  not 
allow  a  detailed  prediction,  it  shows  the  possibility  of  increase 
both  the  saturation  intensity  and  speed  of  the  device  by  means 
of  changes  in  the  structure. 

V.  Conclusion 

We  performed  a  detailed  study  of  the  performance  of 
InGaAs-InGaP  multiple  quantum  well  modulators  under 
different  conditions  of  incident  intensity  and  operating 
voltage.'  We  found  that  the  saturation  intensity  is  relatively 
constant  with  a  value  of  (1.3i0.1)  kW/cm^  up  to  V 
bias  where  the  saturation  begins  to  increase  linearly  with 
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voltage  to  reach  (3.7±0.1)  kW/cm^  at  10  V  bias.  Spectra  at 
different  powers  and  applied  voltages  show  that  the  decrease 
of  modulation  with  power  is  not  due  to  saturation  of  the 
excitonic  absorption.  We  performed  measurements  on  the 
time  response  of  the  devices,  which  not  only  characterizes 
the  modulation  speed  of  the  device,  but  also  confirm  that 
the  most  important  cause  of  the  saturation  is  field  screening 
due  to  charge  build-up.  The  obtained  switching  times  ranging 
between  10  and  90  ns  are  not  limited  by  the  capacitance  of  the 
device,  but  by  the  escape  time  of  the  photogenerated  carriers 
from  the  wells.  In  the  devices  studied  the  escape  time  was 
dominated  by  thermionic  emission,  which  is  slower  than  in 
other  material  systems  due  to  the  larger  difference  in  band  gap 
between  the  well  and  barrier  materials.  In  order  to  improve  the 
speed  of  InGaAs-InGaP  devices  while  keeping  the  advantage 
of  a  transparent  substrate,  the  structure  must  be  carefully 
engineered  to  accelerate  the  escape  of  the  carriers  from  the 
wells,  for  example  by  using  thinner  barriers  that  would  allow 
tunneling.  We  performed  simulations  using  a  simple  model 
of  the  modulator  which  includes  field  screening  effects  that 
confirm  field  screening  to  be  the  dominant  physical  mechanism 
governing  the  saturation  and  time  response  of  these  devices. 
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Photoluminescence  and  Interface  Abruptness  in 
InGaAsP/InGaAsP  Quantum  Wells 

L.M.  WOODS,  P.  SILVESTRE,  P.  THIAGARAJAN,  G.A.  PATRIZI, 
and  G.Y.  ROBINSON 

Department  of  Electrical  Engineering,  Colorado  State  University, 

Fort  Collins,  CO  80523 

KM.  JONES  and  M.  AL-JASSIM 

National  Renewable  Energy  Laboratory,  Golden,  CO  80401 

Quantum  well  (QW)  structures  consisting  of  InGaAsP  wells  and  InG-aAsP 
barriers  grown  by  gas-source  molecular  beam  epitaxy  have  been  examined  by 
low  temperature  photoluminescence  (PL)  in  order  to  evaluate  the  contributions 
of  compositional  fluctuations  in  the  quaternary  alloy  and  of  interface  roughness 
to  the  PL  linewidth.  The  well  material  was  InGaAsP  with  a  bandgap  correspond¬ 
ing  to  a  wavelength  of  1.3  pm  and  the  barrier  material  was  InGaAsP  of  1.15  pm. 
The  theory  for  QW  excitonic  linewidths  as  a  function  of  well  thickness  L^  due  to 
fluctuations  in  alloy  composition  has  been  extended  to  include  the  case  of 
quaternary  InGaAsP  barrier.  If  the  interfaces  are  atomically  abrupt,  the  linewidth 
is  dominated  by  compositional  fluctuations  in  the  well  at  large  L^  and  composi¬ 
tional  fluctuations  in  the  harrier  at  small  L^.  The  theory  predicts  a  weak 
dependence  of  the  linewidth  on  L^  since  the  composition  of  the  well  and  barrier 
are  similar.  For  rough  heterointerfaces,  the  theory  indicates  the  usual  increase 
in  linewidth  with  decreasing  L,.  Photoluminescence  measurements  at  13K  in 
arrays  of  single  InGaAsP/InGaAsP  QWs  with  L^  from  1.0  to  6.0  nm  show  only  a 
weak  variation  of  the  full  width  at  half  maximum  (FWHM)  with  L^,  m  agreement 
with  the  theory  for  smooth  interfaces.  Furthermore,  the  lowest  measured 
FWHM  of  8.9  meV  was  foimd  for  a  narrow  well  of  L^  =  1.8  nm,  indicating  the 
InGaAsP/InGaAsP  interfaces  are  smooth  and  that  the  PL  linewidth  is  domi¬ 
nated  by  compositional  fluctuations. 

Key  words:  Gas  source  molecular  beam  epitaxy  (GSMBE),  InGaAsP, 
photoluminescence 


INTRODUCTION 

The  III-V  materials  system  Inj_^Ga^yPi_y  grown 
epitaxially  on  InP  has  been  used  extensively  in  quan¬ 
tum  well  (QW)  optoelectronic  and  transport  devices. 
The  wide  range  of  bandgap  energies  available  in 
InGaAsP  alloys,  by  varying  the  alloy  composition  x 
and  y,  allows  the  device  desi^er  great  flexibility  in 
optimizing  the  device  for  maximum  performance.  F or 
example,  low  threshold  QW  lasers  operating  at  the 
wavelength  of  1.3  pm  utilize  InGaAsP  wells  of  one 
composition  and  InGaAsP  barrier  and  separate 
confinement  layers  of  different  compositions.'  In  opti¬ 
cal  modulators  using  the  quantum  confined  Stark 
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effect,  the  contrast  ratio  can  be  maximized  and  the 
applied  bias  minimized  by  using  InGaAsP  wells  and 
InGaAsP  barriers  of  the  appropriate  compositions.^ 
The  abruptness  of  QW  heterointerfaces  can  greatly 
influence  the  properties  of  heterostructoe  devices 
derived  from  quantum-size  effects.  Previous  studies 
have  examined  the  interface  abruptness  of  InGaAsP/ 
InP  QWs,3-5  and  InGaAs/InP  QWs,^^  but  there  ap¬ 
pears  to  be  little  information  about  the  abruptness  of 
InGaAsP/InGaAsP  heterointerfaces. 

We  report  here  the  results  of  a  study  of  a  InGaAsP/ 
InGaAsP  QW  interface  using  the  photoluminescence 
:  characterization  technique.  The  InGaAsP  structures 
were  grown  by  gas  source  molecular  beam  epitaxy 
(GSMBE)  on  InP  and  utilized  two  lattice-matched 
quaternaries;  Ino7Gag3AsjgPo4  as  the  well  material 
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with  a  room  temperature  bandgap  emission  at  a 
wavelength  of  approximately  1.3  pm  (here  referred  to 
as  1.3Q),  and  ggGa^  53  as  the  barrier  mate¬ 

rial  with  emission  at  1.15pm(1.15Q).  Quantum  wells 
of  these  compositions  are  vddely  used  in  optical  de¬ 
vices  operating  at  1.3  pm,  but  determination  of  the 
atomic  abruptness  of  the  interface  is  difEcult  by  the 
usual  methods  of  transmission  electron  microscopy 
(TEM)  and  x-ray  diffraction;  the  contrast  between  the 
well  and  barrier  materials  is  small  since  the  composi¬ 
tions  are  so  similar.  However,  the  low  temperature 
photoluminescence  (PL)  line  width  for  excitons  in  quan¬ 
tum  wells  with  well  thickness  less  than  about  3.0  nm 
depends  strongly  on  interface  abruptness.  If  other 
contributions  to  the  linewidth  are  properly  accounted 
for,  the  PL  linewidth  can  be  used  to  estimate  the 
heterointerface  abruptness.  Here,  the  theory  of  PL 
linewidth  for  interface  roughness  was  extended  to 
include  the  contributions  from  compositional  fluctua¬ 
tions  in  both  the  well  and  barrier  InGaAsP  layers.  A 
comparison  of  the  experimental  and  theoretical  PL 
linewidths  as  a  function  of  well  width  was  then  used 
to  ascertain  the  interface  roughness  in  InGaAsP(  1.3Q)/ 
InGaAsPd.lSQ)  QWs. 

THEORY 

Low  temperature  PL,  imder  low  power  excitation, 
of  QW  samples  with  sharp  interfaces  and  uniform 
composition  yields  narrow  luminescence  peaks.  Quan¬ 
tum  wells  can  be  dominated  by  excitons  which  have 
an  intrinsic  emission  spectrum  of  a  lifetime  broad¬ 
ened  delta  function.  Thus,  the  linewidth  of  a  purely 
excitonic  transition  is  very  narrow  making  its  peak 
ideal  for  monitoring  small  inhomogeneous  contribu¬ 
tions  to  the  linewidth.  The  inhomogeneous  contribu¬ 
tions  can  result  from  compositional  variations,  and 
disordered  and/or  graded  heterointerfaces  in  the  QW 
structure.  These  statistical  mechanisms  act  to  ran¬ 
domize  the  discrete  level  of  quantum  confinement, 
thus  spreading  out  the  energy  of  luminescence  and 
adding  to  the  width  of  the  luminescence  peak. 

What  follows  is  an  analysis  that  reviews  the  contri¬ 
butions  to  the  photoluminescence  linewidth  of  QW 
excitons,  so  that  the  contribution  from  interface  rough¬ 
ness  can  be  identified  and  characterized.  If  we  let  T  be 
the  full  width  at  half  maximum  (FWHM)  of  the  PL 
exciton  peak,  then 

n  =  r/  -H  d) 

where  is  the  homogeneous  part  and  is  the 
inhomogeneous  part  of  F.  Linewidths  can  be  added  in 
this  fashion  if  each  component  is  approximated  as  a 
Gaussian  distribution  and  is  statistically  independent, 
r^is  temperature  dependent  and  can  be  considered  to 
result  from  two  contributions: 

=  V  +  (2) 

where  the  contribution  from  longitudinal  opti¬ 
cal  phonon  assisted  broadening  and  is  the  contri¬ 
bution  from  carrier  lifetime  broadening.  The  room 
temperature  F^^^  value  has  been  determined  from 


measurements  of  exciton  absorption  peaks  at  300K  by 
Sugawara  et  al.  as  8.9  meV  for  1.33  jim  InGaAsP.® 
Here  the  value  of  Fpg.j.at  lOK  was  estimated  to  be  0.5 
meV  by  using 

r'phT(T) = Fpg  /[exp(ft(Ogo  /  kT)  - 1]  (3) 

where  ^'pH  is  a  proportionality  constant  determined 
from  the  exciton  absorption  measurements  and 
is  the  longitudinal  optical  phonon  energy  of  32  meV 
for  InGaAsP(1.3Q).®’® 

Carrier  lifetime  broadening  is  due  to  the  time  it 
takes  to  re-establish  the  carrier  distribution  function 
within  the  energy  bands  after  a  recombination  event. 
This  process  does  not  take  place  instantaneously,  as 
represented  by  a  nonzero  lifetime  of  the  final  state. 
The  lifetime  has  been  found  to  be  as  small  as  500  ps 
for  InGaAs/InP  and  InGaAsP/InP  QWs  less  than  20.0 
nm  in  thickness.^  Using  Heisenberg’s  uncertainty 
principle  and  a  lifetime  of 500  ps,  the  maximum  value 
for  Flb  is  estimated  to  be  0.01  meV  for  InGaAsP/ 
InGaAsP  QWs. 

Then  using  Eq.  (2),  we  estimate  F^ = 0.50  meV  at  low 
temperatures.  But  the  measured  PL  FWHM  values 
for  QWs  of  InGaAs  and  InGaAsP  are  typically  5  meV 
or  larger.®"^  Hence,  the  homogeneous  contributions  to 
the  PL  linewidth  can  be  neglected. 

The  inhomogeneous  linewidth,  F^  in  Eq.  (1),  is 
statistical  in  nature  and  is  due  to  random  variation  in 
the  QW  structure  from  point  to  point.  It  can  also  be 
taken  to  consist  of  two  independent  contributions 

Fo2  =  Fj®  +  F/  (4) 

where  F^  is  the  contribution  due  to  the  bandgap 
variation  from  localized  fluctuations  in  alloy  composi¬ 
tion,  and  Fj  is  the  contribution  from  quantum  size 
effects  due  to  interface  roughness. 

Fj  is  the  linewidth  at  absolute  zero  temperature 
that  results  from  random  spacial  fluctuations  in  the 
InGaAsP  composition.  These  spacial  composition  fluc¬ 
tuations  result  in  a  broadening  of  the  PL  emission 
energy  due  to  the  corresponding  fluctuation  of  the 
bandgap  energy.  The  fluctuations  occur  on  a  micro¬ 
scopic  scale  because  of  the  random  arrangement  of 
the  In  and  Ga  atoms  on  the  group-IIII  lattice  sites  and 
the  As  and  P  atoms  on  group-V  sites.  Thus,  even  an 
epitaxial  layer  of  uniform  average  composition  will 
exhibit  compositional  broadening  in  PL.  A  theory  for 
compositional  fluctuations  has  been  developed  by 
Schubert  et  al.^®  for  bulk  ternary  alloys,  and  we  have 
extended  the  theory  to  QWs  with  quaternary  wells 
and  quaternary  barriers  by  using  the  theory  of  ran¬ 
dom  variables  vdth  joint  probability  distributions  and 
a  covariance  of  zero.^®  The  compositional  fluctuations 
are  evaluated  over  the  excitonic  volume,  which  is 
approximated  as  an  ellipsoid  and  given  by: 

V3,  =  47rZ,,pJ/3  (5) 

where  and  are  the  lateral  extension  and  perpen¬ 
dicular  extension,  respectively,  of  the  exciton 
wavefunction  relative  to  the  plane  of  the  QW.  This 
volume  is  quantum  mechanical  in  nature  and  thus. 
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there  is  a  nonzero  probability  of  finding  the  exciton  in 
the  barrier  material.  Schubert  et  al.^*  showed  that  the 
linewidth  is  proportional  to  the  product  of  the  rate  of 
change  of  bandgap  with  composition  and  the  stan¬ 
dard  deviation  of  the  compositional  fluctuations.  Here 
we  take  the  standard  deviations  of  composition  fluc¬ 
tuations  in  to  be  weighted  according  to  the  prob¬ 
ability  of  the  exciton  being  in  the  barrier,  Pj,,  or  the 
well,  P  .  =  1  -  Pb-  Thxis,  the  exciton  emission  linewidth 
due  to  alloy  fluctuations  in  Inj_,jGajjASyP QW  struc¬ 
tures  can  ^  written  as 

r,=  (6) 


VdE,YxJi-xJP,  rdE.y  yji-yjp. 

idxj  UyJ  Kv,,, 

fdE^Y  xb(^-^b)Pb ,  fdE^bYyb(i-yb)Pb 
dx  J  I  dy  J  KV,,,  _ 

In  Eq.  (6),  E^^  is  the  bandgap  energy  of  the  barrier,  K 
=  4a  is  the  cation  or  anion  density  in  the  zinc-blende 
lattice  with  lattice  constant  a„,  and  x^,  y,,,  x^,  y^  are  the 
Ga  and  As  concentrations  in  the  well  and  barrier, 
respectively.  The  PL  peak  is  located  at  the  energy 

E  =  E  +  E  +  E  -  Er  (7) 

"^gw  n,e  m,hh  B 


Here  E  is  the  bandgap  energy  of  the  well,  E„^  and 
E  are^the  quantum  confined  levels  of  the  electrons 
an^'’heavy  holes,  respectively,  and  Eg  is  the  quasi- 
two-dimensional  exciton  binding  energy.  The  depen¬ 
dence  of  E^,  E^„  and  E^^  on  x  and  y  were  found  from 
the  InGaAsP  bandgap  data  of  Moon  et  al.^"*  along  with 
the  solution  to  Schrodinger’s  equation  for  a  finite 
potential  well  of  depth  AE^  for  the  electron  and  AE^  for 
the  heavy  hole,  where  is  the  conduction  band 
offset  energy  and  AE^  is  the  valence  band  offset  energy 
for  the  heterojunction.  The  value  of  dE/dx  (or  dE/dy) 
was  then  evaluated  by  differentiating  each  of  the 
terms  in  Eq.  (7)  with  respect  to  x  (or  y),  neglecting  the 
exciton  binding  energy  dependence  with  composition. 

Equation  (6)  has  been  used  to  calculate  the  de¬ 
pendence  of  r.  on  well  width  L,  for  ground-state 
excitons  inInGaAsP(1.3Q)/InGaAsP(1.15Q)  QWs,  and 
the  results  are  shown  by  the  dashed  curve  in  Fig.  1. 
For  this  combination  of  alloy  compositions,  we  used 
AE„  =  58  meV  and  AE^  =  88  meV  assuming  AE^:  AEv  = 
40:60,^®  and  the  values  of  effective  maisses  from  Restorff 
et  al.^®  The  parameters  Pg,  P,,,  p„,  and  have  been 
calculated  as  a  function  of  well  width  for  InGaAs/InP 
QWs  by  Herman  et  al.^’  and  approximated  here  to  be 
the  same  for  the  quaternary  system  in  Fig.  1 .  varies 
from  the  bulk  value  of  the  well  material  (3.4  meV)  for 
wide  wells,  to  the  bulk  value  of  the  barrier  material 
(4.2  meV)  as  the  exciton  wavefunction  penetrates  into 
the  barrier  for  veiy  narrow  wells.  The  peak  in  P ^  that 
occurs  at  about  =  2  nm  is  when  the  contributions  to 
the  linewidth  from  both  the  well  and  barrier  are 


comparable.  ,  •  i 

Next  consider  P^,  the  contribution  to  the  inho¬ 
mogeneous  linewidth  due  to  interface  roughness  or 


variation  in  well  thickness.  Individual  monolayer 
steps  at  the  interface  will  act  to  randomize  the  thick¬ 
ness  of  the  well,  L^,  thus  randomizing  the  quantum 
confined  levels  and  increasing  the  linewidth.  Theo¬ 
retically,^® 


dL,  dL,  _ 


(8) 


where  AL^  is  the  width  of  the  interface.  The  depen¬ 
dence  of  the  exciton  binding  energy  Eg  with  has 
been  ignored  in  Eq.  (8).  The  solid  curves  in  Fig.  1  show 
the  variation  in  Pj  as  a  function  of  well  width  for  the 

InGaAsP(1.3Q)/InGaAsP(1.15Q)QW system,  with  AL, 

varying  from  0.5  monodiatomic  layers  (ml)  to  3.0  ml. 
For  greater  than  1  nm,  the  linewidth  decreases 
with  increasing  well  width  in  the  usual  manner,  and 
for  L,  less  than  about  1  nm,  the  linewidth  begins  to 
decrease  for  decreasing  as  a  result  of  the  finite 
depth  of  the  electron  and  hole  wells. 

Based  on  the  above  theory,  we  can  expect  the 
measured  linewidth  P  to  be  due  to  P j,  P2,  or  both.  For 
narrow  wells  (i.e.,  less  than  about  1.5  nm),  P  =  P2 
and  the  linewidth  is  dominated  by  interface  rough¬ 
ness  if  AL  >  1  ml.  For  wide  wells  (L^  greater  than 
about  5  nin),  P  =  Pj  and  the  linewidth  is  primarily  due 
to  composition  fluctuations  and  is  relatively  insensi¬ 
tive  to  interface  roughness. 

EXPERIMENTAL  PROCEDURE 

The  samples  were  grown  by  GSMBE  on  (100)  SI  InP 
substrates  at  approximately  500°C  with  elemental 
Ga  and  In  and  thermally  cracked  AsHg  and  PHg  as 
source  materials.^®  Growth  was  paused  at  each 
heterointerface  for  20  s,  to  allow  the  surface  to  anneal 
in  the  group-V  molecular  beams  and  to  allow  stabili¬ 
zation  of  the  gas  flows  before  starting  growth  of  the 
next  layer.  Double  crystal  x-ray  diffraction  and 
photoluminescence  of  companion  samples  were  used 
to  check  lattice  mismatch  (<5  x  10"*)  and  alloy  compo- 


Fig.  1 .  The  theoretical  PL  linewidth  due  to  inhomogeneous  broadening 
from  compositional  fluctuations,  and  interface  ^ 

function  of  well  width,  L,  for  ground  state  excitons  in  lnGaAsP(i  .dU)/ 
lnGaAsP(1.15  Q)  QWs.  AL^  is  the  heterointerface  roughness. 
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Fig.  2.  Transmission  electron  (200)  dark  field  micrograph  of  the  cross 
section  of  an  array  of  InGaAsP/InGaAsP  single  QWs  lattice  matched 
to  InP.  The  wells  are  lnGaAsP(1.28Q)  and  the  barriers  are 
inGaAsP{1.15Q). 


Energy  (eV) 

Fig.  3.  Photoluminescence  spectra  at  13K  of  single  QWs  In  two 
samples  with  lnGaAsP(1 .1 5Q)  barriers:one  with  three  lnGaAsP(1 .28Q) 
wells  and  the  other  with  two  inGaAsP{1 .33Q)  wells.  The  exciton  peaks 
for  each  QW  are  labeled  with  the  corresponding  well  width. 

array  of  single  QWs  with  InGaAsP  wells  of  nominal 
1.3Q  composition  separated  by  1.15Q  barriers  of  12.5 
nm  thickness,  a  thick  ‘Tjulk”  layer  of  the  1.3Q  alloy, 
and  an  InP  cap  layer.  A  TEM  (200)  dark  field  micro¬ 
graph  of  a  cross  section  the  first  sample  is  shown  in 
Fig.  2.  The  sample  contains  three  QWs  with  InGaAsP 
composition  of  1.28Q  and  thicknesses  of  1.2  ±  0.6, 2.4 
±  0.6,  and  5.7  ±  0.3  nm.  The  “bulk”  or  reference  layer 
is  86.3  nm  thick.  The  low  contrast  resulting  from  the 
similar  compositions  of  1.28Q  and  1.15Q  made  it 
difficult  to  resolve  individual  atomic  layers  at  the 
heterointerfaces  in  high  resolution  lattice  imaging 
with  TEM.  Thus,  measurements  of  the  interface 
abruptness  in  this  mode  could  not  be  obtained.  How¬ 
ever,  all  the  QWs  are  clearly  resolvable  in  Fig.  2 
including  the  1.2  nm  well.  Rough  measurements  of 


Fig.  4.  The  shift  in  QW  peak  energies  with  respect  to  the  bulk  layer  as 
a  function  of  well  width,  L^,  from  PL  measurements  at  1 3K.  Also  shown 
is  the  theory  for  (solid  line)  lnGaAsP(1.28Q)/lnGaAsP(1.15Q)  and 
(dashed  line)  lnGaAsP(1.33Q)/inGaAsP(1.15Q)  QWs. 

the  lateral  uniformity  from  dark  field  TEM  images 
provided  an  estimate  of  a  lateral  uniformity  of  the 
interfaces  of  about  2  ml  (0.6  nm)  or  less  over  a  lateral 
distance  of  about  50.0  nm.  The  sample  in  Fig.  2  will  be 
referred  to  as  the  1.28Q/1.15Q  sample. 

The  second  sample  consisted  of  two  QWs  of  1.8  ±  0.6 
and  5.2  ±  0.6  nm  thickness  and  a  reference  layer  of 
99.9  nm  thickness.  The  InGaAsP  composition  in  the 
wells  was  fotmd  to  be  1.33Q  and  the  layer  thicknesses 
were  estimated  from  the  growth  conditions.  This 
sample  will  be  referred  to  as  the  1.33Q/1.15Q  sample. 

Photoluminescence  measurements  were  made  with 
excitation  provided  by  the  488  nm  line  of  an  Ar  ion 
laser  with  incident  power  of  about  10  W/cm^  at  the 
sample.  At  this  power  level,  broadening  of  the  PL 
peak  due  to  band  filling  was  found  to  be  insignificant. 
The  InGaAsP  QW  samples  were  mounted  inside  a 
closed  cycle  refrigerator  for  measurement  in  the  tem¬ 
perature  range  of  10-300K.  The  luminescence  wave¬ 
length  selection  was  accomplished  with  a  0.5  m 
monochrometer,  and  the  intensity  detected  with  LNj- 
cooled  Ge  diode  and  processed  using  a  lock-in  amph- 
fier  and  computer.  Using  1  mm  slit  widths,  the  total 
instrument  line  broadening  (3.4  meV)  resulted  pri¬ 
marily  from  the  monochrometer  and  was  deconvolved 
from  the  measured  PL  linewidths. 

RESULTS 

The  PL  spectra  at  13K  for  the  two  samples  are 
shown  in  Fig.  3.  The  exciton  peaks  for  each  QW  are 
labeled  with  the  corresponding  well  width  and  the 
free  carrier  recombination  peak  for  the  bulk  layer  is 
indicated.  As  the  temperature  of  the  measurement 
was  increased,  the  peaks  shifted  to  lower  energies  in 
a  manner  consistent  with  the  dependence  of  the 
bandgap  of  the  well  material  on  temperature.^®  At 
300K,  only  the  bulk  peak  was  resolvable  in  both 
samples,  with  a  FWHM  of  60.5  meV  for  the  1.28Q/ 
1.15Q  sample  and  42.7  meV  for  the  1.33Q/1.15Q 
sample. 

The  shift  in  QW  peak  energies  with  respect  to  the 
bulk  layer  peak  as  a  function  of  width  L^  are  plotted 
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along  with  the  theory  in  Fig.  4.  Although  the  1.33Q/ 
1.15Q  sample  should  exhibit  a  larger  shift  than  the 
1.28Q/1.15Q  sample,  the  observed  shift  is  about  the 
same  for  both  samples.  However,  the  general  agree¬ 
ment  of  the  calculated  curves  with  the  data  confirms 
the  presence  of  quantum  confinement  and  the  large 
energy  shift  for  the  1.2  nm  QW  is  indicative  of  a 
square-like  well,  as  the  quantum  shift  changes  very 
rapidly  with  well  width  for  narrow  weUs.  Note  the  1.2 
nm  well  produces  an  energy  shift  of  over  90%  of  the 
maximum  possible. 

The  FWHM  linewidths  of  the  QW  peaks  measured 
at  13K  have  been  plotted  as  a  function  of  well  thick¬ 
ness  in  Fig.  5.  Also  shown  in  Fig.  5  is  the  theoretical 
inhomogeneous  linewidth  usingthe  contributions  from 
composition  fluctuation  Fj  and  interface  roughness 
calcidated  in  Fig.  1.  Although  there  is  appreciable 
scatter  in  the  data,  it  is  apparent  that  there  is  little 
variation  in  the  FWHM  with  over  the  range  mea- 
sirred,  in  agreement  with  the  theoretical  trend  for 
smooth  heterointerfaces.  In  particular,  for  small 
where  the  FWHM  is  most  sensitive  to  the  interface 
roughness,  the  measured  hnewidths  are  less  than  10 
meV,  indicating  less  than  1  ml  roughness  at  the 
InGaAsP/InGaAsP  interfaces  for  the  smallest  QWs. 
The  smallest  FWHM  was  8.9  meV  at  =  1.8  nm. 

CONCLUSIONS 

The  interface  abruptness  between  InGaAsP(1.3Q) 
and  InGaAsPd.lSQ)  has  been  studied  using  PL  of 
QW  structures  grown  by  gas  soiuce  MBE.  The  strong 
quantum  shift  in  energy  measured  for  the  narrow 
wells  was  indicative  of  square-hke  wells  with  abrupt 
interfaces.  This  was  further  confirmed  by  TEM  cross- 
sectional  images  that  showed  good  uniformity  and 
values  of  the  well  thicknesses  that  correlated  with  the 
growth  conditions.  The  PL  linewidth  was  studied  as  a 
function  well  size  for  contributions  from  alloy  compo¬ 
sitional  variations  and  interface  roughness.  It  was 
determined  that  the  linewidth  was  dominated  by 
alloy  compositional  variations  and  the  contribution 
from  interface  roughness  was  small.  For  em  interface 
roughness  of  2  ml  or  larger,  we  would  expect  a  large 
increase  in  the  hnewidth  for  the  narrower  wells.  This 
was  not  observed  and  the  linewidth  for  the  narrow 
wells  in  fact  decreased,  consistent  with  alloy  compo¬ 
sition  variations  dominating  the  PL  linewidth  and 
smooth  interfaces  of  1  to  2  monolayer  in  roughness. 
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interface  roughness  of  Fig,  1 . 
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Spectral  transmission,  reflection,  and  photocurrent  absorption  data  obtained  on  gas-source 
molecular  beam  epitaxy  grown  InGaAs/InP  multiple  quantum  well  (MQW)  and  superiattice 
p-i-n  diode  structures  demonstrate,  for  the  first  time  in  this  materials  system,  that  similar 
modulation  to  MQW  structures  can  be  achieved  using  superlattices,  but  at  significantly  lower 
operating  voltages.  Specifically,  we  have  observed  photocurrent  absorption  changes  of  as  much 
as  58%,  transmission  changes  of  8.2%,  and  reflection  changes  of  32%  for  applied  biases  of  only 
4  V,  in  nonresonant  modulators  operating  at  a  wavelength  --1.5  fim.  These  results  encourage 
the  possibility  of  employing  such  devices  in  fast,  high  density  optical  modulator  arrays  operating 
over  the  1. 3-1.6  /tm  range. 


High  speed  III-V  semiconductor  light  modulators  are 
of  considerable  interest  for  a  variety  of  applications  includ¬ 
ing  low  chirp  external  laser  modulators,  photonic  switch¬ 
ing  fabrics,  optical  logic  circuits,  and  smart  pixel  arrays.* 
Most  work  on  such  devices  has  concentrated  on  the  use  of 
the  quantum  confined  Stark  effect  (QCSE)^  to  modify  the 
spectral  characteristics  of  multiple  quantum  well  (MQW) 
stacks  of  various  material  combinations  resulting  in  voltage 
controlled  modulation  of  either  transmitted  or  reflected 
optical  signals.  Absolute  changes  in  reflectivity  as  high  as 
77%  have  been  achieved  using  asymmetric  Fabry  Perot 
resonator  structures  in  strained  InGaAs/AlGaAs^  for  op¬ 
erating  voltages  in  the  range  of  20  V. 

Such  large  voltages,  however,  result  in  very  large 
power  dissipation  densities  when  the  devices  are  switched 
repetitively  at  high  rates.  The  power  P  is  required  to  charge 
and  discharge  the  capacitance  C  associated  with  the  p-i-n 
diodes  and  is  given  by  CFV2t,  where  r  is  the  switching 
time  and  V  is  the  switching  voltage.'^  Reductions  in  the 
operating  voltage  of  such  devices  below  5  V  would  be 
desirable,  both  to  reduce  their  dynamic  power  consump¬ 
tion  and  hence  allow  their  application  in  high  speed  dense 
arrays,  as  well  as  provide  compatibility  with  Si  comple¬ 
mentary  metal-oxide  semiconductor  circuits  in  hybrid 
architectures.^  Attempts  to  reduce  operating  voltage  have 
been  reported  using  shallow  quantum  wells^  and  superiat¬ 
tice  (SL)  structures,^  and  contrast  ratios  as  high  as  40:1  at 
voltages  as  low  as  4  V  have  been  achieved  at  0.76  /xm  in  the 
AlGaAs/GaAs  system.^  Although  superiattice  diodes  have 
been  fabricated  in  the  transparent  substrate  and  lattice 
matched  InGaAs/InP  system,^  and  photocurrent  spectra 
reported,  no  results  have  appeared  on  the  use  of  such 
structures  for  optical  modulation  in  the  1. 3-1.6  /xm  range. 

In  this  communication  we  report  the  results  of  exper¬ 
iments  designed  to  demonstrate  the  reduction  in  operating 
voltage  that  can  be  achieved  in  InGaAs/InP  p-i-n  modu¬ 
lator  diodes  by  employing  the  Wannier-Stark  effect  in  a 
superiattice  stack**  instead  of  the  more  usual  QCSE  MQW 
structure.  A  semiconductor  superiattice  is  a  stack  of  wells 


and  barrier  layers  thin  enough  so  that  the  wells  are  coupled 
by  resonant  tunneling.  Such  a  modulator  achieves  a  change 
in  the  optical  spectrum  by  misaligning  the  resonant  energy 
levels  in  adjacent,  and  coupled,  quantum  wells  by  means  of 
an  applied  electric  field,  thereby  resulting  in  localization  of 
the  initially  delocalized  electrons  and  holes.  Such  localiza¬ 
tion  occurs  when  the  voltage  between  adjacent  wells  is  on 
the  order  of  the  width  of  the  delocalized  carrier  minibands, 
which  is  typically  -30—60  meV.**  Since  the  SL  period  is 
usually  only  —5—10  nm,  this  means  that  the  necessary 
applied  bias  across  a  1  /xm  stack  will  typically  be  :S  6  V. 

Two  samples,  one  a  MQW  structure  and  the  other  a 
SL,  were  grown  by  gas-source  molecular  beam  epitaxy 
(MBE)  on  n'^  (100)  InP  substrates  at  500 ‘'C.  The  gas 
flow  rates  were  6  seem  of  ASH3  for  an  InGaAs  growth  rate 
of  0.5  pm/h.  and  4  seem  of  PH3  for  InP  growth  at  1.0 
/xm/h.  The  dopants  were  Si  for  n-type  and  Be  for  />-type. 
The  MQW  structure,  consisting  of  49.5  periods  of  undoped 
(w  — 10*^  cm”^)  10  nm  InGaAs  wells  and  10  nm  InP  bar¬ 
riers,  was  grown  on  top  of  a  500-nm- thick  n^  InP  buffer, 
resulting  in  an  active  layer  thickness  of  —  1  /xm.  The  struc¬ 
ture  was  completed  by  growing  a  50-nm-thick  undoped 
InP  layer  and  a  final  1  /xm /?'*"(  1 X 10*^  cm*"^)  InP  contact. 
The  second  sample  was  very  similar  except  that  it  consisted 
of  60.5  periods  of  a  4.2  nm  InGaAs  well/2  nm  InP  barrier 
superiattice,  for  a  total  thickness  of  —0.4  /xm,  less  than 
half  that  of  the  MQW  structure,  and  the  final  InP  p"^ 
contact  layer  was  500  nm  thick.  Such  a  SL  should  exhibit 
a  miniband  width  of  — 120  meV  for  the  first  electron  level 
and  a  resulting  bias  of  —7  V  for  full  carrier  localization. 
Double  crystal  x-ray  diffraction  measurements  showed  that 
the  well  composition  in  the  MQW  sample  was 
Ino  5iGa^.49As  whereas  in  the  SL  structure  Ino  53Gap47As 
was  obtained.  Devices  were  fabricated  in  both  samples  by 
wet  etching  in  1  NK2Cr207:lHBr:lH3P04  to  form  mesas 
of  —550  /xmX550  /xm.  Au-Ti  was  then  e-beam  evaporated 
to  define  contacts  to  both  the  p  and  n  regions.  The  final 
step  was  to  polish  the  substrate  back  surface  and  to  apply 
antireflection  (AR)  coatings  of  Si02  to  both  surfaces.  For 
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FIG.  1.  Photocurrent  spectra  of  a  (a)  MQW  and  (b)  superlatticc  diode 
for  various  values  of  reverse  bias. 

reflection  measurements,  the  AR  coating  on  the  back  sur¬ 
face  was  replaced  by  an  evaporated  A1  mirror.  Electrical 
measurements  confirmed  good  diode  characteristics  with 
reverse  breakdown  voltages  —15  V. 

Figure  1  shows  the  photocurrent  absorption  spectra  of 
the  MQW  and  superlattice  devices,  normalized  to  the  same 
value  at  short  wavelength,  for  various  values  of  reverse 
bias.  In  comparison  to  the  MQW  results  in  Fig.  1(a), 
which  clearly  show  excitonic  resonances  at  zero  bias,  the 
spectra  in  Fig.  1(b),  for  the  SL  device,  indicate  a  lack  of 
any  significant  features  at  0  V,  consistent  with  the  expec¬ 
tation  that  the  carriers  are  unlocalized  and  free  to  move  by 
tunneling  between  adjacent  wells.  Application  of  reverse 
bias  localizes  the  carriers  in  the  wells  rapidly  leading  to  the 
appearance  of  excitons  and  the  staircaselike  absorption 
spectra  characteristic  of  2D  confined  carriers,  the  so-called 
Wannier-Stark  ladder."  The  positions  of  the  exciton  peaks 
in  Fig.  1(a)  at  —1.52  fim  and  in  Fig.  1(b)  at  —  1.48/xm, 
using  a  binding  energy  of  5  meV,  are  consistent  with  the 
alloy  compositions  obtained  from  the  x-ray  data  and  the 
well  widths.  Of  most  significance  for  application  of  these 
devices  is  that  the  rapid  distortion  of  the  optical  spectrum 
for  the  SL  sample  occurs  at  low  voltages.  Specifically,  from 
Fig.  1(b)  we  obtain  a  58%  change  in  absorption  at  1.53 
lim  for  just  4  V  of  reverse  bias  corresponding  to  —80 
kV/cm  maximum  field  in  the  diode.  Transmission  mea¬ 
surements  have  shown  only  a  3.4%  absolute  change  in 
transmission  at  the  exciton  peak  position  of  1.52  }im.  in  the 

,•,.5896  J.  Appl.  Phys..  Voi.  74,  No.  9, 1  November  1993 


FIG.  2.  Tran.?mission  spectra  for  the  SL  device  for  various  values  of 
reverse  bias. 

MQW  case  for  as  much  as  15  V  of  bias,  whereas,  in  Fig.  2 
we  show  the  transmission  spectra  for  the  SL  that  yields  an 
8.1%  change  at  only  4  V  of  bias  at  1.48  /im,  despite  the  fact 
that  the  interaction  length  in  the  SL  is  less  than  half  that  in 
the  MQW.  In  both  cases  the  modulation  is  close  to  linear 
over  these  applied  bias  ranges.  Voltages  ^  5  V  across  the 
SL  resulted  in  saturation  of  the  effect,  corresponding  to  full 
localization  of  the  carriers,  and  the  appearance  of  the 
QCSE  manifest  as  a  shift  in  the  peaks  to  longer  wavelength 
and  a  reduction  in  the  strength  of  the  exciton.  In  contrast, 
modulation  in  the  MQW  devices  continues  to  increase 
quadratically  to  voltages  ^15  V  as  is  shown  in  Fig.  3, 
where  we  have  plotted  the  spectral  distribution  of  the  frac¬ 
tional  modulation  at  various  biases,  defined  as  the  change 
in  signal  on  application  of  bias  divided  by  its  zero  bias 
value,  and  measured  using  a  1  kHz  square  wave  bias  signal. 

While  maintaining  low  operating  voltages,  it  is  well 
known  that  increased  levels  of  modulation  with  little 
change  in  spectral  bandwidth  can  be  achieved  by  operating 
such  modulators  in  reflection,  where  the  signal  passes  twice 
through  the  active  region.  Figure  4  shows  the  reflection 
modulation  characteristics  of  a  SL  device  where  the  rear 
AR  coating  was  replaced  by  an  evaporated  A1  reflector  and 
where  the  total  interaction  length  for  two  passes  through 


FIG.  3.  Fractional  transmission  modulation  measured  at  1  kHz  for  a 
MQW  device. 

Chen  et  ai  5896 


54 


Wavelength  (A) 


FIG.  4.  Fractional  reflection  modulation  measured  at  1  kHz  for  a  SL 
device. 


the  SL  stack  is  --0.8  ptm,  about  the  same  as  for  the  single 
pass  transmission  MQW  modulator  of  Fig.  3.  Maximum 
fractional  modulation  at  4  V  is  32%  at  1.48  ^m,  corre¬ 
sponding  to  a  contrast  ratio  CR  of  1.5:1,  which  should  be 
compared  to  the  10%  change  seen  in  Fig.  3  for  the  MQW 
device  at  this  same  bias.  The  best  reflection  results  we  have 
observed  on  these  nonresonant  mode  SL  devices  at  this 
same  voltage  is  a  CR  of  1.8:1. 

Given  that  the  capacitance  of  these  SL  structures  is 
^15  pF,  we  estimate  a  switching  energy  at  4  V  of  600  pJ 
and  an  associated  dynamic  switching  power  of  99  W/cm^ 


of  the  active  area  for  switching  in  2  ns,  a  value  well  within 
the  generally  accepted  range  for  continuous  nonthermally 
limited  operation."*  Lower  residual  doping  levels  in  the  SL 
would  result  in  a  reduction  in  C  and  a  further  improvement 
in  power  dissipation. 
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A2.2  Gas-Source  MBE 
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A2.2.0  INTRODUCTION 

In  recent  years  the  use  of  gaseous  feed  slocks  to  replace  conventional  solid  sources  in  molecular  beam 
epitaxy  (MBE)  has  provided  a  means  of  greatly  improving  the  control  during  growth  of  III~V  semiconductor 
heterostructures.  Several  variations  of  this  new  technology  have  appeared,  including  gas-source  MBE 
(GSMBE)  where  gaseous  hydrides  are  used  for  the  group  V  elements,  metal-organic  MBE  (MOMBE)  where 
the  metal  alkyls  are  used  for  the  group  III  elements,  and  chemical  beam  epitaxy  (CBE)  where  gases 
are  used  for  both  group  III  and  group  V.  A  review  of  GSMBE,  as  applied  to  the  growth  of  the  III-V 
semiconductors  primarily  containing  phosphorous,  is  presented  here  with  emphasis  on  growth  mechanisms, 
vacuum  technology  and  safety.  A  few  III-V  materials  and  their  associated  heterostructures  will  be 
described  as  illustrations  of  the  GSMBE  technique;  however,  details  of  how  to  grow  a  specific  material 
will  be  left  to  section  E.  A  comprehensive  description  of  GSMBE  and  its  application  to  the  InGaAsP/InP 
materials  system  and  devices  can  be  found  in  the  book  by  Panish  and  Temkin  [1]. 

The  technique  of  GSMBE  is  illustrated  by  the  schematic  diagram  of  figure  A2.2.1.  As  in  conventional 
solid-source  MBE,  a  heated  substrate  in  a  UHV  chamber  is  exposed  to  molecular  beams  of  the  group  III 
elements,  each  beam  supplied  by  a  separate  effusion  oven,  called  a  Knudsen  cell.  Unlike  conventional 
MBE,  however,  the  group  V  molecular  beam(s)  is  produced  by  the  thermal  dissociation  (i.e.  cracking) 
of  the  appropriate  hydride  gas  as  the  gas  enters  the  UHV  growth  chamber.  The  gases  are  usually  arsine 
(ASH3)  and  phosphine  (PH3),  which  at  the  temperatures  of  800-1  lOO'C  necessary'  for  efficient  cracking, 
produce  the  dimer  molecules  As2  and  P2,  respectively.  At  the  usual  MBE  growth  rate  of  1  ^tm  hr“^  a 
significant  amount  of  the  by-product  H2  is  produced  and  is  removed  from  the  growth  chamber  by  suitable 
vacuum  pumps.  Since  arsine  and  phosphine  are  extremely  toxic,  great  care  is  exercised  in  gas  handling 
and  storage  in  order  to  ensure  a  safe  environment  for  the  crystal  grower  and  the  community.  Finally, 
the  phosphors  deposited  in  the  growth  chamber  and  vacuum  pumps  can  produce  hazardous  situations  and 
special  precautions  are  exercised  during  system  maintenance.  It  is  these  features  of  GSMBE  technology  that 
are  unique  and  will  be  discussed  in  more  detail  in  the  next  section  on  equipment  considerations. 

The  physical  mechanisms  controlling  film  growth  in  GSMBE  are  very  similar  to  that  operating  in  solid- 
source  MBE  but  very  different  from  that  in  MOMBE  and  CBE.  For  example,  in  the  growth  of  III-III-V  ternary 
alloys,  with  an  excess  of  the  group  V  flux  present,  the  growth  rate  and  film  composition  in  GSMBE  are 
determined  by  the  rate  of  arrival  of  the  group  III  species  and  largely  independent  of  substrate  temperature. 
Since  it  is  difficult  with  present  technology  to  control  substrate  temperature  in  a  reproducible  manner,  the 
insensitivity  to  run-to-run  variations  in  substrate  temperature  is  a  major  advantage  of  GSMBE.  In  MO.MBE 
and  CBE,  the  rates  of  adsorption,  incorporation,  and  desorption  of  the  metal  alkyl  molecule  are  strongly 
temperature  dependent  and,  thus,  the  film  growth  rate  and  composition  can  be  significantly  affected  by 
a  modest  change  in  substrate  temperature.  Thus  controlled  growth  of  complex  multilayer  heteroepitaxial 
structures  is  much  easier  in  GSMBE  than  in  MOMBE  and  CBE.  Furthermore,  the  carbon  released  by  the 
pyrolysis  of  the  metal  alkyls  can  produce  large  and  undesirable  levels  of  carbon  doping  during  MOMBE 
and  CBE  growth,  while  no  carbon  by-product  is  produced  during  GSMBE.  In  fact,  the  molecular  fragments 
produced  during  cracking  of  ASH3  in  gsmbe  may  act  to  suppress  incorporation  of  low  levels  of  carbon 
(from  CO  and  other  residual  gases  in  the  growth  chamber)  during  growth  [2]. 

In  comparison  to  conventional  MBE  where  solid  AS4  and  P4  are  used  in  heated  effusion  cells,  GS.MBE 
has  several  significant  advantages.  Using  a  gas  source,  fast  and  accurate  control  of  the  group  V  beam 
is  possible.  Varying  the  V/III  ratio  to  check,  for  example,  the  transition  from  an  As-stabilized  to  a  Ga- 
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Figure  A2.2.1.  Schematic  diagram  showing  the 
major  features  of  gas-source  molecular  beam 
epitaxy.  The  group  III  molecular  beams  are 
supplied  by  mbe  effusion  ovens  and  the  group  V 
molecular  beams  are  produced  by  thermal  cracking 
of  the  gases  PH3  and  AsHj. 


stabilized  condition  during  GaAs  growth  takes  just  a  few  minutes,  whereas  changing  the  flux  from  a  large 
capacity,  solid-As4  oven  can  take  hours.  The  use  of  PH3  rather  than  solid  phosphorus  makes  the  MBE 
growth  of  phosphide  compounds  feasible  for  the  reasons  outlined  in  the  next  section.  The  dimer  molecules 
produced  in  GSMBE  have  a  much  larger  sticking  coefficient  than  the  tetramer  molecules  produced  from 
the  elemental  group  V  materials,  and  thus  more  efficient  use  of  source  materials  is  achieved.  Feeding  the 
group  V  source  material  from  outside  the  UHV  chamber  can  significantly  extend  the  time  between  venting 
of  the  growth  chamber  since  group  V  effusion  cells  can  be  rapidly  depleted  at  the  high  V/III  ratios  often 
needed  in  mbe  growth. 

Although  the  above  advantages  are  also  shared  with  MOMBE  and  CBE,  GSMBE  offers  some  unique  advantages 
relative  to  MOMBE  and  CBE.  GSMBE  avoids  the  carbon  incorporation  problem  inherent  in  the  use  of  organic 
source  materials.  It  is  also  advantageous  to  retain  from  conventional  MBE  the  use  of  effusion  cells  for  the 
group  III  sources  since  the  current  generation  of  commercial  effusion  cells  provides  long  life,  excellent 
film  uniformity,  and  adequate  group  III  flux  stability  for  most  applications.  The  use  of  gaseous  feed 
stocks  increases  the  cost  and  complexity  over  that  of  a  solid-source  MBE  system  because  of  the  additional 
equipment  needed  to  control  and  store  the  gases.  However,  a  fully  automated  GSMBE  system  can  be 
implemented  for  substantially  less  cost  than  a  CBE  system. 

An  important  advantage  of  the  MOMBE  and  CBE  techniques  is  the  selective-area  growth  that  results  from 
the  catalytic  decomposition  of  the  metalorganic  compounds  which  can  occur  on  some  surfaces.  The 
resulting  in  situ  patterning  is  a  promising  approach  to  obtain  lateral  device  definition  in  highly  complex 
heteroepitaxial  structures.  GSMBE  does  not  usually  exhibit  growth  selectivity,  although  selective-area 
growth  of  InGaAs/InP  has  recently  been  reported  by  using  GSMBE  with  a  high  group  V  flux  combined 
with  atomic  hydrogen  exposure  [2]. 

We  have  found  that  the  primary  advantage  of  CS.MBE  is  the  capability  to  grow  with  relative  ease  a  wide 
variety  of  III-V  materials  in  one  growth  chamber  and  retain  the  precise  control  of  layer  thicknesses  on 
the  atomic  scale  characteristic  of  elemental  MBE  which  is  absolutely  necessary  for  modem  heteroepitaxial 
devices.  The  primary  disadvantage  of  GSMBE  (and  CBE  as  well)  is  the  extreme  toxicity  of  the  hydrides 
which  results  in  a  substantial  increase  in  the  cost  of  equipment,  installation,  and  operation  and  creates 
a  not  insignificant  amount  of  additional  stress  for  the  crystal  grower.  This  issue  is  expected  to  become 
less  important  in  the  future  as  less  hazardous  methods  of  generating  and  storing  the  hydrides  on-site  are 
developed  and  as  suitable  less  toxic,  alternative  group  V  source  materials  become  available. 
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Vacuum  deposition  using  AsHs  and  PH3  was  first  reported  in  1974  by  Morris  and  Fukui  at  Bell  Laboratories 
who  formed  polycrystalline  layers  of  GaAs,  GaP,  and  GaAsP  [3],  Panish  (1980)  was  the  first  to  report 
MBE  growth  with  the  gases  AsHj  and  PH3  [4],  and  a  year  later  Calawa  at  MIT  Lincoln  Laboratories 
reported  GaAs  with  very  high  electron  mobility  using  ASH3  instead  of  elemental  As  in  MBE  [5].  Panish 
pioneered  the  GSMBE  technique  by  elucidating  the  mechanisms  underlying  thermal  disassociation  of  AsHs 
and  PH3  [4],  by  developing  several  forms  of  hydride  cracking  ovens  [6],  and  by  demonstrating  that  the 
quaternary  alloy  InGaAsP  could  be  readily  grown  using  GSMBE  [7].  GSMBE  has  since  been  used  to  produce 
very  high  quality  layers  of  GaAs  [8]  and  InP  [9],  and  GSMBE  heterostructures  have  been  used  to  fabricate 
a  variety  of  high  performance  devices,  including  InGaAsP/InP  lasers  [10-12],  InGaAs/GaAs/InGaP  lasers 
[13],  InGaAs/InP  bipolar  transistors  [14],  InGaAsP/InP  photodiodes  [15],  and  multiple  quantum  well 
optical  modulators  [16],  The  gases  ASH3  and  PH3  have  been  widely  used  to  produce  films  of  GaAs, 
InP,  GaP,  AIP,  InGaAs,  InGaP,  InGaAsP,  and  InGaAlP,  while  NH3  has  been  used  to  grow  GaN  films 
[17]  and  GaPN  films  [18].  Laboratories  besides  our  own  now  using  gsmbe  to  produce  III-V  films 
include  ATT  Bell  Laboratories  (USA),  British  Telecom  Research  Laboratories  (UK),  CNET  (France), 
Alcatel  Alsthom  Recherche  (France),  Hughes  Research  Laboratories  (USA),  Matsushita  Research  Institute 
(Japan),  Mitsubishi  Kasei  Co.  (Japan),  Sharp  Laboratories  (England),  Sophia  University  (Japan),  Tampere 
University  of  Technology  (Finland),  University  of  California  at  San  Diego  (USA),  University  of  Illinois 
(USA),  and  University  of  Minnesota  (USA)  [19]. 


A2.2.1  EQUIPMENT  CONSIDERATIONS 
Generic  GSMBE  System 

An  MBE  system  designed  for  GSMBE  growths  is  shown  in  figure  A2.2.2  and  consists  of  three  UHV  chambers 
and  a  gas  handling  system.  The  UHV  system  is  a  standard  commercial  MBE  system  (Perkin-Elmer  430) 
with  the  single  3-inch  wafer  capability.  The  introduction  and  preparation  chambers  are  used  for  loading, 
removal,  and  in  situ  storage  of  wafers.  The  N2-purged  glove  box  minimizes  particulate  contamination 
of  the  wafer,  speeds  pump  down  of  the  introduction  chamber  by  reducing  the  noble  gases  incorporated 
in  the  sorption  roughing  pumps,  and  provides  an  extra  degree  of  safety  for  the  operator.  The  growth 
chamber  contains  seven  effusion  cells  and  a  single  gas-cracking  oven  for  both  ASH3  and  PH3,  as  well  as  a 
reflection  high  energy  electron  diffraction  (RHEED)  system  for  in  situ  characterization  of  epitaxial  surfaces. 
A  quadrapole  mass  analyser  (QMA)  is  used  for  determination  of  ASH3  and  PH3  cracking  patterns  as  well 
as  for  leak  checking  and  residual  gas  analysis.  The  gas  handling  system  provides  controlled  introduction 
of  ASH3  and  PH3  at  low  flow  rates  into  the  UHV  growth  chamber  from  high  pressure  gas  storage  cylinders. . 
A  single  computer  provides  integrated  control  of  gas  flows,  effusion  cell  temperatures,  and  shutters.  All 
of  the  effluent  from  the  vacuum  pumps  and  the  gas  handling  system  is  passed  through  a  gas  scrubber  for 
safe  removal  of  toxic  gases. 

The  issues  unique  to  using  the  gaseous  hydrides  in  a  MBE  system,  including  the  choice  of  vacuum  pumps, 
operation  of  the  cracking  oven,  various  schemes  for  hydride  delivery,  and  safety  concerns  are  discussed 
in  the  following  sections. 


Phosphorus  in  MBE 

The  vapour  pressure  of  phosphorus  is  much  higher  than  that  of  arsenic  and  thus  control  of  phosphorus 
in  a  UHV  MBE  system  is  correspondingly  more  difficult.  To  further  complicate  the  situation,  there  many 
allotropic  forms  of  solid  phosphorus,  each  with  a  different  dependance  of  vapour  pressure  on  temperature 
(see  figure  A2.2.3)  [20,21].  White  phosphorus  has  a  very  high  vapour  pressure  (10^  toix  at  200°C),  and 
spontaneously  ignites  in  air  producing  a  dense  white  smoke  that  is  toxic  and  corrosive.  Red  phosphorus 
has  a  lower  vapour  pressure  (IC^  torr  at  200°C),  is  relatively  stable  in  air  at  temperatures  less  than  250°C, 
and  is  less  reactive  than  white  phosphorus. 

Early  MBE  studies  using  solid  red  phosphorus  [22, 23]  in  a  conventional  effusion  cell  produced  a  molecular 
beam  containing  primarily  P4  molecules,  which  have  a  low  incorporation  coefficient  at  the  substrate 
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Figure  A2.2.2.  A  gsmbe  system  with  a  single  growth  chamber,  two  auxiliary  UHV  chambers, 
and  a  gas  handling  system. 
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Figure  A2.2.3.  The  vapour  pressures  as  functions  of 
temperature  for  solid  AS4  and  the  allotropic  forms  of  black,  red 
and  white  P4  (based  on  the  vapour  pressure  curves  summarized 
in  [20,21]). 


Hnn^hnnlr  nfThin  Film  PrOCCSS  TechnolORV 


§0 


©  1995  lOP  Publishing  Ltd 


Gas^Source  MBE 


temperatures  used  in  MBE.  Thus  very  high  P4  fluxes  are  needed  for  growth  and  substantial  accumulations 
of  phosphorus  are  deposited  on  the  walls  of  the  MBE  chamber.  Furthermore,  the  wall  deposits  are  a  mixed 
phase  of  red  and  white  phosphorus,  with  the  amount  of  each  varying  from  run  to  run.  Reproducible  growth 
is  very  difficult  to  achieve  because  of  the  varying  background  phosphorus  pressure,  and  the  large  amount 
of  white  phosphorus  greatly  complicates  safe  procedures  for  chamber  cleaning.  In  addition,  Panish  has 
pointed  out  that  the  molecular  flux  obtained  from  heating  red  phosphorus  is  expected  to  depend  on  the 
thermal  history  of  the  phosphorus  effusion  cell,  and  thus  achievement  of  precise  control  of  the  phosphorus 
flux  is  unlikely  [1]. 

On  the  other  hand,  employing  a  source  of  P2  molecules  can  be  a  practical  means  of  producing  phosphide 
films  in  MBE.  The  P2  species  exhibits  a  much  higher  incorporation  coefficient  S\  than  the  P4  species, 
and  thus  requires  a  much  lower  phosphorus  beam  flux.  We  have  found  that  for  GSMBE  of  InP  at  520°C 
and  1  /imhr*‘^  si  is  about  0.2-0.3  for  P2,  which  is  about  50  times  larger  than  si  for  growth  with  P4. 
Furthermore,  P2  vapour  condenses  primarily  as  red  phosphorus  on  the  surfaces  of  the  growth  chamber, 
producing  a  much  lower  phosphorus  background  pressure  [20].  This  permits  reproducible,  controlled 
epitaxial  growth  of  phosphide  films  and  the  growth  of  arsenide  films  (i.e.,  GaAs,  InGaAs)  with  low 
phosphorus  contamination,  in  the  same  growth  chamber.  The  smaller  accumulation  of  white  phosphorus 
deposits  extends  the  time  between  venting  of  the  system  for  maintenance  and  reduces  the  fire  hazard  during 
chamber  cleaning. 

A  P2  molecular  beam  can  be  produced  by  several  methods.  Wright  and  Kroemer  [24]  have  used  heated 
GaP  to  provide  a  P2  beam  for  MBE  growth,  but  the  beam  contained  Ga  and  the  P2:Ga  flux  ratio  varied 
with  use.  In  GSMBE,  phosphine  is  thermally  decomposed  using  special  gas  cracking  ovens,  described  in  a 
following  section,  to  produce  a  well  controlled  P2  beam.  In  solid-source  MBE,  a  valved  effusion  oven  for 
generation  of  a  As2  beam  from  solid  AS4  has  recently  been  developed  [25],  and  versions  of  this  cracker 
may  make  the  generation  of  a  P2  molecular  beam  from  solid  phosphorus  practical  [26, 27],  A  large  canister 
of  solid  red  phosphorus,  mounted  outside  the  MBE  chamber,  is  heated  and  the  resulting  P4  vapour  is  passed 
through  a  needle  valve  and  then  thermally  cracked  at  a  higher  temperature  to  produce  a  P2  beam.  The 
valve  controls  the  P2  flux  while  the  canister  is  kept  at  a  constant  temperature.  With  proper  design,  the  valve 
can  be  kept  free  of  clogging,  and  III-V  phosphide  films  can  be  grown  [27,28].  The  primary  advantage  of 
the  valved  cracking  oven  is  safety;  the  highly  toxic  PH3  is  avoided.  Although  high  efficiencies  have  been 
reported  for  a  valved  cracker  [25],  the  efficiency  for  creation  of  dimer  molecules  from  tetramer  molecules 
in  a  solid-source  cracker  may  not  as  high  as  for  the  hydride-gas  cracker  discussed  below.  We  have  found 
the  AS2/AS4  ratio  to  be  about  6-8  times  higher  for  our  hydride-gas  cracker  than  for  a  solid-source  cracker 
measured  under  similar  conditions. 


Vacuum  Pumping 

During  growth,  the  As2  and  P2  species  are  readily  pumped,  because  of  their  high  vapour  pressures,  by  the 
liquid-nitrogen-cooled  panels  mounted  on  the  walls  of  the  growth  chamber.  Thus  the  primary  gas  load 
for  the  vacuum  pumps  is  the  large  amount  of  H2  produced  during  hydride  cracking.  For  a  typical  III-V 
layer,  the  pressure  at  the  sample  should  be  less  than  1  x  10“"^  torr  for  molecular  flow  conditions,  and  for 
a  film  growth  rate  of  1  fim  hr"\  a  gas  flow  rate  of  about  2-10  seem  is  needed,  depending  on  the  design 
of  the  growth  chamber.  This  will  produce  a  gas  load  of  up  to  6  litres  of  H2  in  a  10  hour  deposition.  The 
pressure  at  the  nozzle  of  the  gas  cracking  oven  will  be  about  1  x  10”^  torr,  where  the  pumping  rate  needed 
is  about  130  1  s"^  This  rate  translates  to  a  pumping  speed  about  1000  1  s""*  or  more  at  the  vacuum  pump, 
depending  on  the  configuration  of  the  growth  chamber.  To  maintain  lower  pressures  or  to  operate  at  higher 
growth  rates  requires  even  higher  pumping  rates.  Also,  the  group  V  element  is  preferentially  desorbed 
from  a  lU-V  surface  at  elevated  temperatures  and,  thus,  the  maximum  substrate  temperature  is  determined 
by  the  rate  at  which  the  group  V  molecules  can  be  supplied  to  replace  the  lost  group  V  surface  atoms. 
Hence,  the  pumping  rate  limits  the  maximum  substrate  temperature  as  well  as  the  maximum  growth  rate. 

Three  different  types  of  UHV  pumps — diffusion,  cryogenic  and  turbomolecular — have  all  been  used 
successfully  alone  or  in  combination  in  GSMBE.  The  diffusion  pump  offers  high  pumping  rates,  low 
cost  and  excellent  reliability,  and  large  diffusion  pumps  have  been  used  in  GSMBE  commercial  systems 
from  Vacuum  Generators  [1]  and  Perkin-EImer  [29]  to  grow  III-V  films  of  high  quality.  The  major 
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disadvantages  of  the  diffusion  pump  are  the  cold  trap  which  can  greatly  reduce  the  pumping  speed  and 
the  possible  contamination  of  the  growth  chamber  by  the  pumping  fluid.  With  proper  operation,  carbon 
contamination  from  the  diffusion  pump  can  be  insignificant.  Using  a  diffusion-pumped  GSMBE  system, 
Stillman  and  co-workers  have  obtained  GaAs  with  a  residual  free  carrier  concentration  as  low  as  9  x  10'^ 
cm“^  [29],  which  is  comparable  to  highest  purity  GaAs  grown  by  conventional  solid-source  MBE. 

Turbomolecular  pumps  have  high  pumping  rates  for  Ha  and  the  newer  turbomolecular  pumps  are  chemically 
resistant  to  the  metal  alkyls  used  in  MOMBE.  However,  turbomolecular  pumps  are  very  expensive  and  failure 
of  the  pump  can  be  catastrophic.  If  metalorganic  gases  are  to  be  used,  turbomolecular  pumps  are  best 
suited  to  the  task,  but  if  the  only  gases  to  be  used  are  the  hydrides,  then  either  diffusion  or  cryogenic 
pumps  will  work  well. 

Closed  cycle  He-cooled  cyropumps  are  well  suited  for  GSMBE,  since  they  can  provide  high  pumping  rates 
for  Ha  at  the  required  pressures  through  the  use  of  activated  charcoal  to  adsorb  hydrogen.  Cyropumps 
cost  substantially  less  than  turbomolecular  pumps  and  produce  little  or  no  contamination  to  the  growth 
chamber.  But  since  the  pumped  gas  is  retained,  the  pump  must  be  frequently  regenerated.  Figure  A2.2.4 
shows  our  GSMBE  growth  chamber  with  two  CTI  model  CT-8  cryopumps,  each  with  2200  1  s  ’  (Ha) 
pumping  rate,  that  share  a  single  compressor.  The  pumps  are  mounted  close  to  the  growth  region  in  order 
to  maximize  the  available  pumping  speed  and  in  a  horizontal  position  to  minimize  particle  collection  on 
the  seat  of  the  UHV  gate  valves.  The  pressure  relief  valves  on  the  CT-8  pumps  are  vented  directly  to 
a  gas  scrubber.  A  mechanical  pump  containing  a  fluorocarbon  fluid  (Fomblin)  and  flushed  with  N2  is 
used  for  regeneration  of  the  He  cryopumps.  No  in-line  traps  are  used  and  no  evidence  of  phosphorus 
accumulation  in  the  roughing  pump  or  lines  has  been  observed  in  over  600  days  of  growth.  We  find  we 
can  grow  GaAs  or  InP  at  1  /zm  hr"'  using  a  new  CT-8  pump  for  3-5  days  of  growth  before  regeneration 
is  needed.  However,  for  the  growth  of  InGaAsP  alloys  where  both  As2  and  P2  beams  are  used  at  high 
fluxes,  daily  regeneration  is  required.  Two  pumps  then  allow  uninterrupted  growth  by  alternating  pumps 
between  growth  and  regeneration. 

We  have  found  that  the  H2  capacity  of  the  CT-8  cryopump  degrades  considerably  with  use  in  our  GSMBE 
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Figure  A2.2.4.  Vacuum  pumping  system  for  gsmbe  growth.  One  cryopump  is  used 
during  growth  while  the  other  can  be  regenerated.  The  ion  pump  is  used  only  during  idle 
periods. 
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system.  As  shown  in  figure  A2.2.5,  we  have  measured  the  pump  capacity  by  noting  the  H2  dose  where 
the  pressure  suddenly  increases  due  to  pump  regurgitation.  The  tests  were  conducted  at  a  constant  flow 
rate  of  32  seem  using  Ha,  not  ASH3  or  PH3,  through  the  gas  cracker,  and  thus  do  not  duplicate  the  exact 
conditions  present  during  growth.  However,  it  is  apparent  that  the  capacity  decreases  markedly  with  use, 
dropping  from  over  1 1  standard  litres  of  H2  to  under  5  litres  in  about  24  months  of  use.  Replacement  of  the 
activated  charcoal  restores  the  pumping  capacity.  Apparently,  phosphorus  is  accumulating  on  the  surfaces 
of  the  charcoal  and  inhibiting  H2  pumping.  The  replacement  of  the  activated  charcoal  is  a  straightforward 
procedure,  costs  about  5%  of  a  new  cryopump,  and  does  not  require  opening  the  growth  chamber.  Note 
that  care  must  exercised  when  growing  with  a  nearly  full  cryopump.  The  pressure  rises  very  rapidly  when 
the  pump  reaches  its  current  capacity,  which  can  result  in  tripping  of  pressure-sensitive  interlocks  that 
control  the  main  power  to  the  MBE  system  and  damage  to  the  effusion  cells  during  the  ensuing  abrupt  cool 
down. 


Hydrogen  Capacity  of  Cryopump 


Figure  A2.2,5.  Measurement  of  the  hydrogen  capacity  of  a  CTI  CT-8 
cxyopump  with  increasing  use  in  a  GSMBE  system.  The  time  at  which  the 
pressure  abruptly  increases  is  a  measure  of  the  pump  capacity. 


Hydride  Delivery  Systems 

The  flow  of  the  hydride  gases  can  be  controlled  using  either  of  two  methods;  mass  flow  control  or  pressure 
control.  Mass  flow  control  systems  use  mass  flow  controllers  based  on  the  temperature  change  that  occurs 
in  a  gas  flowing  through  a  restricted  volume.  This  approach  offers  direct  control  of  gas  flow  rates  and 
precise  mixing  of  AsHs  and  PH3  for  growth  of  such  materials  as  GaAsP  or  InGaAsP.  The  currently 
available  mass  flow  controllers  lack  adequate  reproducibility  at  the  low  flow  rates  (about  1-2  seem)  used 
in  GSMBE  and  thus  can  require  frequent  recalibration.  However,  mass  flow  control  systems  are  relatively 
simple  to  implement  and  have  been  used  in  several  commercial  GSMBE  systems. 

Pressure  control  systems  employ  a  precision  leak  valve  to  adjust  the  pressure  of  a  gas  which  is  flowing 
through  a  fixed  aperture.  A  pressure  transducer  (usually  a  capacitance  manometer)  and  a  servo-valve  are 
used  to  obtain  a  highly  reproducible  flow  rates  over  a  very  wide  dynamic  range,  with  very  good  control 
at  low  flow  rates.  Calibration  of  the  flow  rate  in  terms  of  the  indicated  pressure  is  required  for  each  gas, 
since  the  flow  rate  will  depend  on  the  conductance  of  the  gas  in  the  manifold  and  cracking  oven.  Pressure 
control  works  extremely  well  for  GSMBE  and  has  also  been  employed  in  MOMBE  for  precise  control  of  low 
vapour  pressure  gases  without  the  need  for  a  carrier  gas  [30]. 

Figure  A2.2.6  is  a  schematic  diagram  of  the  pressure-control  hydride  delivery  system  in  use  in  our 
laboratory.  Only  the  components  for  PH3  flow  control  are  shown;  an  identical  system  for  ASH3  flow 
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control  is  mounted  in  the  same  cabinet  This  system  is  based  on  the  work  of  Panish  and  Sumsld  at 
Bell  laboratories  [7],  designed  in  our  laboratory,  and  assembled  by  Precision  Flow  Devices  (Beaverton, 
Oregon).  The  run/vent  valve  coniiguradon  minimizes  flow  transients  when  the  P2  (or  Asa)  molecular  beam 
is  switched  on  or  off  by  keeping  the  pressure  constant  against  the  conductance  of  the  gas  cracking  oven  in 
the  growth  chamber  downstream  of  the  run  valve  and  the  gas  tubing  down  stream  of  the  vent  valve.  The 
cooled  sorption  pump  provides  a  low  pressure  to  mimic  the  growth  chamber  and  adequate  capacity  for 
temporary  storage  of  any  uncracked  hydride.  At  the  end  of  the  growth  day,  the  sorption  pump  is  allowed 
to  warm  up  and  the  effluent,  diluted  with  N2,  is  pumped  into  a  gas  scrubber.  In  this  system  designed  to 
operate  with  a  low  pressure  gas  cracking  oven,  an  indicated  gas  pressure  of  10  torr  corresponds  to  a  flow 
rate  of  approximately  10  seem.  The  flow  rate  can  controlled  with  a  precision  of  ±0.02  seem. 


Figure  A2.2.6.  Schematic  diagram  of  a  hydride  delivery  system  for  GSMBE  growth.  A 
computer  controls  the  opening  and  closing  of  the  mn  and  vent  valves  as  well  as  the  gas  flow 
rate  via  the  precision  leak  valve.  All  components  which  contain  toxic  gas  at  pressures  above 
atmospheric  pressure  are  enclosed  in  the  vented  cabinet. 


The  growth  of  complex  III-V  heterostructures  with  many  layers  of  different  composition  usually  requires 
the  hydride  gases  to  be  switched  on  and  off  at  each  heterointerface.  In  order  to  achieve  rapid  gas 
switching,  the  volume  of  the  gas  lines  down  stream  of  the  run  valve  in  figure  A2.2.6  should  be  small. 
Also  the  conductance  of  the  gas  tubing,  which  depends  on  the  pressure  as  well  as  the  mbe  diameter  and 
length,  should  be  large  to  allow  rapid  removal  of  the  residual  gas  via  the  MBE  chamber  after  the  run 
valve  closes  and  rapid  charging  of  the  gas  lines  after  the  run  valve  opens.  Furthermore,  rapid  switching 
requires  the  gas  cracking  oven  not  to  impede  gas  flow.  Using  the  low  pressure  gas  cracker  described 
below,  our  system  exhibits  a  gas  switching  transient  of  about  4-6  seconds;  other  GSMBE  systems  have 
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achieved  gas  switching  times  of  less  than  1  second  [1,31].  It  should  be  noted,  however,  that  in  order  to 
obtain  atomically  abrupt  interfaces  of  InGaAs/InP  or  InGaP/GaAs,  a  growth  pause  of  10  or  more  seconds 
is  required  at  each  interface  [32-35],  Hence,  faster  gas  switching  times  are  not  usually  needed. 

Pressure  control  can  be  readily  used  to  grow  quaternary  alloys  of  In;tGai«;cAs3,P]-y  with  excellent  control 
of  alloy  composition.  As  shown  in  figure  A2.2.7(a),  the  gases  AsHs  and  PH3  are  mixed  before  entering  the 
MBE  system.  The  flow  rate  F  of  each  gas  is  controlled  by  adjusting  the  pressure  P  against  the  conductances 
G  and  kG  of  the  tubing  in  the  gas  mixing  manifold.  Under  normal  operating  conditions,  the  total  flow 
rate  F\  +  F2  is  approximately  a  linear  function  of  the  total  inlet  pressure  P\  +  P2.  as  long  as  the  pressure 
at  node  N  is  larger  than  the  smaller  of  Pi  or  P2.  On  the  other  hand,  the  flow  ratio  P1/P21  which  controls 
the  group  V  composition  of  the  epitaxial  film,  is  a  nonlinear  function  of  the  pressure  ratio  P2/P1.  For  our 
gas  manifold,  the  conductance  ratio  k  is  about  2.0  and  the  calculated  dependance  of  P1/F2  on  P1/P2  is 
shown  in  figure  A2.2.7(b)  along  with  the  operating  points  for  three  different  InGaAsP  compositions.  With 
this  system,  we  can  easily  control  flow  ratios  as  high  as  15:1  with  better  than  ±0.5%  reproducibility  at 
a  total  flow  rate  of  10  seem.  The  maximum  ratio  of  AS2/P2  needed  at  the  surface  of  the  InGaAsP  was 
found  to  only  be  about  6:1  under  typical  growth  conditions.  High  quality  InGaAsP  layers  with  band  gap 
emission  at  wavelengths  from  1.15  to  1.67  mm  and  lattice  matched  to  InP  substrates  have  been  obtained 
in  a  reproducible  manner  using  this  hydride  delivery  system  [36]. 


(b) 

Figure  A2JZ.7,  Pressure  control  of  hydride  gas  flow  in  gsmbe.  (a) 
Schematic  diagram  of  gas  manifold,  (b)  Flow  ratio  as  a  function  of 
pressure  ratio  for  a  manifold  with  gas  conductance  ratio  of  k  =  2.0. 
The  operating  points  for  three  different  InGaAsP  compositions  are 
noted  for  PH3  as  gas  1  and  AsHs  as  gas  2. 
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Hydride  Cracker 

It  is  found  that  the  hydrides  do  not  decompose  on  m-V  surfaces  during  MBE  growth,  even  though  the 
substrate  temperatures  are  typically  500-600“C.  Thus,  the  hydrides  are  thermally  cracked  in  specially 
designed  high  temperature  gas  cracking  oven  (i.e.  cracker)  mounted  in  the  MBE  chamber  in  place  of 
a  conventional  effusion  cell.  The  ideal  cracker  should  reduce  the  hydride  to  a  beam  containing  only 
elemental  group  V  species  and  hydrogen  in  order  to  avoid  introduction  of  the  volatile,  toxic  hydrides  into 
the  vacuum  system.  As  discussed  above,  the  cracker  should  also  produce  dimer  molecular  species,  since 
the  incorporation  coefficient  of  the  dimers  Asa  and  P2  are  much  higher  than  the  tetramers  AS4  and  P4. 
Furthermore,  the  use  of  P2  minimizes  white  phosphorus  accumulation  in  the  growth  chamber. 

Two  general  types  of  crackers  have  evolved  for  use  in  GSMBE.  Panish  first  used  a  thermal  cracker  operating 
at  inlet  pressures  of  200-1500  torr  to  thermally  decompose  AsHs  and  PH3  at  temperatures  of  900- 
lOOO^C  [6].  The  cracker  contained  a  high  pressure  and  a  low  pressure  region,  separated  by  a  carefully 
constructed  flow  constriction  in  AI2O3  tubing  as  shown  in  figure  A2.2.8(a).  By  examining  the  equilibrium 
thermodynamics  for  hydride  decomposition,  Panish  showed  that  in  the  high  pressure  region,  both  ASH3  and 
PH3  should  readily  form  the  tetramer  molecules  AS4  and  P4,  but  dimerization  required  the  lower  pressure 
region.  Both  processes  take  place  efficiently  at  900-1200®C.  A  version  of  the  high  pressure  cracker  was 
manufactured  by  VG  Semicon. 


Figure  A2JS.8.  Gas  sources  for  gsmbe.  (a)  High-pressure  gas  cracking  oven,  (b)  Low-pressure  gas 
cracking  oven.  (Reprinted  with  permission  from  [58]). 


The  second  and  more  widely  used  type  of  cracker  is  the  low  pressure  cracker.  Here  the  inlet  pressure 
varies  from  2  mtorr  to  10  torr,  depending  on  the  design  of  the  cracking  region.  Early  low  pressure  CTackers 
were  constructed  of  tantalum,  since  several  workers  demonstrated  that  heated  Ta  acts  as  a  catalysis  in  the 
decomposition  of  ASH3  and  PH3  [5,7,37].  The  use  of  Ta  baffles  in  the  form  of  crushed  foil  and  pellets 
can  increase  the  efficiency  of  cracking  and  dimerization,  albeit  at  the  expense  of  flow  conductance  and  gas 
switching  speed.  A  low  pressure  cracker  of  this  type  is  shown  in  figure  A2.2.8(b).  We  have  successfully 
used  a  low  pressure  all-Ta-constructed  cracker  supplied  by  Ulvac  in  Japan  (model  ZB61-  9610)  for  GSMBE. 
The  cracker  consists  of  a  heated  region  with  11  Ta  disks  each  with  several  small  holes,  inside  a  12  cm 
long,  1.5  cm  diameter  Ta  tube.  The  cracker  operates  at  an  inlet  pressure  of  2-10  torr  and  at  temperatures  of 
8(X)-1000°C.  At  cracker  temperatures  above  1(K)0°C,  the  total  ionized  impurity  level  is  found  to  increase 
in  the  sample,  apparently  as  a  result  of  contamination  from  the  hot  surfaces  in  the  cracker.  Below  800  C, 
the  cracking  efficiency  falls  off  rapidly,  as  evidenced  by  AsH  and  AsHz  molecular  species,  detected  by 
QMA  measurements,  produced  by  the  incomplete  pyrolysis  of  ASH3.  We  have  found  that  the  molecular 
fragments  created  by  partially  decomposed  ASH3  may  reduce  the  incorporation  of  residual  carbon  during 
GS^E  of  GaAs  [38].  However,  the  advantage  of  a  lower  carbon  contamination  in  the  sample  is  greatly 
offset  by  the  increased  risk  of  AsHs  accumulation  in  the  growth  chamber  during  long  periods  of  growth 
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with  the  cracker  temperature  below  8(X)®C.  Although  an  early  study  [37]  indicated  that  the  Ta  catalytic 
action  degraded  with  use  through  the  formation  of  tantalum  phosphide  on  the  surfaces  of  the  cracker,  in 
over  six  years  of  operation  with  our  Ulvac  cracker  using  both  AsHs  and  PH3,  we  have  observed  no  change 
in  the  efficiency  of  cracking  or  dimerization. 

In  order  to  increase  the  conductance  of  the  low  pressure  cracker,  recent  cracker  designs  have  reduced  or 
removed  the  baffles  without  appreciable  loss  of  cracking  efficiency.  Figure  A2.2.9  compares  the  cracking 
efficiency  r?cr  and  dimerization  efficiency  r/dim  for  Ta-  and  W-based  high  conductance  hydride  crackers 
used  by  Jackson  and  co-workers  at  the  University  of  Illinois  [31].  The  efficiencies  were  measured  using 
a  QMA  to  determine  the  relative  amounts  of  each  molecular  species  emanating  from  the  cracker.  Both 
crackers  approach  100%  efficiency  at  high  cracking  temperatures,  but  the  Ta-based  cracker  can  operate 
about  200°  C  lower  in  temperature,  indicating  that,  even  in  a  cracker  with  a  straight  tube  without  baffles, 
the  catalytic  action  of  Ta  can  be  very  effective.  Other  high  conductance  crackers  appear  to  use  large-area, 
inert  surfaces  without  Ta  to  achieve  useful  cracking  efficiencies,  such  as  the  proprietary  design  by  the 
Intevac  Corporation  which  employs  Mo  construction.  In  this  type  of  cracker,  an  appreciable  concentration 
of  monomers  may  be  produced  [1].  Recently,  As  monomers  were  found  to  be  more  readily  incorporated 
than  As2  which  in  turn  was  more  readily  incorporated  than  AS4,  in  the  GS.MBE  growth  of  GaAS;cPi~x  alloys 
[39].  Thus  monomer  generation  via  cracked  AsHs  may  prove  to  be  very  useful. 


Gas  Cracker  Temperature  (°C) 


Figure  A22S.  Comparison  of  cracking  and  dimerization  efficiencies  for 
Ta-based  and  W-bas^  crackers  using  ASH3.  Dashed  curves  are  for 
dimerization  efficiency  and  solid  curves  are  for  cracking  efficiency.  (Reprinted 
with  permission  from  [31]). 


It  should  be  noted  that  accurate  measurement  of  the  value  of  the  cracking  efficiency  is  fraught  with 
difficulties  since  the  QMA  will  detect  molecular  species  reflected  from  hot  surfaces  in  the  growth  chamber, 
the  QMA  can  be  nonlinear  at  the  pressures  used  in  growth,  and  the  QMA  can  crack  the  molecules  before 
they  are  detected.  The  comparison  of  the  cracking  behaviour  in  figure  A2.2.9  is  meaningful  since  the 
measurements  were  carried  out  in  an  identical  fashion,  but  the  value  of  the  efficiencies  given  may  be  in 
error. 

In  order  to  determine  the  optimum  gas  switching  performance  of  a  gas  cracking  oven  and  gas  handling 
system,  Jackson  et  al  [31]  have  also  have  examined  the  composition  depth  profiles  of  InP/InGaAs 
heterostructures  grown  under  different  gas  flow  conditions.  The  heterostructures  were  grown  without 
growth  pause  at  the  interfaces  in  a  GSMBE  system  fitted  with  a  5000  Is”*  He  cryopump  and  a  2200  I  s”* 
turbomolecular  pump.  Using  secondary  ion  mass  spectroscopy  (SIMS),  the  composition-depth  profiles  of 
figure  A2.2.10  were  obtained  for  two  different  growths.  In  figure  A2.2. 10(a),  a  single  low  pressure  cracker 
with  Ta  baffles  was  used  for  both  AsHs  and  PH3  and  the  gas  switching  manifold  was  located  about  1  m 
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from  the  cracker  inlet  port.  A  long  P  tail  is  seen  to  extend  into  the  InGaAs  layer  and  the  As  appears  to 
penetrate  into  the  InP  layer  as  well.  In  figure  A2.2.10(b),  separate  high  conductance  Ta  crackers  without 
baffles  were  used  with  the  run/vent  switching  valves  mounted  at  the  gas  inlet  flange  on  the  MBE  growth 
chamber.  The  layer  thicknesses  and  growth  conditions  were  the  same  as  for  the  sample  of  figure  A2.2. 10(a). 
The  P  concentration  in  the  InGaAs  has  dropped  to  the  detection  limit  of  the  SlMS  instrument  and  the  P 
and  As  tails  have  been  greatly  reduced.  Clearly,  in  this  GSMBE  configuration,  significant  improvement  in 
heterointerface  abruptness  was  achieved  by  using  high  conductance  hydride  crackers,  reducing  the  volume 
between  the  gas  switching  manifold  and  the  cracker,  and  utilizing  separate  flow  paths  for  AsHs  and  PH3. 


Profile  Depth  (pm) 


Figure  A2^.10.  SIMS  depth  profiles  for  InGaAs/InP 
heterostructures  grown  by  gsmbe.  (a)  Grown  with  a  single 
low  conductance  hydride  gas  cracker  and  remotely  mounted 
gas  switching  manifold,  (b)  Grown  with  separate  high 
conductance  gas  crackers  and  closely  mounted  gas  switching 
manifold.  (Reprinted  with  permission  from  [31]). 


System  Operation 

The  placement  of  the  eight  sources  on  the  growth  flange  in  our  GSMBE  is  shown  in  figure  A2.2.11.  The 
dopant  cells  for  Be  and  Si  have  a  maximum  capacity  of  2  cc  while  the  remaining  effusion  cells  have  a 
maximum  capacity  of  60  cc  each.  The  inclination  of  the  axis  of  each  cell  from  the  horizontal  varies  from 
less  than  10  degrees  near  the  top  of  the  flange  to  46  degrees  at  the  bottom.  In  and  Ga  are  used  at  the 
highest  rate  and  thus  the  cells  for  these  elements  are  positioned  near  the  bottom  of  the  flange  to  provide 
to  largest  usable  volume  of  the  cell  crucibles  without  spillage.  Two  In  and  two  Ga  cells  are  used  to  allow 
growth  of  at  least  two  InGaAsP  layers  of  different  group  III  composition  in  the  same  growth  run  without 
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changing  cell  temperatures.  The  placement  of  the  cells  are  chosen  to  minimize  the  lateral  variation  of  the 
group  ni  alloy  composition  across  the  wafer. 


Rgure  A2.2.11.  Placement  of  effusion  ovens  and  gas 
cracker  on  gsmbe  source  flange. 


Each  effusion  cell  exhibits  a  momentary  overshoot  followed  by  a  slow  decrease  in  the  molecular  flux  when 
the  mechanical  shutter  is  opened  because  the  melt  surface  temperature  will  drop  slightly  due  to  radiation 
loss  after  shutter  opening.  The  flux  transient  depends  on  cell  construction  as  well  as  on  the  interval  for 
which  the  shutter  has  been  closed.  The  In  and  Ga  cells  are  used  in  pairs  to  match  the  flux  transients  and 
thus  minimize  changes  in  layer  composition  with  depth  at  the  start  of  the  growth  of  alloys  such  as  InGaAs 
or  InGaAsP.  However,  the  growth  rate  will  vary  with  time  and  in  the  growth  of  the  thin  layers  required 
in  quantum  well  (QW)  heterostructures,  the  flux  transients  can  significantly  alter  the  layer  thickness  from 
that  for  a  constant  growth  rate  [40].  The  placement  of  the  cells  is  chosen  to  minimize  the  lateral  variation 
of  alloy  composition  across  the  finished  wafer.  The  placement  of  the  gas  cracker  is  not  critical,  since  the 
As2  and  P2  beams  emanate  from  the  same  location  on  the  source  flange  when  using  a  single  gas  cracker. 

During  a  typical  GSMBE  growth  the  chamber  pressure  is  (2-8)xl0~^  torr,  which  is  primarily  the  partial 
pressure  of  hydrogen.  Chilled  water  is  used  to  cool  the  shroud  surrounding  the  effusion  cells  and  liquid 
nitrogen  is  used  to  cool  the  shrouds  on  the  walls  of  the  growth  chamber.  The  hot  filaments  in  the  effusion 
cells,  RHEED,  QMA,  and  ionization,  gauges  operate  in  the  normal  fashion.  To  prolong  ionization  gauge  life, 
the  emission  current  for  the  gauge  should  be  kept  at  its  lowest  setting.  The  useful  operating  life  of  the 
electron  multiplier  in  the  QMA  is  about  six  months.  We  find  the  current  generation  of  commercial  effusion 
cells  are  long  lived  and  provide  a  stable  and  reproducible  flux  even  in  the  H2  ambient  of  GSMBE.  No 
adverse  effects  on  the  life  of  the  oven  heater  filaments  have  been  observed,  unlike  in  MOMBE  where  the 
metal  alkyls  can  rapidly  corrode  and  destroy  hot  filaments.  However,  we  have  found  a  white  crast  forms 
on  the  surface  of  the  Be  source  material  over  a  period  of  months  in  our  GSMBE  system.  The  crust  will 
eventually  prevent  the  sublimation  of  the  Be  and  the  cell  is  no  longer  useable.  The  Be  can  be  removed, 
the  Be  metal  etched  in  HF:H20  1:100,  and  then  replaced  without  the  loss  of  the  expensive,  high  purity 
Be  material.  We  speculate  that  the  crust  is  a  result  of  the  presence  of  phosphorus  in  the  growth  chamber, 
but  the  composition  of  the  crust  has  not  been  identified. 

Commercial  phosphine  is  currently  available  with  total  impurity  levels  in  the  range  2-5  ppm.  In  some 
applications,  higher  levels  of  purity  are  needed  and  on-site  purification  can  be  used  to  improve  the  properties 
of  GSMBE  films  grown  with  commercial  phosphine  [41].  We  have  used  a  gas  filter  consisting  of  synthetic 
zeolite  molecular  sieve  to  purify  100%  PH3  prior  to  thermal  cracking  during  the  growth  of  InP,  InGaP, 
and  InGaAlP.  The  filter  contained  150  cc  of  3A  molecular  sieve,  which  has  a  molecular  structure  that  is 
particularly  effective  in  trapping  water.  SIMS  data  indicated  that  a  major  source  of  oxygen  contamination 
was  the  as-supplied  PH3  gas.  The  oxygen  had  no  measurable  effect  on  the  InP  and  very  little  effect  on  the 
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InGaP.  However,  the  addition  of  only  1%  A1  to  form  InGaAlP  had  the  dramatic  effect  of  greatly  increasing 
the  oxygen  content  of  the  epitaxial  layer.  The  A1  was  veiy  effective  in  aiding  the  incorporation  of  oxygen 
in  the  InGaAlP  epitaxial  surface  during  growth.  The  combination  of  the  A1  and  oxygen  led  to  a  significant 
decrease  in  electron  concentration  in  n-type  InGaAlP,  producing  non-conducting  material.  However,  the 
molecular  sieve  was  found  to  reduce  the  oxygen  content  and  permitted  the  growth  of  conducting  InGaAlP. 
No  tests  were  conducted  to  determine  the  effectiveness  of  on-site  purification  on  ASH3. 

At  end  of  a  day  of  growth,  the  temperatures  of  the  effusion  cells  and  the  gas  cracker  are  reduced  to  200°C, 
with  a  slow  cooling  rate  of  about  10°C  min"'  used  for  the  A1  cell  to  avoid  breakage  of  the  crucible.  The 
liquid-nitrogen-cooled  shrouds  are  allowed  to  warm  up  to  room  temperature,  which  produces  a  pressure 
increase  in  the  growth  chamber  from  1  x  10"’  torr  to  about  1  x  10"^  torr  in  a  period  of  about  one  hour. 
Within  about  five  hours  the  pressure  has  dropped  to  about  1  x  10"*  torr,  mostly  due  to  residual  P2,  and  the 
system  is  ready  for  another  day  of  growth.  By  cycling  the  shrouds  to  room  temperature  between  growth 
days,  the  consumption  of  liquid  nitrogen  is  greatly  reduced.  The  total  volume  of  gas  evolved  from  the  wall 
shrouds  is  1-2  si,  which  can  be  a  significant  fraction  of  the  cryopump  capacity,  as  seen  from  figure  A2.2.5. 
Thus,  the  cycling  also  avoids  the  large  gas  load  that  would  accumulate  during  many  days  of  growth  and 
which  could  easily  exceed  the  capacity  of  the  He  cyropumps,  leading  to  a  catastrophic  pressure  increase 
during  warm  up  of  the  wall  shrouds.  Also,  incomplete  dimerization  of  P4  to  P2  in  the  gas  cracker  can  lead 
to  a  build  up  of  white  P4  on  the  shrouds,  which  can  produce  an  increasing  background  pressure  with  each 
growth  and  lead  to  unwanted  P  incorporation  in  arsenide  compounds.  Frequent  cycling  of  the  shrouds 
avoids  this  problem. 

For  maintenance,  the  growth  chamber  is  opened  to  the  atmosphere  about  once  in  every  six  months  of  use. 
Before  opening,  the  MBE  chamber  is  baked  for  4—5  days  to  a  temperature  of  180— 200°C  to  stabilize  the 
phosphorus  coated  walls,  reducing  the  chance  of  a  white-  phosphorus  fire  and  minimizing  contamination 
fi-om  reaction  of  the  phosphorus  with  the  air.  Upon  opening,  we  usually  detect  toxic  gas  (e.g.  100  parts 
per  billion)  at  open  ports  on  the  growth  chamber;  thus,  self-contained  breathing  apparatus  must  be  worn 
during  servicing.  The  toxic  gas  levels  at  the  open  ports  usually  drop  to  below  detectable  concentrations 
within  about  24  hours  if  room  air  is  continually  drawn  through  the  chamber  and  exhausted.  Maintenance 
includes  cleaning  of  the  RHEED  gun  and  screen,  view  ports  and  effusion  cell  shutters,  and  replacement 
of  filaments  in  the  ion  gauges,  RHEED,  and  QMA.  We  have  found  that  most  of  the  phosphorus  deposited 
in  the  growth  chamber  is  on  the  surfaces  of  the  wall  shrouds  and  is  mixed  with  arsenic  deposits.  Little 
phosphorus  is  found  in  the  He  cryopumps  and  none  has  been  found  in  the  mechanical  pump  used  for 
cryopump  regeneration.  After  reloading  the  effusion  cells,  the  chamber  is  evacuated  and  again  baked 
for  4-5  days  with  the  pressure  held  constant  at  about  3  x  10"*  torr  until  the  wall  temperatures  reach 
180-200°C.  Upon  cooling,  the  base  pressure  is  in  the  10"'°  torr  range  and  the  system  is  carefully  leak 
checked  before  introducing  P2  or  As2.  After  re-establishing  the  calibration  of  the  group  V  and  dopant 
fluxes,  the  residual  free  carrier  concentration  is  measured  in  unintentionally  doped  material  to  check  the 
‘purity’  of  the  GSMBE  system.  We  have  used  this  up-to-air  procedure  for  more  than  15  openings  of  the 
growth  chamber,  and  we  find  no  increase  in  the  residual  earner  concentrations  in  either  GaAs  (p-type, 
p  <2x  10'*  cm"*)  or  InP  (n-type,  n  <  6  x  lO'*  cm"*).  Also,  we  have  found  no  degradation  in  the 
photoluminescence  intensity  of  InGaAsP  layers  after  opening  the  growth  chamber.  Thus,  with  the  proper 
operating  and  maintenance  procedures  the  quality  of  the  III-V  material  grown  in  a  phosphorus-based 
GSMBE  system  does  not  deteriorate  with  system  use.  Furthermore,  the  properties  of  HI— V  layers  are  as 
reproducible  as  those  grown  in  a  arsenic-based  solid-source  MBE  system. 


Safety 

Arsine  and  phosphine  are  extremely  toxic  and  great  care  must  be  exercised  in  the  design  and  operation 
of  a  GSMBE  laboratory  to  ensure  it  is  a  safe  working  environment  for  the  laboratory  personnel  and  that 
the  GSMBE  facility  does  not  present  a  major  health  hazard  to  the  community  [21,42,43].  Figure  A2.2.12 
summarizes  the  critical  concentrations  of  ASH3  and  PH3  relevant  to  health  hazards;  note  that  even  a  30 
min  exposure  to  only  30  parts  per  million  (ppm)  of  AsHs  can  be  lethal.  The  OSHA  Permissible  Exposure 
Limit  (PEL)  and  the  Total  Limit  Value  (ILV),  defined  as  the  maximum  concentration  that  a  worker  can  be 
exposed  to  in  an  8  hour  period  in  a  40  hour  work  week,  is  currently  50  parts  per  billion  (ppb)  for  ASH3 
and  300  ppb  for  PH3  [44].  The  gases  are  stored  in  pressurized  bottles  and  even  a  very  small  leak  can 
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result  in  concentrations  that  exceed  TLV  levels.  Thus,  to  provide  a  safe  environment,  the  laboratory  should 
be  designed  so  that  no  gas  is  allowed  to  escape  and  all  the  AsHs  and  PH3  that  is  consumed  is  converted  to 
a  safe  form  for  proper  hazardous  waste  disposal.  A  sufficient  volume  of  air  flow  must  be  provided  in  the 
laboratory  areas  for  ventilation  and  in  the  room  air  exhaust  system  in  order  to  dilute  the  toxic  gases  to  safe 
levels  in  case  of  a  leak  in  the  gas  handling  system.  The  hydrides  are  pyrophoric  [21]  and  toxic  PO;c  oxides 
are  released  during  a  PH3  fire;  thus,  provision  for  controlling  toxic  fires  must  also  be  made.  Hydride  gas 
detectors  [45]  must  be  used  to  continuously  monitor  the  laboratory  areas  and  gas  storage  cabinets  and  must 
be  connected  to  an  independent  alarm  system.  Self-contained  breathing  apparatus  must  be  used  during 
system  maintenance  and  be  readily  available  for  use  in  an  emergency.  Most  importantly,  all  personnel 
must  be  well  informed,  well  trained  in  the  proper  operation  of  the  equipment,  and  well  exercised  in  the 
emergency  procedures  in  case  of  an  accident  [46]. 
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Figure  A2.2.12.  Critical  concentrations  of  the  toxic  gases 
arsine  AsHa  and  phosphine  PH3. 


In  our  laboratory  the  hydride  gas  bottles  and  associated  gas  purge  system,  the  gas  handling  system  of  figure 
A2.2.6,  and  the  gas  scrubber  are  housed  in  a  single  vented  steel  cabinet.  The  cabinet  immediately  adjoins 
the  UHV  MBE  system  and  the  hydride  gas  lines  outside  of  the  cabinet  are  under  vacuum.  A  minimum 
amount  of  gas  is  stored  on  site  (0.5— 1.0  lbs  of  ASH3  and  1—2  lbs  of  PH3)  in  the  same  cabinet.  The  hydride 
flow  rates  are  low  (typically  <10  seem)  during  growth  and  thus  each  bottle  valve  is  left  closed  if  there 
is  no  limited-flow  aperture  supplied  with  the  gas  bottle,  minimizing  the  loss  of  toxic  gas  in  case  of  a  leak 
downstream  of  the  gas  bottles.  Two  toxic  gas  monitors  are  permanentiy  connected  to  an  integrated  alarm 
system:  a  colorimetric  paper  tape  monitor  (MDA  Scientific,  Inc.,  model  7100)  continuously  checks  the 
laboratory  area  and  a  solid-state  sensor  (1ST  Corp,  model  AGS  100)  monitors  the  exhaust  air  from  the  gas 
cabinet.  The  alarm  system  also  provides  an  indication  of  loss  of  the  main  electrical  power,  loss  of  vacuum 
in  the  MBE  system,  insufficient  flow  of  exhaust  air  from  the  room  and  gas  cabinet,  fire  in  the  gas  cabinet, 
insufficient  nitrogen  gas  pressure  for  purging  of  the  toxic  gases,  and  excessive  flow  of  the  toxic  gases 
from  their  storage  bottles.  The  alarm  system  is  connected  to  an  automatic  telephone  dialing  system  which 
can  alert  key  personnel  on  a  24  hour  basis.  The  functioning  of  the  alarm  system  is  tested  periodically. 
The  MDA  7100  monitor  and  the  microcomputer  for  control  of  the  MBE  system  receive  ac  power  from  an 
uninterruptible  power  supply.  A  separate  battery-operated  toxic  gas  monitor  (MDA  model  TLD-1)  is  used 
during  gas  bottle  changes  and  during  maintenance  of  the  growth  chamber. 

All  procedures  for  using  the  hydride  gases,  whether  for  epitaxial  growth  or  for  changing  the  gas  bottles, 
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are  carefully  documented  in  a  step-by-step  procedure.  The  same  procedure  is  used  by  each  crystal  grower 
at  the  start  of  each  day  of  growth  to  bring  the  hydride  gases  into  the  gas  handling  system.  In  normal 
operation,  all  AsHs  and  PH3  is  purged  from  the  gas  handling  system  and  bottle  pressure  regulators  at  the 
end  of  each  day  of  growth,  leaving  the  only  AsHs  and  PH3  present  in  the  laboratory  contained  in  the 
storage  bottles.  The  amount  of  gas  withdrawn  from  the  bottles  is  noted  and  used  to  predict  when  the  gas 
scrubber  will  need  changing  and  when  to  order  a  replacement  ASH3  or  PH3  bottle.  The  toxic  gas  botdes 
are  usually  changed  every  4-6  months. 

The  most  unpredictable  safety  hazard  of  GSMBE  is  a  phosphorus  fire.  This  is  most  likely  to  occur  during 
cleaning  of  the  growth  chamber  or  the  vacuum  pumps  connected  to  the  growth  chamber.  As  discussed 
above,  white  phosphorus  can  ignite  in  air,  and  the  resulting  fire  produces  a  dense  white  smoke  which 
contains  toxic  POj;  oxides.  Thus,  each  worker  must  wear  self-contained  breathing  apparatus  while  the 
growth  chamber  is  open  to  air.  A  fire  may  start  spontaneously,  but  more  likely  it  will  be  initiated  by 
mechanical  scraping  on  a  phosphorus-coated  metal  surface.  The  chance  of  a  phosphorus  fire  can  be 
minimized  by  using  low  P2  fluxes  during  growth,  baking  the  growth  chamber  prior  to  opening  to  air,  and 
flushing  the  chamber  with  dry  cold  N2  gas  during  repairs.  Liquid  nitrogen  is  a  very  effective  and  clean  fire 
extinguisher.  Large  phosphorus  fluxes  can  lead  to  substantial  phosphorus  build-up  in  mechanical  vacuum 
pumps,  with  fires  occurring  during  pump  maintenance  [44].  Furthermore,  the  oxides  in  the  smoke  can 
react  with  the  water  in  the  air  to  form  phosphoric  acid,  which  in  turn  is  highly  corrosive.  The  corrosive 
smoke  can  result  in  serious  damage  to  electronic  equipment  and  adequate  air  exhaust  in  the  laboratory 
must  be  provided  to  limit  smoke  damage.  An  ‘elephant’s  trunk’  connected  to  an  exhaust  duct  can  be  very 
effective  in  removing  the  smoke  (and  toxic  gases)  when  placed  in  the  immediate  vicinity  of  the  working 
area  or  loosely  connected  to  an  open  port  on  the  growth  chamber  during  repairs. 

As  shown  in  figure  A2.2.6,  a  gas  scrubber  is  used  to  convert  vented  ASH3  and  PH3  from  the  gas  handling 
system  and  any  uncracked  ASH3  and  PH3  from  the  growth  chamber  into  a  stable  form  for  safe  disposal.  For 
the  small  amount  of  hydrides  used  in  GSMBE,  we  have  found  a  dry  filter  containing  potassium  permanganate 
KMn04  to  be  a  simple,  cost-effective  means  of  neutralizing  the  toxic  gases.  KMn04  is  a  dry  solid 
that  reacts  with  the  hydrides  to  form  non-volatile  oxides.  The  scrubber  consists  of  100  lbs  of  Purafil 
Chemisorbant  (KMn04  in  an  alumina  binder)  placed  in  a  15  gallon,  sealable  polyethylene  drum  which 
also  serves  as  the  final  shipping  container.  The  lid  of  the  drum  is  fitted  with  two  pipes:  an  inlet  pipe  for 
introduction  of  the  effluent  from  the  GSMBE  system  at  the  bottom  of  the  Purafil  and  an  outlet  pipe  near 
the  top  of  the  Purafil  Chemisorbant  for  connection  to  the  exhaust.  When  the  scrubber  is  full,  the  pipes  are 
disconnected  and  the  drum  is  sealed  for  shipment;  no  handing  of  the  contaminated  Purafil  is  necessary. 
We  have  found  the  capacity  of  the  Purafil  scrubber  to  be  very  reproducible,  with  about  0.5  lbs  of  hydride 
adsorbed  for  every  100  lbs  of  Purafil.  We  change  the  scrubber  when  the  toxic  gas  monitor  first  indicates 
that  hydride  is  present  in  the  exhaust  from  the  scrubber,  which  occurs  at  intervals  of  6-9  months.  The 
scrubber  material  is  then  disposed  of  in  a  manner  consistent  with  governmental  regulations. 

Although  only  relatively  small  quantities  of  ASH3  and  PH3  are  needed  for  GSMBE,  their  extreme  toxicity 
presents  serious  safety  issues.  The  possibility  of  more  stringent  environmental  regulations  could  lead  to 
elimination  of  these  gases  for  commercial  purposes  in  the  future.  Several,  less  toxic,  direct  replacements 
for  the  hydrides  have  appeared  recently  for  use  in  MOMBE  and  MOCVD,  and  methods  that  reduce  the 
hazards  associated  with  on-site  storage  of  the  AsHs  and  PH3  under  pressure  have  been  developed  in  recent 
years.  The  compounds  tertiarybutylarsine  (TEA)  and  tertiarybutylphosphine  (TBP)  are  now  available 
as  substitutes  for  arsine  and  phosphine,  respectively  [47,48].  Both  TBA  and  TBP  are  liquids  at  room 
temperature  and  exhibit  vapour  pressures  that  make  delivery  with  a  earner  gas  feasible  in  both  MOMBE 
and  atmospheric  MOCVD.  Although  TBA  and  TBP  are  still  quite  toxic,  the  spillage  of  a  liquid  presents 
substantially  fewer  hazards  than  the  release  of  gaseous  ASH3  and  PH3  from  high  pressure  gas  cylinders. 
The  major  disadvantage  of  the  organic  materials  TBA  and  TBP  is  carbon  contamination  of  the  epitaxial 
layers  and  their  high  cost. 

One  method  to  avoid  the  hazards  during  the  transportation  and  storage  of  ASH3  and  PH3  in  pressurized 
cylinders  is  to  use  on-site  generation  of  the  hydride  gas  from  less  hazardous  materials.  Buckley  et  al 
[49]  have  shown  that  high  purity  arsine  can  be  electrochemically  generated  at  the  point  of  use  using  an 
arsenic  cathode  in  an  aqueous  electrolyte.  The  solid  arsenic  used  in  the  electrolytic  cell  is  considerably 
less  hazardous  than  ASH3  in  compressed  gas  cylinders. 

_ _ _ Z2 - 


Handbook  of  Thin  Film  Process  Technology 


0  1995  lOP  Publishing  Ltd 


Gas-Source  MBE 


Finally,  an  alternative  method  for  safer  transportation  and  storage  of  the  hydrides  is  to  use  adsorption  of 
the  hydrides  in  a  molecular  seive  at  atmospheric  pressure  [50],  The  molecular  sieve  consists  of  a  synthetic 
zeolite  in  a  sealed  container  which  is  baked  prior  to  absorbing  either  AsHs  or  PH3  upon  cooling  to  room 
temperature.  The  pressure  at  room  temperature  in  the  container  is  about  one  atmosphere  and  can  be  safely 
transported  to  the  point  of  use.  Upon  heating  to  about  200®C,  the  ASH3  or  PH3  is  released;  the  recovery 
can  be  almost  complete  and  thus  little  ASH3  or  PH3  is  lost.  The  zeolite  is  reusable  and  adequate  volumes 
of  gas  can  be  easily  stored  [50],  making  the  molecular  sieve  an  attractive  method  of  handing  AsHs  and 
PH3.  Arsine  absorbed  on  zeolite  has  been  used  recently  to  grow  GaAs  by  GSMBE  and  MOMBE  [51]. 


Growth  Procedure 

We  have  used  GSMBE  to  grow  a  variety  of  high  quality  ni-V  films  on  GaAs,  GaP,  InP,  and  Si  substrates. 
Preparation  of  the  ni-V  substrates  followed  the  same  procedures  as  used  in  conventional  solid-source  MBE 
[52].  Both  GaAs  and  InP  substrates  with  *epi  ready’  and  chemically  etched  protective  oxide  layers  have 
been  successfully  used.  Each  substrate  is  mounted  with  indium  solder  on  a  molybdenum  block,  or  in  the 
case  of  a  full  3-inch  diameter  wafer,  mounted  in  a  molybdenum  ring  without  In  solder.  The  protective 
oxide  is  then  desorbed  in  the  growth  chamber  by  heating  the  block  or  ring  under  the  appropriate  group  V 
molecular  beam.  Epitaxial  growth  then  proceeds  in  a  fashion  identical  to  solid-source  MBE. 

In  the  case  of  growth  on  Si  substrates,  we  have  found  the  use  of  phosphorus  in  the  growth  chamber 
can  influence  the  procedure  for  thermal  cleaning  of  Si  substrates.  (100)  Si  substrates  were  etched  in 
5:1  H2S04:H202,  then  transferred  to  a  N2  atmosphere  and  briefly  etched  in  a  HF:H20  solution  before 
loading  into  the  MBE  system.  The  substrates  could  then  be  cleaned  in  the  growth  chamber  at  800°C 
for  20  min,  using  RHEED  to  monitor  the  surface  cleanliness.  The  process  was  very  reproducible  and  a 
number  of  InP-on-Si  heterostructures  were  successfully  grown  [53,54].  However,  after  about  one  year  of 
growth  in  the  same  MBE  chamber,  it  became  increasing  difficult  to  obtain  a  RHEED  pattern  typical  of  a 
clean  well  ordered  Si  surface,  and  the  structural  quality  of  the  InP  epitaxial  films  deteriorated  noticeably. 
Measurements  indicated  that  a  small  amount  (i.e.  a  few  atomic  per  cent)  of  P  was  being  deposited  on 
the  hot  wafer  surface  from  the  walls  of  the  chamber.  The  P  apparently  reacted  with  any  oxide  on  the  Si 
surface  to  form  a  thin  phosphosilicate  layer  whose  vapour  pressure  was  high  enough  to  prevent  adequate 
sublimation  to  occur  at  8CX)®C.  By  exposing  the  surface  to  a  Ga  molecular  beam  at  a  low  flux  (about  0.1 
^m  hr”^)  for  30  min  at  800°C,  a  clean  Si  surface  could  be  produced  as  evidenced  by  a  well-defined  (2  x  2) 
RHEED  reconstruction  pattern.  Presumably  the  Ga  decomposes  the  phosphosilicate  layer,  producing  volatile 
reaction  products  at  800®C  [55].  The  Ga  beam  cleaning  technique  has  consistently  produced  atomically 
clean  Si  surfaces  in  our  phosphorus-based  MBE  system. 

One  advantage  of  GSMBE  over  sold-source  MBE  is  the  capability  to  rapidly  alter  the  group  V  flux  and 
produce  abrupt  changes  in  arsenic  and  phosphorus  composition.  This  capability  has  been  exploited  in  the 
growth  of  several  quantum  well  heterostructures  consisting  of  alternating  arsenide  and  phosphide  layers. 
Examples  include  QWs  of  GaAs  with  barriers  of  InGaP  and  QWs  of  InGaAs  with  barriers  of  InP.  The 
growth  is  usually  paused  at  each  heterointerface  to  increase  the  smoothness  of  the  growth  front  and  to 
obtain  atomically  abrupt  interfaces.  Because  of  the  high  vapour  pressure  and  low  sticking  coefficients 
of  the  group  V  elements,  the  switching  of  the  group  V  molecular  beams  during  the  growth  pause  may 
produce  an  interfacial  region  with  a  composition  that  differs  significantly  from  that  of  the  surrounding  well 
and  barrier  layers.  A  grading  in  bandgap  and  localized  strain  at  the  interface  can  result  and  can  adversely 
affect  device  operation.  The  method  of  switching  the  AsHs  and  PH3  gas  flows  and  the  pumping  speed  in 
the  GSMBE  chamber  can  greatly  influence  the  interfacial  compositional  abruptness. 

A  typical  example  of  the  sequencing  of  the  molecular  beams  during  the  growth  of  GaAs/InGaP  QWs  is 
shown  in  figure  A2.2.13,  where  the  Ino.48Gao.52P  alloy  composition  is  chosen  to  be  lattice  matched  to 
the  GaAs  substrate.  The  growth  pause  at  the  InGaP-to-GaAs  interface  is  of  duration  t\  + 12  and  at  the 
GaAs-to-InGaP  interface  of  duration  +  ^2-  E)uring  the  interval  t\  the  InGaP  growth  surface  is  allowed 
to  anneal  in  a  protective  P2  overpressure.  The  gases  are  then  switched  and  growth  recommences  after 
interval  tz-  A  similar  sequence  is  preformed  at  the  GaAs-to-InGaP  interface  but  using  different  intervals 
and  ^2-  A  systematic  study  of  the  effect  of  varying  t\,  t2,  t[  and  on  the  structural  and  compositional 
abruptness  of  each  interface  was  carried  out  using  double-crystal  x-ray  (DCXR)  diffraction,  transmission 
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electron  microscopy  (TEM),  photoluminescence  (PL),  and  RHEED  [33,34].  It  was  found  that  t\  =  10  s  and 
/]'  =  30  s  were  long  enough  at  the  growth  temperature  to  provide  atomically  abrupt  interfaces.  The  InGaP 
surface  required  very  little  annealing,  but  the  smoothness  of  GaAs  surface  improved  with  longer  times  for 
t[.  The  time  tz  (or  must  be  chosen  short  enough  to  minimize  intermixing  on  the  group  V  sublattice 
resulting  from  incorporation  of  As  (or  P),  but  long  enough  to  allow  the  residual  P2  (or  AS2)  to  be  pumped 
from  the  grovrth  chamber  and  to  allow  the  As2  (or  P2)  flux  to  stabilize  before  starting  the  growth  of  GaAs 
(or  InGaP).  With  the  pumping  system  of  figure  A2.2.4  and  the  gas  handling  system  of  figure  A2.2.6,  we 
found  r2  =  /j  =  6s  or  less  to  be  suitable.  Using  this  gas  switching  sequence  GaAs/InGaP  QWs  have  been 
grown  in  a  reproducible  manner  with  interfacial  regions  of  1-2  ml  in  thickness  and  with  minimal  As— P 
intermixing  [34]. 


Similar  considerations  have  been  applied  to  the  growth  of  InGaAs/InP  QWs  by  GSMBE,  where  the 
Ino.53Gao.47As  alloy  composition  is  lattice-matched  to  InP.  Panish  and  co-workers  [32,56]  had  previously 
used  a  sequence  similar  to  that  of  figure  A2.2.13  where  the  group  V  beams  are  switched  under  conditions 
such  that  a  group  V  overpressure  is  always  present  during  each  growth  pause.  Pausing  for  5—15  s  at 
each  interface  produced  InGaAs-to-InP  and  InP-to-InGaAs  interfaces  which  exhibited  equal  and  opposite 
strains.  The  strained  interfacial  regions  occur  as  result  of  asymmetric  ordering  inherent  in  the  atomic 
arrangement  of  perfectly  abrupt  interfaces.  Thus  the  interfaces  obtained  were  atomically  abrupt  and  ideal¬ 
like  in  composition.  Atomically  abrupt  interfaces  in  InGaAs/InP  heterostructures  have  also  been  obtained 
by  Moy  etal  [35]  but  using  a  different  gas-handling  and  pumping  arrangement  and  a  different  gas  switching 
sequence.  Here  growth  is  paused  for  17-20  s,  but  P2  and  As2  is  rapidly  pumped  from  the  growth  chamber 
and  no  group  V  overpressure  is  provided,  except  for  a  10  s  interval  prior  to  recommencing  growth  of 
InGaAs.  Apparently  the  lack  of  protective  group  V  overpressure  at  the  InGaAs-to-InP  interface  causes  no 
adverse  effects  as  long  as  the  growth  pause  is  kept  short,  but  the  As2  beam  prior  to  growth  of  InGaAs  is 
needed  to  obtain  an  abrupt  InP-to-InGaAs  interface. 


Another  important  issue  is  cross  contamination  of  the  group  V  elements  during  GSMBE  growth.  For 
example,  P  can  be  incorporated  in  GaAs  as  result  of  incomplete  removal  of  P  from  the  growth  chamber  by 
the  vacuum  pumps  and  cyropanels  or  from  the  evaporation  of  P  from  the  walls  of  the  growth  chamber  onto 
the  growth  surface.  The  degree  of  contamination  is  strongly  dependent  on  the  incorporation  coefficient 
of  the  unwanted  group  V  element;  the  design  of  the  growth  chamber  and  associated  vacuum  pumps,  gas 
delivery  system,  and  gas  cracker;  and  the  recent  growth  history  of  the  MBE  system.  Hgure  A2.2.14  is  a 
DCXR  spectrum  of  a  GaAs  epitaxial  film  of  1.86  /rm  in  thickness  grown  on  a  GaAs  substrate  one  day  after 
InP  was  grown  in  the  same  vacuum  chamber  and  with  the  same  gas  cracker.  Two  peaks  are  can  be  resolved 
in  the  spectrum:  the  peak  on  the  left  is  from  the  GaAs  substrate  and  the  peak  on  the  right  is  from  an 
epitaxial  film  of  composition  GaAso.999Po.001  •  We  have  found  As  contanunation  in  InP  is  somewhat  higher 
than  P  in  GaAs  probably  because  the  incorporation  coefficient  for  As  is  larger  than  for  P  in  InGaAsP  [36]; 
the  composition  is  typically  InPo.998Aso.oo2'  Fot  tnost  applications,  this  level  of  group  V  contamination  is 
insignificant 


The  degree  of  cross  contamination  is  undoubtedly  higher  in  layers  of  heterojunction  samples  where  the 
group  V  species  from  the  previous  layer  has  not  been  completely  pumped  out  of  the  chamber  before 
growth  of  the  next  layer  begins.  However,  with  proper  sequencing  of  the  AsHs  and  PH3  gas  flows, 
rapid  gas  switching,  and  adequate  pumping  in  the  growth  chamber,  cross  contamination  can  be  reduced 
to  acceptable  levels.  Jackson  and  coworkers  [31]  have  examined  the  degree  of  cross  contamination  in 
InGaAs/InP  heterostructures  and  their  results  are  seen  in  the  SIMS  depth  profiles  in  figure  A2.2.10.  No 
growth  pause  was  employed  so  that  the  speed  of  response  of  the  hydride  gas  delivery  system  and  gas 
cracker  could  be  clearly  evaluated.  For  the  sample  of  figure  A2.2. 10(a),  a  single  low  conductance  cracker 
and  a  slow  gas  switching  system  was  used.  A  long  P  tail  penetrates  deep  in  the  InGaAs  layer  and  the 
residual  P  is  substantial  (>  0.01%).  In  figure  A2.2.10(b),  separate  high  conductance  crackers  were  used 
for  AsHj  and  PH3  along  with  a  fast  gas  switching  system.  The  P  tail  is  greatly  reduced  and  the  residual 
P  in  the  InGaAs  is  below  the  SIMS  detection  limit  (<  0.0003%).  Thus,  cross  contamination  of  the  group 
V  elements  in  heterostructures  can  be  reduced  to  insignificant  levels. 
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Figure  A2J2A3.  Sequence  of  molecular  beams  during  the  growth 
of  InGaP/GaAs  heterostructures.  Growth  is  interrupted  by  closing 
the  group  III  shutters  at  the  InGaP-to-GaAs  interface  for  the 
period  ti  +  ti  and  at  the  GaAs-to-InGaP  interface  for  the  period 
Annealing  of  the  growth  front  occurs  during  the  intervals 
ti  and  r[. 
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Figure  A2.2.14.  High  resolution  double-crystal  x-ray 
rocking  curve  of  a  GaAs  film  grown  on  a  GaAs 
substrate  using  a  gsmbe  growth  chamber  with 
residual  phosphorus  vapour.  The  composition  of  the 
film  is  determined  to  be  GaAs0.999P0.001  from  the 
separation  of  the  film  peak  and  the  substrate  peak. 
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A2^2  SUMMARY 

We  have  previously  stated  that  the  primary  advantage  of  GSMBE  is  the  capability  to  grow  a  wide  variety  of 
ni-V  materials  in  one  growth  chamber  with  relative  ease.  This  capability  is  illustrated  in  figure  A2.2.15 
where  PL  spectra  are  shown  for  a  variety  of  III-V  samples,  all  grown  in  the  same  GSMBE  chamber.  Each 
sample  consisted  of  an  epitaxial  layer  of  1-2  iim  in  thickness  and  lattice  matched  to  the  appropriate 
substrate  material.  The  materials  range  from  the  InGaAs  with  emission  in  the  infirared  region  of  the 
spectrum  to  alloys  of  InGaAlP  with  emission  in  the  visible  region.  Portions  of  the  spectrum  not  covered  in 
figure  A2.2.15  have  been  addressed  with  lattice  mismatched  (i.e.  strained)  alloys  and  QW  heterostructures 
of  the  same  materials.  Growth  of  any  of  the  materials  shown  can  take  place  in  any  order  and  in  the  same 
day,  without  compromising  film  quality.  Thus,  GSMBE  allows  great  flexibility  in  sample  design  as  well  as 
scheduling  of  growth  runs. 

Photoluminescence  at  300 K  of 

GSMBE  ill'V  Semiconductors 


Energy.  (eV) 

Figure  A2J2.15.  Photoluminescence  spectra  of  different 
III-V  epitaxial  samples  grown  in  the  same  gsmbe 
growth  chamber.  Note  the  wide  range  of  energies  for  the 
emitted  light.  The  InGaAs,  InGaAsP,  and  InP  samples 
were  lattice  matched  to  InP  substrates,  and  the  InGaP, 

InGaAlP,  and  GaAs  samples  were  lattice  matched  to 
GaAs  substrates. 

When  combined  with  the  precise  control  of  layer  thicknesses  on  the  atomic  scale  of  conventional  MBE, 
GSMBE  is  a  powerful  tool  for  fabricating  a  wide  range  of  ni-V  heteroepitaxial  devices.  This  feature  is 
exemplified  by  the  high  resolution  TEM  lattice  image  in  figure  A2.2.16  of  the  cross  section  of  a  GaAs/InGaP 
QW  sample  grown  by  gsmbe  [57].  Here  a  film  of  GaAs  only  two  monolayers  (0.57  nm)  in  thickness  is 
bounded  on  both  sides  with  thick  layers  of  lattice-matched  InGaP.  Heterostructures  with  similar  interface 
abruptness  and  layer  thickness  control  have  been  routinely  obtained  in  other  III— V  material  combinations 
using  GSMBE  technology 

In  summary,  we  have  provided  the  technological  details  for  the  design  and  operation  of  a  GSMBE  system 
when  using  the  hydride  gases  as  source  materials  for  arsenic  and  phosphorus.  With  the  procedures  outlined, 
reproducible  epitaxial  growth  of  complex  heterostructures  can  be  achieved  in  a  safe  and  reliable  manner. 
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Figure  A22.16.  tem  lattice  image  of  a  [010]  cross  section  of  a  InGaP/GaAs  quantum  well  grown  by  gsmbe.  The 
GaAs  well  layer  is  the  thin  darker  band  across  the  centre  of  the  micrograph.  Tlie  {200}  atomic  planes  are  resolved 
and  the  GaAs  QW  is  two  monolayers  thick  with  atomically  abrupt  interfaces,  tem  photo  by  N  Otsuka  and  D  Li  at 
Purdue  University, 
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